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Abstract 

Precision Rosenbluth Measurement of the Proton Elastic 
Electromagnetic Form Factors and Their Ratio at Q 2 = 2.64, 

3.20, and 4.10 GeV 2 
Issam A. Qattan 

Due to the inconsistency in the results of the M ^P gp ratio of the proton, as extracted 
from the Rosenbluth and recoil polarization techniques, high precision measurements 
of the e-p elastic scattering cross sections were made at Q 2 = 2.64, 3.20, and 4.10 
GeV 2 . Protons were detected, in contrast to previous measurements where the scat- 
tered electrons were detected, which dramatically decreased e dependent systematic 
uncertainties and corrections. A single spectrometer measured the scattered protons 
of interest while simultaneous measurements at Q 2 = 0.5 GeV 2 were carried out us- 
ing another spectrometer which served as a luminosity monitor in order to remove 
any uncertainties due to beam charge and target density fluctuations. The absolute 
uncertainty in the measured cross sections is ~ 3% for both spectrometers and with 
relative uncertainties, random and slope, below 1% for the higher Q 2 protons, and 
below 1% random and 6% slope for the monitor spectrometer. The extracted electric 
and magnetic form factors were determined to 4%-7% for Ge p and 1.5% for G Mp . 
The ratio was determined to 4%-7% and showed « 1.0. The results of 

this work are in agreement with the previous Rosenbluth data and inconsistent with 
high-Q 2 recoil polarization results, implying a systematic difference between the two 
techniques. 
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Chapter 1 Introduction 



1.1 Overview 

The modern concept of the atom as an object composed of a dense nucleus surrounded 
by an electron cloud came to life as a result of Rutherford's experiments in 1911. As 
time went on, it became clear that the nucleus itself is composed of even smaller 
objects or particles which we call nucleons. The fact that the proton's magnetic 
moment is 2.79/itv (/^n = 2lTc)> ra ther than the 0.5/itv expected for a point-like spin 
1/2 charged particle demonstrates that the proton posses a structure. Similarly, the 
neutron's magnetic moment of -1.91/itv is very different from the 0.0/ztv that a point 
neutral particle would have. Thus, nucleons, protons and neutrons, which were once 
thought to be the fundamental building blocks of matter in fact are complex particles 
that turn out to be composed into even smaller objects which we call quarks and 
gluons. The whole picture boils down to electrons (leptons), quarks, and gluons as 
the fundamental types of building blocks. 

Several challenging questions and issues arose when physicists tried to understand 
nucleons. For example, what are these particles made off, how do they interact, and, 
most importantly, how do they bind together? These are fundamental and compli- 
cated questions that must be addressed, answered, and resolved by the nuclear and 
particle physics communities. The strong force is the reason why protons and neu- 
trons bind together to form the nuclei. It is the theory of Quantum Chromodynamics 
(QCD) that describes the strong interaction between the quarks and gluons, which 
in turn make up the protons and neutrons. 

The strong coupling constant, a s , as defined by the theory of Quantum Chromo- 

1 



2 

dynamics P, is a measure of the strength of the strong interaction between quarks 
and gluons and is a function of the four momentum transfer squared Q 2 : 

(33-2^)1^)' 

where Nf is the number of flavors active for QCD renormalization scale Aqcd 0EI 
When Q 2 is large (short distance scale), equation lj) tells us that a s is small. 
Therefore, at large momentum transfer, the quarks within a hadron (neutron or 
proton) are weakly interacting and behave almost as if they are free particles. On 
the other hand, when Q 2 is small (large distance scale), a s is large and the quarks 
are strongly interacting particles and hence form hadronic matter. 



1.2 Electron Scattering 

In order to reveal the underlying structure of the nucleon, experimental techniques 
such as electron scattering have been developed. Electron scattering in particular has 
proven to be a powerful tool for studying the structure of the nucleus. The reason 
lies in the fact that the electron is a point-like particle and has no internal structure, 
making it a clean probe of the target nucleus. In this case, the information extracted 
such as the differential cross section reflects the structure of the target without any 
contribution from the projectile. The incident electron is scattered off a nuclear target 
(single particle in the target) by exchanging a virtual photon. The electron-photon 
scattering vertex is known and understood within the theory of QED. 

In electron scattering experiments, highly relativistic electrons are used. At low 
energy transfers, the virtual photon interacts with the entire nucleus. The nucleus 
stays intact and the electron scatters elastically or excites nuclear states. In this case, 
the virtual photon is dominantly interacting with the nuclear target by coupling to a 



3 

vector meson resonance states or qq jlj |3] . 

As the energy transfer increases and becomes larger than the nuclear binding 
energy, we enter the quasi-elastic region. The virtual photon becomes more sensitive 
to the internal structure of the nuclear target and the scattering process is described 
in terms of photon-nucleon coupling rather than photon-meson coupling. In this case, 
the target is viewed by the virtual photon as a set of quasi-free nucleons. The electron 
scatters elastically from the nucleon which in turn is ejected from the nucleus. 

By increasing Q 2 and the energy transfer further, we enter the resonance region. 
In this case, the virtual photon becomes sensitive to the internal structure of the 
nucleon. The quarks inside the nucleon absorb the virtual photons and form excited 
states which we refer to as nucleon-resonances. 

If we increase Q 2 and the energy transfer even further, we enter the deep inelastic 
region. The virtual photon probes smaller distance scales and has enough energy 
to resolve the constituents of the nucleon, or partons. Perturbative QCD models 
(pQCD) E] seem to give the correct description of the scattering process at high 

Q 2 . 

In electron scattering experiments, see Figure ITTT| the initial electron's energy and 
momentum k = (E, k) are known. Elastic electromagnetic scattering of the electron 
can be described, to lowest order in the electromagnetic coupling constant a, as the 
exchange of a single virtual photon 7* of momentum q and energy u. The virtual 
photon then interacts with the target. After the scattering, the final electron's energy 
and momentum k! = (E', k') are measured. This allows us to determine the energy 
uj = (E — E') and momentum q={k — k') of the virtual photon. Usually we refer to 
the energy of the virtual photon, u = (E — E') , as the energy transfer. 

We describe the scattering process using two main variables, the energy transfer 
to the virtual photon, u = E — E', and the square of the 4-momentum transfer, 
Q 2 = —q^ = —q 2 = \k — k'\ 2 — (E — E') 2 . We define the Bjorken variable x = ^fp^ 
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Figure 1.1: Feynman digram for electron-nucleus scattering. 

where M p is the mass of the nucleon. The Bjorken variable x varies between and 
1.0, where the x = 1.0 case corresponds to elastic scattering, and the x < 1.0 case 
corresponds to inelastic scattering. 

The 4-momentum transfer Q 2 = —q 2 is given in terms of the electron's kinematics 
(initial energy E, final energy E', and the scattering angle 9 e ) as: 

Q 2 = 4££'sin 2 - . (1.2) 

We refer to the scattering as an exclusive scattering if the final state is fully 
identified. In this case we can also express the scattering kinematics in terms of the 
knocked out or recoil target (nucleon) rather than the electron. If the initial and the 
final four-momentum of the nucleon are defined as P = (E p ,p) and P' = (E' p ,p / ), 

— * 

where E p and E' p are the initial and final energy of the nucleon and p and p' are 
the initial and final 3-momentum of the nucleon, we can write the the 4-momentum 
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transfer Q 2 as: 

Q 2 = -q 2 = -(P> - Pf = {{E' p - E p ) 2 - (P - p) 2 ) . (1.3) 

For the case where the final hadronic state contains a single nucleon only, the 
missing mass squared of the scattered hadron W 2 must equal the mass of the nucleon 
squared, M 2 , where we define the missing mass squared as: 

W 2 = M 2 + 2M p u - Q 2 . (1.4) 

Examining equation (jl.4)l . we see that if the final state contains only a single 
nucleon, then W 2 = M 2 and Q 2 = 2M p u yielding a value of x = 2 ® - = 1 for the 
Bjorken variable. So by choosing scattering kinematics where x — 1, we can isolate 
elastic scattering. 

If the nucleus is knocked into an excited state, there is some additional energy 
transfer, and x will decrease as the energy transfer increases. At somewhat higher 
energy transfer, where quasielastic scattering is the dominant process, the electron 
knocks a single nucleon out of the nucleus. This corresponds to scattering near x = 1. 
At higher energy transfers, corresponding to x < 1, the scattering is inelastic and the 
struck nucleon is either excited into a higher energy state (resonance scattering), or 
temporarily fragmented (deep inelastic scattering). At very high energy transfers, 
where deep inelastic scattering dominates, the electron is primarily interacting with 
a single quark via the virtual photon 7*. In this short wavelength virtual 

photon 7* is needed in order to resolve the quarks within the proton. 
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1.3 Elastic Electron-Proton Scattering 

The electron-photon vertex as shown in Figure 11.11 is well understood and described 
by QED. On the other hand, the photon-proton vertex is complicated and the detail 
of such interaction, 7* + P(p) — > P(p') , cannot be calculated from first principles. 
This is due to the fact that the proton is not a point-like particle but rather a particle 
with internal structure. We introduce two Q 2 -dependent functions that contain all 
the information about the internal structure of the proton. We refer to these two 
functions as the proton electromagnetic form factors which parameterize the internal 
structure of the proton. 

For the single-photon exchange diagram for the electron-proton elastic scattering 
shown in Figure ITTT| the Lorentz invariant transition matrix element is given by: 

M fi = (j^k,k')^J»(p,p')) , (1.5) 

where and J M are the electromagnetic currents for the electron and proton, respec- 
tively: We can express the electromagnetic currents as: 

jp = -eu^k'^^k) , (1.6) 



;i.7) 



where u(k), u(k'), v(p), and v(p') are the four-component Dirac spinors for the initial 
and final electron and proton, respectively, which appear in the plane-wave solutions 
for Dirac equation, and T M and 7^ are the Dirac 4x4 matrices. Therefore, the Lorentz 
invariant transition matrix element can be expressed as: 



iM 



fi 



iev(p')T' 1 (p, p')v{p) ieu(k')^iju{k) 



Because of the composite nature of the proton, T M rather than 7^ is used to 
describe the proton current since T M contains all the information about the internal 
electromagnetic structure of the proton. It is worth mentioning that if the proton 
were to be treated as a point-like particle, then T M — > 7 M as in the case with the 
electron. The proton current J M is a Lorentz four-vector that is Lorentz invariant 
and satisfies both parity and current conservations in electromagnetic interactions, 
i.e. <9J M = 0. With this in mind, the proton current can be written as: 

J" = v{p'){fM 2 )i, + ^F 2 (^)a^)v{p) > (1-9) 

where k p = 1.793/xjv is the proton anomalous magnetic moment and it is expressed in 
the unit of nuclear magneton (/zat), M p is the mass of the proton, a^ v = f [7^, 7"], and 
the factor has been included as a matter of convention. The functions Fi(q 2 ) and 
F 2 (q 2 ) are functions of q 2 and known as Dirac and Pauli form factors, respectively. 
The Dirac form factor, Fi(q 2 ), is used to describe the helicity-conserving scattering 
amplitude while the Pauli form factor, F 2 (q 2 ), describes the helicity-flip amplitude. 
In the limit that q 2 — ► 0, the structure functions Fi(q 2 = 0) = F 2 (q 2 = 0) = 1.0, and 
in this limit, the virtual photon becomes insensitive to the internal structure of the 
proton which is viewed as a point-like particle. 

The elastic differential cross section in the lab frame for the e(k) + P(p) — > 
e(k') + P'{p') reaction is given by: 



da 



l M /*l 2 (27r )W (k' + v-k- A d3k ' d3p ' 

4((ifc • p) 2 - mlM 2 ) V n(2nf2E\2nf2{M p + uj 



;i.io) 

By using the scattering amplitude Mfi expressed in terms of the electron and 
proton currents, integrating out the 5-function that imposes momentum conservation, 
and finally expressing the proton current in terms of the Dirac and Pauli form factors, 



the differential cross section for an unpolarized beam and target can be written as: 



da 



a 



E' 



dVtdE' 



4£ 2 sin 4 (f ) E 



F 2 + 



4Q 2 



4M^ W& 



Q 2 , „ oAs 

-^(F 1 + S F 2 ) 2 sin 2 (^) 



5(E - E' 



2M„ 



1.11) 



where we have averaged over initial spins and summed over final spins. It is un- 
derstood that the 5-function is used to assure elastic scattering, that is, at a given 
energy E and angle 6 e , the elastic differential cross section is a ^-function in E' at 
E' = E — 7^-. If we integrate over E' in equation (jl.lljl above and divide the 
numerator and denominator by cos 2 ^); we can write the elastic cross section as: 



da 

dn 



a,, 



2M„ 2 2 



Q 2 , „ oAs 

^(F 1 + S F 2 ) 2 tan 2 (^) 



;i.i2) 



where a ns in equation (jl. 12)1 above is known as the non-structure cross section and is 
given by: 



a„ 



c^cos 2 (f) E' 
4E 2 sin 4 (|)¥ 



;i.i3) 



where a is the fine structure constant. 

It is worth mentioning that a ns , equation (|1.13|) . is nothing but the famous Ruther- 
ford formula for the elastic electron-proton scattering modified to account for the 
proton's recoil and the spin-orbit coupling effects: 

1. The recoil effect: The term which is relativistic in nature, is due to the recoil 
of the proton. Although the protons are massive, their recoil effect can not be 
neglected at high momentum transfer squared Q 2 . We can write ^ in terms of 
the electron's kinematics as: 



1 

1 + Hsin 2 (!)' 



(1.14) 



where it should be clear that the term (^) — > 1.0 as Q 2 — > 0. 

2. The spin-orbit coupling: The electron is a spin-| particle and has a magnetic 
moment /i e which interacts with the magnetic field B p of the proton as it is 
felt by the electron in its own frame of reference. This is referred to as the 
spin-orbit coupling effect which results in the cos 2 (y) term or more precisely, 
1 — (3 2 sin 2 ( y), with j3 = (|) — > 1.0 for extremely relativistic electrons. 

We refer to the cross section without taking into account the recoil effect as the 
Mott cross section and it is given by: 

^ _(^) 2 (l-/3 2 sin 2 (|)) 
UQJmou 4Psin 4 (|) ' ( ' 5) 

where k is the initial momentum of the incident electron and Z is the atomic number. 
In the case of a high momentum electron scattered off a spin- 1/2 point-like proton 
(Z = 1), we recover equation ()1.13|) without accounting for the recoil effect i.e. — ► 
!)■ 

In the non-relativistic limit where (3 — > 0, E — > k, and ^ — > 1, equation (jl,13j) 
becomes: 



da\ a 2 



/da 

Rutherford Ak 2 SU1 4 ( y 



1.16) 



and this is the famous Rutherford cross section formula. 
To put things in more perspective we can write: 



fda\ 

\d~QJ Rutherfo 



2 



E'\ fda\ (E'\ 



i-^ra(iJ = ^-l¥J- (L17) 



which describes the proton as a spin-1/2 point-like particle without any internal 
structure. 



1.4 Form Factor Interpretations 
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As shown earlier the concept of the nucleus as a point-like particle with a static 
charge distribution g(r) = Ze5(r) is not the correct one. This is due to the fact that 
the nucleus has an internal substructure in the form of electric charge and magnetic 
moment distribution. This substructure modifies the non-structure scattering cross 
section as given by equation ()1.13j) to: 

where F(Q 2 ) is the form factor which accounts for the fact that the target nucleon 
possess an internal structure. It should be mentioned that F is a function of Q 2 alone 
in the case of an elastic scattering. 

It can be shown pQ that in the non-relativistic limit, the form factor F(Q 2 ) can 
be expressed in terms of the charge distribution g(f) of the target nucleus as: 

Fnr(Q 2 ) = J gtfe&W? , (1.19) 

where NR stands for Non-Relativistic. Equation ()1.19|) is the Fourier transform of 
the electric charge distribution of the target. In the case of the proton, F NR {Q 2 ) 
corresponds to the electric form factor, Ge p {Q 2 )- In the same way, if the target has 
an extended magnetic moment distribution, the magnetic form factor is the Fourier 
transform of the magnetic moment distribution of the target and is referred to as 
Gmp(Q 2 )- 

If we scatter off a point-like particle with static charge distribution g{f) = ZeS(r), 
equation ()1.19|) gives Fnr(Q 2 ) = Z, and in the case of an extended spatial charge 
distribution g(f) = po e_ ^"\ equation ()1.19j) gives the well-known dipole form factor 
Fnr{Q 2 ) = Gd(Q 2 ) = (1 + a 2 Q 2 )~ 2 where a is the scale of the proton radius. Ex- 
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per iment ally, it was observed that both the electric and magnetic form factors could 
be described to a good approximation by the dipole form factor Gd(Q 2 ) with a 2 = 
(O.TlGeV 2 )' 1 . Figures \l.2\ IO| and 11.41 show the world data on the electric and 
magnetic form factors. Figure 11.41 also shows ^ p q P d with the fits of Bosted jS] and 
Brash |U] to the proton magnetic form factor. The Bosted and the Brash fits will be 
discussed in detail in the next chapter. 




Figure 1.2: Ge p world data. The black solid line is the dipole form factor fit to data. 
Data from references [TOJ HH H21 CHI HH HS1 HEl HZJ - 

If we assume that the magnetic moment of the proton has the same spatial de- 
pendence as the charge distribution, then the electric and magnetic form factors (in 
the non-relativistic limit) are related by l^poff- = 1.0 and this is referred to as form 
factor scaling or: 

G Ep (Q 2 ) = ^lA^l . (i.20) 
Hp 

The electric and magnetic form factors seem to approximately follow the dipole 



12 



3 1.75 



1.5 
1.25 
1 

0.75 
0.5 
0.25 




* Andivahis et al. Price et al. 

• Walker et al. * Berger et al. 
Janssens et al. Bartel et al. 

A Hanson et al. 



10 



Q (GeV/c) 



12 

2 



Figure 1.3: Ge p world data normalized to the dipole form factor Go- Data from 
references HH H31 HH US IIEl UZ1 UH] 



p 1.15 

O 

i a 1.1 

"""""a 

S 1.05 
1 

0.95 
0.9 

0.85 
0.8 

0.75 
0.7 

0.65 



_p -p 

Kl • ir 

fib* IT _L I J_ \ 






















1 '' 




* 


Andivahis et al. 


J- - - 


• 


Walker et al. 


i- • - 




Bartel et al. 






Janssens et al. 






Price et al. 




▲ 


Hanson et al. 




* 


Berger et al. 




• 


Sill et al. 






Bosted Fit (Bosted et al.) 






Brash Fit (Brash et al.) 





10 



12 



(T (GeV/c) 



14 

2 



Figure 1.4: Gmp world data normalized to the dipole form factor Go- The solid and 
the dotted lines are the Bosted and Brash fits, respectively. Data from references 

nn n2i csi us nni nzi chi nni 



13 

form factor as shown in Figures 11.31 and 11.41 This indicates that both form factors 
have the same Q 2 dependence although Gm p seems to significantly deviate from dipole 
form factor at high Q 2 . 

For \Q\r « 1, equation ()1.19|) can be expanded in terms of the root-mean-square 
charge radius as [2D]: 

tt (r^\ 1 l^| 2 ^ ^ ^r.m.s \Q\ 4 < ?" 4 >r.m.s /1 oi\ 

*nr[Q ) = 1 g 1 1 ' ^ ' 

where the root-mean-square charge radius is given by: 

< r 2 > r . m . s = J g(f)r 2 d 3 r . (1.22) 

Finally, it is important to realize that the interpretation of the form factors (elec- 
tric and magnetic) as the Fourier transform of the electric charge and magnetic mo- 
ment distribution is only valid in the non-relativistic limit. 



1.5 Rosenbluth Separations Technique 

A linear combinations of the Dirac and Pauli form factors Fi(q 2 ) and F 2 (q 2 ) can be 
used to define the Sachs form factors 21J, Ge p and Gm p , the electric and magnetic 
form factors of the proton: 

G Ep (Q 2 ) = F X (Q 2 ) - k p tF 2 (Q 2 ) , (1.23) 



G Mp (Q 2 ) = F 1 (Q 2 ) + k p F 2 (Q 2 ) , (1.24) 



where r = jjjjy. In the limit Q — > 0, where the virtual photon becomes insensitive 
to the internal structure of the proton, equations (jl.23j) and (jl.24j) reduce to the 



14 



normalization conditions for the electric and magnetic form factors respectively: 



GeM = Fi(0) = 1 



1.25) 



G Mp (0) = [Fi(0) + K p F 2 (0)j = (l + «p) =Hp = 2.793/^r , (1-26) 
where \x v = 2.793 is the proton magnetic moment in units of the nuclear magneton 

If we express Fi{Q 2 ) and F 2 {Q 2 ) in terms of Ge p {Q 2 ) and Gm p (Q 2 ), we can write: 



Fi(Q 2 



G Ep (Q 2 ) + tG Mp (Q< 
1 + t 



;i.27) 



Gm p {Q 2 ) — Ge p {Q* 



1.28) 



«p(l + t) 

By substituting for i*i(Q 2 ) and i^(Q 2 ) i n equation (|1.12|) . and dropping the inter- 
ference term between G Ep {Q 2 ) and Gm p {Q 2 ) (keeping leading order terms in pertur- 
bation theory with single-photon exchange), we can express the elastic electron-proton 
cross section in terms of the Sachs form factors as: 



da G\ v m + rG\ 1P m n2 Q 

m=°A (TT7) + 2rG ^ iQ } tan ( -2 ] 



1.29) 



If we express the non-structure cross section in terms of the Mott cross section, 
we can write equation (jl.29|) as: 



da (da\ 1 

In ~ \!n) Mott \e) 1 + t 



(-1- \G 2 Ep (Q 2 ) + r 



l + 2(l+r) tan 2 (-) 



G 2 Mp (Q 2 ) 



. (1.30) 
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The virtual photon longitudinal polarization parameter e is defined as: 



l + 2(l + r) tan 2 



1.31) 



which simplifies equation ()1.30j) to: 



1.32) 



which is the Rosenbluth formula 



Due to the r = ^2 factor in equation (jl.32j) that multiplies G 2 Mp (Q 2 ) but not 



G 2 Ep (Q 2 ), two cases of interest arise: 



1. Small Q 2 (r << 1): the magnetic form factor Gm p {Q 2 ) of the proton is sup- 
pressed and the cross section is dominated by the contribution of the electric 
form factor Ge p {Q 2 ) except as e gets close to zero. 

2. Large Q 2 (r >> 1): the electric form factor G 2 Ep (Q 2 ) of the proton is suppressed 
and the cross section is dominated by the contribution of the magnetic form 
factor G 2 M JQ 2 ). Figure IT31 shows the fractional contribution of Ge p {Q 2 ) to 
the cross section as a function of e for several Q 2 values assuming form factor 
scaling. 

The fact that the magnetic form factor Gm p {Q 2 ) dominates at large Q 2 (small 
distances), makes it difficult to extract Ge p {Q 2 ) with high accuracy from the measured 
cross section at large Q 2 . On the other hand, highly accurate Gm p {Q 2 ) at small Q 2 
is also difficult to extract. 

If we pull out the - factor in equation (jl.32j) . we can write: 



da 



(da\ 



r / e 



dVt \dVtJ Mott V E J 1 



r e \ r 



G 2 Ep (Q 2 ) + G 2 Mp (Q 2 ) 



1.33) 



Finally, we can define the reduced cross section a R as: 
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da (1 + r)e 



rG 2 Mp (Q 2 ) + eG 2 Ep (Q 2 ) 



1.34) 



which is the measured cross section multiplied by a kinematic factor. By measuring 
the reduced cross section <jr at several e points for a fixed Q 2 , a linear fit of ctr to e 
gives tG 2 m p (Q 2 ) as the intercept and G 2 Ep (Q 2 ) as the slope. 

This is known as the Rosenbluth separations method and an example is shown in 
Figure rOl Having extracted Gm p and Ge p for a fixed Q 2 , the ratio of the electric to 
magnetic form factors of the proton can be determined for that Q 2 point. 
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Figure 1.5: Fractional contribution of Ge p (Q 2 ) to the cross section — as a func- 
tion of the virtual photon polarization parameter e for several Q 2 values assuming 
H p Ge p (Q 2 ) = Gm p (Q 2 )- 
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Figure 1.6: Reduced cross section or as a function of the virtual photon polarization 
parameter e at fixed Q 2 = 2.5 GeV 2 . Data is taken from the SLAC NE11 experiment 

ESI- 

1.6 Recoil Polarization Technique 

In recoil polarization experiments [23 123 123 EH]- a longitudinally polarized beam 
of electrons is scattered by unpolarized protons. This results in a transfer of the 
polarization from the electrons to the recoil protons. It must be mentioned that any 
observed polarization of the struck proton was transfered from the electron since the 
target is unpolarized. In elastic electromagnetic scattering in the single-photon ex- 
change approximation, there is no induced polarization. In the l H(e,e',p) reaction, 
see Figure ITTfT the only non vanishing polarization transfer observables are the trans- 
verse, P t , and the longitudinal, Pi, components of the transfered polarization. The 
normal component, P n , does not exist in elastic scattering in single-photon exchange. 

In the single-photon exchange approximation, it can be shown (23123121] that the 
the transverse, Pt, and the longitudinal, Pi, components of the transfered polarization 
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Figure 1.7: Schematic diagram for the recoil polarization spin transfer in the 
l H(e, e',p) reaction. 

are related to the Sachs form factors by: 

IoPi = V^OTTX^Q 2 ) tan 2 (|) , (1.35) 

IoPt = -2y/r(l + T)G Ep (Q 2 )G Mp (Q 2 ) tan(^) , (1.36) 

where E, E', and 9 e are the incident energy, final energy, and scattered angle of the 
electron, and Io is defined as: 

Io = G 2 Ep (Q 2 ) + -G 2 Mp {Q 2 ) . (1.37) 

In the recoil polarization experiments, a Focal Plane Polarimeter (FPP) jSUj is 
used to simultaneously measure both the transverse and longitudinal polarization 
components Pt and Pi. By scattering the proton off a secondary target inside the FPP 
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and measuring the azimuthal angular distribution, both P t and Pi can be determined 
at the same time for a given Q 2 value. 

______ ______ /~i 

By dividing equation (jOfij) by ()1.35|) and solving for we can write: 

^ = §^4^ tan (-), (1-38) 
Gm p Pi 2M P 2 

which is the ratio of the electric to magnetic form factors of the proton as ex- 
tracted by direct and simultaneous measurement of the transverse and longitudinal 
polarization components of the recoiling proton. An interesting point is that while 
the Rosenbluth separations method measures absolute cross sections and then ex- 
tracts Ge p and Gm p from these cross sections, the recoil polarization methods gives 
directly without any measurements of cross sections. The world data for — 
as extracted from the two techniques will be discussed in more detail in the next 
chapter. 



Chapter 2 Previous Form Factor Data 

2.1 Overview 

Understanding the internal structure of the proton is a fundamental problem of 
strong-interaction physics. It was the famous experiment of Hofstadter and collabo- 
rators [SI] at Stanford that first measured the internal structure of the proton. Since 
that time the structure of hadrons has become one of the most important topics in 
nuclear physics and received the attention of experimentalists and theorists world- 
wide. The virtual photon-proton vertex cannot be calculated from first principles. 
Therefore, the internal structure of the proton has been parameterized in terms of 
electric Ge p {Q 2 ) and magnetic Gmp(Q 2 ) form factors. The electromagnetic form fac- 
tors of the protons are used to describe the deviations of the proton from a point-like 
particle in elastic electron-proton scattering. In the non-relativistic limit, the form 
factors are interpreted as the Fourier transform of the spatial distributions of the 
charge and magnetic moment. 

2.2 Previous Measurements 

Several experiments have been conducted over the last fifty years to measure the elas- 
tic electron-proton cross section. Some of these experiments were able to extract the 
electric and magnetic form factors of the proton, Ge p and Gm p , using the Rosenbluth 
separation technique. The form factor ratio M g^ E? has also been measured using the 
recoil polarization technique where the elastic electron-proton cross section are not 
feasible. In this section I will briefly summarize the work done and the results quoted 
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by these two techniques in a chronological order. In particular, a brief description of 
the technique used, kinematics range covered, experimental details, radiative correc- 
tions applied (if applicable), and uncertainties quoted in the measured cross sections 
and in the overall normalization is presented. A comparison between the results of 
the two techniques will be made and a discussion of the global analyses and fits of 
the world data will be presented. Because the focus of the present work is on M ^P Ep , 
the stress is placed on the values of this ratio obtained in the various experiments. 

2.2.1 Elastic e-p Cross Sections Measurements 

The following experiments measured the elastic e-p cross sections for different Q 2 
range. Experiments are listed under the name of the first author: 

• Janssens et al, 1966 \17$ : 

The Stanford Mark III linear accelerator was used to produce an incident electron 
beam of energies in the range of 0.25< E <1.0 GeV. A liquid hydrogen target 0.953 
cm thick with 0.0254 mm stainless steel target walls was used to scatter electrons. 
The scattered electrons were detected using the 72 inch double focusing magnetic 
spectrometer. Elastic e-p cross sections were measured at 25 different Q 2 in the range 
of 0.15< Q 2 <0.86 GeV 2 covering an angular range of 45°< 9 <145° with uncertainty 
never more than 0.08°. Of the 25 Q 2 measured, 20 Q 2 points had enough e coverage 
to do an L-T separation (Rosenbluth separation). Typically 3-5 e points per Q 2 value. 
Measurements for a single e point at Q 2 = 1.01, 1.09, and 1.17 GeV 2 were made at 
constant spectrometer angle of 9 = 145°, and that of Q 2 = 0.49, and 0.68 GeV 2 were 
made at constant spectrometer angle of 9 = 75°. The internal radiative corrections 
were calculated using the method of Tsai j32] and the external radiative corrections 
(bremsstrahlung in the target) were calculated using the formulas of Schwinger [33] 
and Bethe [SI]. The quoted uncertainty in the absolute elastic cross section was on 
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the 4.0% level with 1.6% as an overall normalization uncertainty. L-T extractions 
of the proton form factors were performed. Figure 12.11 shows the ratio of electric to 
magnetic form factor of the proton from this work. 




Q 2 (GeV/c) 2 



Figure 2.1: The ratio of electric to magnetic form factor [i p Ge p /Gmp by Janssens et 
al HZ|. 



Litt et al, 1970 L 35]: 



An electron beam from the Stanford Linear Accelerator Center (SLAC) was pro- 
duced in the range of 4.0< E <10.0 GeV. A liquid-hydrogen target 23 cm in length 
was used to scatter electrons. The scattered electrons were detected using the SLAC 
8-GeV/c magnetic spectrometer. Six Q 2 -values in the range of 1.0< Q 2 <3.75 GeV 2 , 
corresponding to a scattering angle in the range of 12.5°< 9 <41.4°, were covered. On 
the average, 3-5 e points were taken for each Q 2 setting allowing for an L-T extraction 
of the proton form factors. Data was corrected for radiation loss due to straggling of 
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the electrons in the target using Eyges Internal radiative corrections were ap- 
plied using Tsai [32] • Measured cross sections were determined to within a (1.5-2.0)% 
point-to-point uncertainty. An overall normalization uncertainty of 4.0% was quoted. 
Figure 12.21 shows the ratio of electric to magnetic form factor of the proton from this 
work. 
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Figure 2.2: The ratio of electric to magnetic form factor ^ p Ge p /Gmp by Litt et al 
• Price et al, 1971 J73|/ : 

This experiment is an extension of the forward angle experiment done by Goitein 
et al [37J. Incident electron beams from the Cambridge Electron Accelerator (CEA) in 
the range of 0.45< E <1.6 GeV were used and directed on a 3.3 cm liquid-hydrogen 
target in length. Electrons were scattered at large angle in the range of 80°< 9 <90° 
and detected using the 14% total momentum acceptance and 0.83 msr solid angle 
magnetic spectrometer. Cross sections were measured in the range 0.25< Q 2 <1.75 
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GeV 2 . Radiative corrections were applied using the equivalent radiators method of 
Mo and Tsai [3H] with modifications added using Meister and Yennie . The cross 
sections from the large angle measurements were known with uncertainties of (3.1- 
5.3)% including both statistical and systematic uncertainties. There is also an overall 
normalization uncertainty of 1.9%. It should be mentioned that the large angle data 
by itself was not sufficient to do an L-T extraction. Rather, elastic cross sections from 
this work were combined with several e-p scattering experiments and a correction for 
normalization difference between the several experiments was applied. Therefore, no 
estimation of the uncertainties in the cross sections was given due to the difference in 
the normalization procedures used. Global fits were performed using Ge p and Gm p 
as parameters of the fit. The results of combining e-p cross sections from several 
experiments showed deviations from form factor scaling. Figure 12.31 shows the ratio 
of electric to magnetic form factor of the proton from this work. 
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Figure 2.3: The ratio of electric to magnetic form factor [x p Ge p /Gm p by Price et al 
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• Berger et al, 1971 JTfX[ : 

Measurements of the elastic e-p cross sections were made at the Physikalisches 
Institute at the University of Bonn in Germany. Incident electron beams of energies 
in the range of 0.66< E <1.91 GeV were directed on a liquid-hydrogen target 5 
cm in diameter. Cross sections measurements were made for 0.10< Q 2 <1.95 GeV 2 
covering angular range of 25°< 9 <111° with low Q 2 data taken at 9 = 30° in order to 
normalize to other experiments (DESY, Bartel et al 0=0 J. With the goal to combine 
the data and extract form factors, 3-14 e data points were taken for each Q 2 value. 
Internal radiative correction were applied using the method of Meister and Yennie 
|39j . The external bremsstrahlung corrections were applied using Heitler Cross 
sections were determined to within (2-6)% with an overall normalization uncertainty 
of 4%. It was concluded that form factor scaling was not valid at least in the region of 
0.39< Q 2 <1.95 GeV 2 . Figure EU shows the ratio of electric to magnetic form factor 
of the proton from this work. 

• Kirk et al, 1972 U^: 

Electron beams from the Stanford Linear Accelerator Center (SLAC) in the range 
of 4.0< E Q <17.31 GeV were directed on a five condensation-type liquid-hydrogen 
target cells of different sizes (8-32 cm diameter with 25-75- /i-thick stainless steel 
walls). The scattered electrons were detected using the SLAC 1.6-GeV/c and SLAC 
8-GeV/c magnetic spectrometers. With the assumption that the electric form factor 
contribution to the elastic cross section is small, the experiment aimed to extract the 
magnetic form factors over a large Q 2 range of 1.0< Q 2 <25 GeV 2 , covering three 
main angular regions to a maximum of 180°. The 20-GeV/c spectrometer was used 
for small angle measurements 0°< 9 <20°, the 8-GeV/c spectrometer was used for 
intermediate angle measurements 12°< 9 <105°, and the 1.6-GeV/c spectrometer was 
used for backward angle measurements 25°< 9 <165°. The main data was taken in the 



26 



* Berger et al. 




: l l 




iii 



Q 2 (GeV/c) 2 



Figure 2.4: The ratio of electric to magnetic form factor ^ p Ge p /Gm p by Berger et al 

HE- 

range of 12°< 9 <35° while the data at low Q 2 was taken to provide cross calibration 
with other experiments. Some of the higher Q 2 data (Q 2 = 5 and 10 GeV 2 ) was 
taken to be combined with large angle measurements from different experiments to 
provide an upper limit value for Ge p - Internal radiative corrections were applied using 
Tsai j22]; and Eyges [HE] for the external corrections. Cross sections were reported 
with uncertainties of 2.0% and with an overall normalization uncertainty of 4.0%. 
An L-T extraction was impossible since there were not enough e data points taken. 
Although no plot of the ratio of electric to magnetic form factor of the proton from 
this work is possible, the cross sections measured in this experiment were combined 
with cross sections from other experiments for global extractions of the proton elastic 
form factors. 



Murphy et al, 1974 EH'- 
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The University of Saskatchewan Linear Accelerator was used to provide electron 
beams in the range of 0.057< E <0.123 GeV. The beam was directed on a gaseous- 
hydrogen target in a circular cylinder 2.54 cm in radius and 3.50 cm in high. In 
this work, the recoil protons rather than the scattered electrons were detected by a 
double-focusing magnetic spectrometer. Measurements were made for 11-values of 
Q 2 in the range of 0.006< Q 2 <0.031 GeV 2 covering only two angles for the recoiled 
protons, 6 p = 30° and 45°. Not enough e data points were covered to perform an 
L-T extraction. Radiative corrections for the protons were applied using the method 
of Meister and Yennie The experiment extracted the values of Ge p with total 
uncertainty in the range of (0.3-0.9)%. No plot of the ratio of electric to magnetic 
form factor of the proton from this work is possible. Cross sections measured in 
this experiment were combined with cross sections from other experiments for global 
extraction of the proton elastic form factors. 

• Bartel et al, 1973 [L^: 

Elastic e-p scattering cross sections were measured using the Deutsches Elecktronen- 
Synchrotron (DESY) in Hamburg Germany. An electron beams in the range of 
0.8< Eq <3.0 GeV were directed on a cylindrical vessel 5 cm in diameter and 6 
cm long liquid- hydrogen target. Using the high resolution magnetic spectrometer 
(small angle spectrometer), measurements of cross sections were made at electron 
scattering angle in the range of 10°< 8 e <20° and by detecting protons at forward an- 
gles (corresponding to 8 e = 86°) using the recoil nucleon detector. Electrons scattered 
at large angle, 9 e = 86°, were measured using the large angle spectrometer. Cross 
sections were measured at 7-values of Q 2 in the range of 0.67< Q 2 <3.0 GeV 2 with 
uncertainty on the 2-4% level and an overall normalization uncertainty of 2.1%. Typ- 
ically 2-3 £ points were taken per Q 2 value allowing for an L-T extraction. Internal 
radiative corrections were calculated using Meister and Yennie jSH] , and the external 
radiative contribution were calculated using Mo and Tsai jHHl- Electron-proton cross 
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sections in this work and several other experiment were compiled and analyzed. A 
global fit was performed in order to extract the form factors. Deviations from form 
factor scaling were reported. Figure ESI shows the ratio of electric to magnetic form 
factor of the proton from this work. 
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Figure 2.5: The ratio of electric to magnetic form factor h p Ge p /Gm p by Bartel et al 



• Stein et al, 1975 ;43j: 

The Stanford Linear Accelerator Center (SLAC) was used to produce an electron 
beams in the range of 4.5< Eq <20.0 GeV. Elastic cross sections measurements were 
made in the range of 0.004< Q 2 <0.07 GeV 2 at a scattering angle of 9 = 4° providing 
one e data point for each Q 2 value measured. Clearly not enough e data points 
for an L-T extraction. The beam was directed on a vertical cylinder with 0.0076 
cm aluminum walls liquid-hydrogen target. The scattered electrons were detected 
using the SLAC 20-GeV/c spectrometer. Radiative corrections were applied to the 
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measured cross sections using the procedure of Mo and Tsai |2Hj. The uncertainties 
in the cross sections are on the 3.1% level with an overall normalization uncertainty 
of 2.8%. No plot of the ratio of electric to magnetic form factor of the proton from 
this work is possible. Cross sections measured in this experiment were combined with 
cross sections from other experiments for global extractions of the proton elastic form 
factors. 

• Borkowski et al, 1975 \10}'- 

Elastic e-p scattering cross sections have been measured using the 300 MeV Elec- 
tron Linear Accelerator at Mainz. An electron beam was produced in the range of 
0.15< Eq <0.30 GeV and was directed on a 10 mm diameter thin- walled cylindri- 
cal cell filled with liquid-hydrogen. Elastic cross sections measurements were made 
in the range of 0.005< Q 2 <0.183 GeV 2 covering a scattering angle in the range of 
28°< 9 e <150°. Measurements were made using two spectrometers. The first is the 
double-focusing 180° spectrometer and was set to different scattering angles. The 
second consisted of two quadrupoles and 12° b ending-magnet spectrometer and was 
fixed at scattering angle of 28° as a monitor. Radiative corrections were applied using 
Mo and Tsai [SB] • Random uncertainties in the cross sections were reported and were 
on the (1.0-2.0)% level. Figure l2~o1 shows the ratio of electric to magnetic form factor 
of the proton from this work. 

• Bosted et al, 1990 

The primary goal of this experiment was to measure the magnetic structure func- 
tion B(Q 2 ) of the deuteron to the largest Q 2 -value possible. The elastic cross sections 
of the e-p scattering were measured near 180° and the value of the magnetic form 
factor Gmp was determined assuming form factor scaling with small Ge p contribution 
to the cross sections at high Q 2 points. Incident electron beam from the Stanford 



30 



* 1 - 3 I 1 

o 

jj w 1 .2 — * Borkowski et al. 

a 
=L 

1.1 

T T 

— i ill--- 

0.9 
0.8 

0.7 I i i i i i i i i i i_l 

-2 -1 

10 10 1 

Q 2 (GeV/c) 2 

Figure 2.6: The ratio of electric to magnetic form factor H p Ge p /Gm p by Borkowski 
et al P2|. 

Linear Accelerator Center (SLAC) of energies in the range of 0.5< Eq <1.3 GeV was 
directed on a liquid-hydrogen target with nominal lengths of 40, 20, 10, and 5 cm. 
Cross sections were measured at 11- values of Q 2 in the range of 0.49< Q 2 <1.75 GeV 2 
for electrons backscattered near 180° in coincidence with protons recoiling near 0° in 
a large solid-angle double-arm spectrometer. An L-T extraction is impossible from 
these measurements. The equivalent radiator approximation method of Tsai [32] was 
used for the internal radiative corrections. The effect of the bremsstrahlung and Lan- 
dau straggling in the target external radiative corrections were combined with that 
of the internal corrections. Cross sections with uncertainty of 3% were quoted with 
an overall normalization uncertainty of 1.8%. No plot of the ratio of electric to mag- 
netic form factor of the proton from this work is possible. Cross sections measured in 
this experiment were combined with cross sections from other experiments for global 
extraction of the proton elastic form factors. 
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• Rock et al, 1992 J^j : 

The main goal of this experiment was to extract the elastic neutron cross sections 
but elastic e-p cross sections were also measured since they were needed for the 
analysis. The Stanford Linear Accelerator Center (SLAC) electron beam with energies 
in the range of 9.761< E <21.0 GeV was directed on a 30-cm long liquid-hydrogen 
cell. The elastic e-p cross sections were measured at 5 values of Q 2 in the range of 
2.5< Q 2 <10.0 GeV 2 at a fixed scattering angle of 9 = 10°. An L-T extraction is 
impossible with this data. The scattered electrons were detected using the SLAC 20- 
GeV/c spectrometer. The cross sections were radiatively corrected using the method 
of Tsai [S2] ■ The uncertainties in the cross sections were on the 1-4% level in addition 
to an overall normalization uncertainty of 3%. No plot of the ratio of electric to 
magnetic form factor of the proton from this work is possible. Cross sections measured 
in this experiment were combined with cross sections from other experiments for 
global extraction of the proton elastic form factors. 

• Sill et al, 1993 

Elastic e-p cross sections were measured using the beam line at the Stanford Linear 
Accelerator (SLAC). An electron beams of energies in the range of 5.0< Eq <21.5 
GeV were directed on two liquid-hydrogen targets with different lengths. The 25-cm 
target was used for normalization of the acceptance of the 65-cm target and to provide 
a test for the low Q 2 data. The 65-cm target provided higher counting rate and was 
used to take the majority of the elastic data. The scattered electrons were detected 
using the SLAC 8-GeV/c spectrometer. Cross sections were measured for 13-values 
of Q 2 for the range of 2.9< Q 2 <31.3 GeV 2 covering three-angle settings of 9 e = 21°, 
25°, and 33°. Due to the limited angular range, an L-T extraction was impossible 
to perform. Radiative corrections were applied to the data using Mo and Tsai |38j . 
Cross sections uncertainties were on the 3-4% level with an overall normalization 
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uncertainty of 3.6%. No plot of the ratio of electric to magnetic form factor of the 
proton from this work is possible. Cross sections measured in this experiment were 
combined with cross sections from other experiments for global extraction of the 
proton elastic form factors. 

• Walker et al, 1994 

The Stanford Linear Accelerator Center beam line was used to produce incident 
electron beams of energies in the range of 1.594< E Q <8.233 GeV. The beam was di- 
rected on cylindrical liquid-hydrogen target 20-cm in length and 5.08-cm in diameter. 
The scattered electrons were detected using the SLAC 8-GeV/c spectrometer. Elastic 
e-p cross sections were measured for 4-values of Q 2 in the range of 1.0< Q 2 <3.007 
GeV 2 covering angular range of 11.714°< 9 <45.221° with an average of 3-8 e points 
per Q 2 point allowing for an L-T extraction. Internal radiative corrections were done 
using Mo and Tsai [38j. Also, improvements were made to the internal corrections 
using the equivalent radiator approximation. The external corrections were applied 
using the work of Tsai to account for bremsstrahlung in the target material and 
the effects of the Landau tail of the ionization energy loss spectrum. Uncertainty in 
the cross sections was on the 1% level with an overall normalization of 1.9%. Figure 
12.71 shows the ratio of electric to magnetic form factor of the proton from this work. 
Cross sections from several experiments including this work were combined and a fit 
for a global extraction of the form factors was performed. Results indicate good con- 
sistency between the different data sets. The form factors extracted from this work 
supported form factor scaling. See Figure l2~T8l in section ()2.4|) for more detail. 

• Andivahis et al, 1994 

The primary goal of this experiment was to minimize both statistical and sys- 
tematic uncertainties and extend the measurements of the form factors to a high Q 2 
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Figure 2.7: The ratio of electric to magnetic form factor ^ p Ge p IGm p by Walker et al 

values. The Stanford Linear Accelerator Center (SLAC) was used to produce electron 
beams with energies in the range of 1.511< E <9.8 GeV. The beam was directed 
on a 15-cm liquid-hydrogen target. The scattered electrons were detected using both 
the SLAC 1.6 and 8-GeV/c spectrometers simultaneously. The 8-GeV/c spectrometer 
was used to measure 5-values of Q 2 in the range of 1.75< Q 2 <5.0 GeV 2 where an L-T 
extraction is possible with an average of 3-6 e points per Q 2 . In addition, two extra 
points at Q 2 = 6.0 and 7.0 GeV 2 were measured as a single e points. The angular 
range was 13.25°< 9 <90.066°. The 1.6-GeV/c spectrometer was used to measure 
8-single-points of e at Q 2 = 1.75, 2.50, 3.25, 4.0, 5.0 6.0, 7.0, and 8.83 GeV 2 at a 
constant angle 9 ~ 90° so data could be combined with previous forward-angle cross 
sections |1HJ d] collected at SLAC at the same Q 2 . The 1.6-GeV/c cross sections 
were normalized to the 8-GeV/c results. The procedure of the radiative corrections 
applied was similar to the procedure introduced by Walker et al [T3] reported above. 
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Average uncertainties in the cross sections were less than 2% with an overall normal- 
ization uncertainty of 1.77%. Figure I2~51 shows the ratio of electric to magnetic form 
factor of the proton from this work. 
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Figure 2.8: The ratio of electric to magnetic form factor ii p Ge p /Gm p by Andivahis 
et al |TS]. 

• Christy et al, 2004 

The focus of this experiment was to separate the longitudinal and transverse 
unpolarized proton structure functions in the resonance region using Rosenbluth sep- 
aration technique, but elastic e-p cross sections were also measured at 28-values of 
Q 2 in the range of 0.4< Q 2 <5.5 GeV 2 covering an angular range of 12.5°< 9 e <80°. 
This range covered 3 e points per Q 2 allowing for L-T separations at 7 Q 2 -values. The 
experiment was done at the experimental Hall C of the Thomas Jefferson National 
Laboratory (JLAB) in Newport News Virginia. An electron beam in the range of 
1.148< E <5.494 GeV was directed on a tuna-can shaped cryogen liquid hydrogen 
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cell. The cell has an inside diameter of 40.113 mm when warm and 39.932 mm when 
cold with cylindrical wall thickness of 0.125 mm. The scattered electrons were de- 
tected using the High Momentum Spectrometer. Radiative corrections were applied 
using the procedure described in Walker et al ^3] and based on the prescription of 
Mo and Tsai [SHI- Uncertainties in the cross sections were on the 1.96% level with an 
overall normalization uncertainty of 1.7%. The uncertainties in the ratio in the 



Mp 



Q 2 region of interest are quite large and do not represent an improvement in accuracy 
over previous L-T determinations. However, when the polarization transfer results 
became available and were at odds with the previously known Rosenbluth results, 
Christy et al carefully analyzed their elastic e-p data and concluded that their 
were not consistent with the ratio reported by Jones et al. Figure l2~9"l shows the ratio 
of electric to magnetic form factor of the proton from this work. 




Q 2 (GeV/c) 2 



Figure 2.9: The ratio of electric to magnetic form factor [i p Ge p I 'Gmp Christy et al 
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2.2.2 Polarization Measurements 

The following experiments represent the proton polarization measurements (recoil 
polarization and polarized target measurements) to date for different Q 2 range. 

• Alguard et al, 1976 

The main idea of this experiment was to measure the antiparallel-parallel asym- 
metry A in the differential cross section which in turn is related to the form factors of 
the proton. The polarized electron source (PEGGY) at the 20-GeV Stanford Linear 
Accelerator Center (SLAC) was used to produce a polarized electron beam which was 
directed on a polarized proton target polarized by the method of dynamic nuclear ori- 
entation in a butanol target doped with 1.4% porphyrexide. The scattered electrons 
were detected using the 8-GeV/c spectrometer. Data were taken at Q 2 = 0.765 GeV 2 
with incident electron energy of Eq = 6.473 GeV and scattering electron angle of 9 e 
= 8.005°. The beam and target polarization P e and P p were measured. The exper- 
imental asymmetry A was determined in order to solve for the antiparallel-parallel 
cross section asymmetry A where A = P e P p FA and F is the fraction of the elastically 
scattered electrons within the elastic missing- mass region. Figure Yl. 1 01 shows the ratio 
of electric to magnetic form factor of the proton from this work. 

• Milbrath et al, 1998 L 23]: 

This is the first experiment to demonstrate the feasibility of the recoil polariza- 
tion as a technique to extract the ratio of electric to magnetic form factors of the 
proton. Measurements of the recoil proton polarization observables in the reactions 
pie, e'p) and d(e, e'pjn were made. The MIT-Bates Linear Accelerator Center was 
used to produce a longitudinally polarized electron beam of energy E = 0.58 GeV. 
The beam was directed on an unpolarized cryogenic target of liquid hydrogen and 
deuterium cells of 5 and 3 cm in diameter respectively. The scattered electrons were 
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Figure 2.10: The ratio of electric to magnetic form factor \x v G ' e p i 'G 'm p by Alguard et 
al HE]. 

detected using the Medium Energy Pion Spectrometer, and the scattered protons 
were detected using One-Hundred Inch Proton Spectrometer. Two Q 2 -values of 0.38 
and 0.5 GeV 2 were measured corresponding to electron-scattering angles of 9 e = 82.7° 
and 113° respectively. The recoil proton polarization was measured in the focal plane 
polarimeter (FPP). Figure 12. Ill shows the ratio of electric to magnetic form factor of 
the proton from this work. 

• Jones et al, 2000 \2~$: 

Recoil proton polarization measurements were carried out at the experimental 
Hall A of Thomas Jefferson National Laboratory (JLAB) in Newport News Virginia. 
Polarized electron beam of energy in range of 0.934< Eq <4.090 GeV was directed on 
a 15-cm-long unpolarized liquid hydrogen target. The elastically scattered electrons 
and protons were detected in coincidence using the two identical high resolution 
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Figure 2.11: The ratio of electric to magnetic form factor ^ p Ge p /Gmp by Milbrath 
et al 



spectrometers (HRS) of Hall A. The ratio ^|r~f was determined at 9 Q 2 - values in the 
range of 0.49< Q 2 <3.47 GeV 2 covering an angular range of 22.26°< 9 e <79.88° for 
the electrons and 29.98°< P <46.03° for the protons. The recoil proton polarization 
was measured in the focal plane polarimeter (FPP). External radiative corrections 
were not applied. The internal radiative corrections such as the hard photon emission 
and higher order contributions were calculated and found to be on the order of 
a few percent and were not applied. The results of this work showed the decline of 
M Gmp P w ^ n increasing Q 2 deviating from form factor scaling and indicating for the 
first time a definite difference in the spatial distribution of charge and magnetization 
currents in the proton. Figure 12". 1 21 shows the ratio of electric to magnetic form factor 
of the proton from this work. 



Pospischil et al, 2001 L 50]: 
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Figure 2.12: The ratio of electric to magnetic form factor \x v Ge p IGm v by Jones et al 

Recoil proton polarization measurements were carried at the 3-spectrometer setup 
of the Al-Collaboration at the Mainz microtron MAMI. Longitudinally polarized 
electron beam of energy Eq = 0.8544 GeV was directed on a 49.5-mm-long Havar 
cell filled with unpolarized liquid hydrogen. Data were taken at 3 Q 2 -values in the 
range of 0.373< Q 2 <0.441 GeV 2 covering an angular range of 48.2°< 9 e <54.4° for 
the electrons and 45.5°< 8 P <49.5° for the protons. The recoil proton polarization 
was measured in the focal plane polarimeter (FPP). Radiative corrections were not 
applied. Figure 12.131 shows the ratio of electric to magnetic form factor of the proton 
from this work. 

• Gayou et al, 2001 Jl?gj/ : 

Recoil proton polarization measurements were carried out at the experimental Hall 
A of Thomas Jefferson National Laboratory (JLAB) in Newport News Virginia using 
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Figure 2.13: The ratio of electric to magnetic form factor IL p Ge p /Gm p by Pospischil 
et al |3g. 

the same experimental hardware used by Jones et al [21]. Polarized electron beam 
of energy in range of 1.0< E <4.11 GeV was directed on a 15-cm-long unpolarized 
liquid hydrogen target. The elastically scattered electrons and protons were detected 
in coincidence using the two identical high resolution spectrometers (HRS) of Hall 
A. Although the goal of this experiment was to study the D(j,p)n and H(ff, p)7r° 
reactions, 13 measurements of coincidence ep — > ep polarizations were performed to 
calibrate the focal plane polarimeter (FPP) used to study the reactions defined above. 
Data were taken in the range of 0.32< Q 2 <1.76 GeV 2 covering an angular range of 
18.012°< 6 e <48.65° for the electrons and 37.446V 6 P <55.984° for the protons. The 
recoil proton polarization was measured in the focal plane polarimeter. Radiative 
corrections were not applied. Figure |2~. 141 shows the ratio of electric to magnetic form 
factor of the proton from this work. 
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Figure 2.14: The ratio of electric to magnetic form factor (J> p Ge p /Gmp by Gayou et 
al (first measurements). 

• Gayou et al, 2002 L 26]: 

This work in an extension of that of Jones et al 24 j where measurements were 
taken to higher Q 2 value of 5.54 GeV 2 . Recoil proton polarization measurements 
were carried out at the experimental Hall A of Thomas Jefferson National Laboratory 
(JLAB) in Newport News Virginia also using the same experimental hardware used 
by Jones et al. Polarized electron beam of energy in range of 4.59< E <4.607 
GeV was directed on a 15-cm-long unpolarized liquid hydrogen target. Data were 
taken in the range of 3.5< Q 2 <5.54 GeV 2 . At Q 2 = 3.5 GeV 2 , the elastically 
scattered electrons and protons were detected in coincidence using the two identical 
high resolution spectrometers (HRS) of Hall A with fixed angles of 8 e = 30.6° for the 
electrons and 8 P = 31.785° for the protons. At higher Q 2 values of 3.97, 4.75, and 
5.54 GeV 2 and at a fixed beam energy E = 4.607 GeV, the electrons scattered at a 
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higher angles than the protons and were detected using a calorimeter in coincidence 
with the protons covering an angular range of 19.275°< 9 e <28.587° for the electrons 
and 34.5°< 6 P <51.44° for the protons. The recoil proton polarization was measured 
using a focal plane polarimeter. Radiative corrections were not applied. Figure 12.151 
shows the ratio of electric to magnetic form factor of the proton from this work. 
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Figure 2.15: The ratio of electric to magnetic form factor [i p Ge p /Gmp by Gayou et 
al (second measurements). 



2.3 Summary of Previous e-p Measurements 

A summary of section 12.21 for the world's data on elastic e-p scattering cross sec- 
tion measurements is given in Table 12.11 Experiments are listed under the principal 
author's name, laboratory at which they were performed, energy, Q 2 , number of e 
points measured at each Q 2 value N E , cross section uncertainty Act, and the overall 
normalization uncertainty in the cross section Acr/v- 
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Author 


Laboratory 


Energy 


Q 


N e 




Aa N 






(GeV) 


(CpV 2 ) 


















% 


% 


JanssensjTTj 


Mark III 


0.250-1.00 


0.15-0.86 


3-5 


4 


1.6 


Litt0Sl 


SLAC 


4.00-10.00 


1.00-3.75 


3-5 


1.5-2.0 


4 


Price [IE! 


CEA 


0.45-1.60 


0.25-1.75 


1 


3.1-5.3 


1.9 


Berger [To] 


Bonn 


0.66-1.91 


0.10-1.95 


3-14 


2-6 


4 


Kirk [12] 


SLAC 


4.00-17.31 


1.00-25.00 


1 


2 


4 


Murphy [TT] 


Saskatchewan 


0.057-0.123 


0.006-0.031 


1 






Bartel 14 


DESY 


0.80-3.00 


0.67-3.00 


2-3 


2-4 


2.1 


Stein US] 


SLAC 


4.50-20.00 


0.004-0.07 


1 


3.1 


2.8 


BorkowskijlOj 


Mainz 


0.15-0.30 


0.005-0.183 




1-2 




BostedHU 


SLAC 


0.50-1.30 


0.49-1.75 


2-7 


3 


1.8 


Rockg3] 


SLAC 


9.761-21.00 


2.50-10.00 


1 


1-4 


3 


Sill [IH] 


SLAC 


5.00-21.50 


2.90-31.30 


1 


3-4 


3.6 


Walker [IH1 


SLAC 


1.594-8.233 


1.00-3.007 


3-8 


1 


1.9 


Andivahis|15j 


SLAC 


1.511-9.80 


1.75-5.00 


3-6 


< 2 


1.77 


Christy |Hj 


JLAB 


1.148-5.50 


0.40-5.50 


3 


1.96 


1.7 



Table 2.1: Summary of selected world data on e-p elastic scattering cross section 
measurements. 
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2.4 Discussion 

The world data on the ratio of electric to magnetic form factor of the proton as ex- 
tracted from elastic e-p cross sections (Rosenbluth separation) and recoil polarization 
measurements in the single-photon exchange approximation have been summarized in 
section l2~2l Figure l2~. 161 shows the M ^P gp world data as determined by the Rosenbluth 
separation method (notice the logarithmic scale for the Q 2 axis). With the excep- 
tion of the data of Berger ^H] and Bartel [Tl], which show a decrease in the ratio 
with increasing Q 2 especially for the Q 2 > 1.0 GeV 2 , the data show approximate 
form factor scaling or = l.o with some fluctuation at high Q 2 points where 

the uncertainties become large. This is due to the fact that the experimental cross 
section is small and dominated by Gm p which makes it difficult to extract Ge p with 
high precision. The difficulty in extracting Ge p at high Q 2 values with high precision 
using the Rosenbluth separation technique was the motivation behind developing the 
recoil polarization technique, which is more sensitive to Ge p at large Q 2 . 

Figure 12.171 shows the world data on ^± as determined by recoil polarization 
technique. The data agree with form factor scaling for the region Q 2 < 1.0 GeV 2 . 
However, for the region Q 2 > 1.0 GeV 2 , the data decrease with increasing Q 2 (notice 
the logarithmic scale for the Q 2 axis) deviating significantly from form factor scaling as 
suggested by Rosenbluth separation. The data from recoil polarization measurements 
are more precise at high Q 2 and maybe less sensitive to systematic uncertainties than 
the Rosenbluth data. The dashed line in Figure 12.171 is the recoil polarization fit to 
the data or: 

^^ = l-0.13(Q 2 -0.04) . (2.1) 

Gm p 

We can see from Figures f2. 161 and 12. 171 that the two techniques give different re- 
sults in the region Q 2 > 1.0 GeV 2 . The values of ^t^ 2 from the two techniques 
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Figure 2.16: The world data on the ratio of electric to magnetic form factor 
^ p Ge p /Gm p by Rosenbluth separation. 
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Figure 2.17: The world data on the ratio of electric to magnetic form factors 
H p Ge p /Gm p by recoil polarization. The dashed line is the recoil polarization pub- 
lished fit. 
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differ almost by a factor of three at the high Q 2 points. This discrepancy between 
the results of the two techniques is becoming to be known as the G P E crisis. This dif- 
ference implies uncertainties in our knowledge of the form factors of the proton and 
raises several questions that must be answered. For example, the high precision data 
provided by the recoil polarization technique and the larger uncertainties in the L-T 
data (small Ge p contribution at high Q 2 ) have led people to believe that the previous 
Rosenbluth extractions are inconsistent. If this is the case, then all the form factors 
extracted using Rosenbluth separation technique which are supposed to parameterize 
the deviation of the proton's structure from point-like particle are unreliable. Also, 
if there is a significant error in the elastic e-p cross section measurements (Q 2 > 1.0 
GeV 2 ), then there could be errors in all previous experiments that require normaliza- 
tion to elastic e-p cross sections, or require the use of the elastic cross sections [SB El] 
or form factors [53 as an input to the analysis. If the cross section measurements and 
hence the Rosenbluth extractions are incorrect, that still will not solve the problem 
since the recoil polarization technique provides the ratio of Ge p to Gm p and not the 
actual values for the individual form factors. 

At this point, the following questions are paramount: 

1. Why do the two techniques disagree? 

• Is there a missing correction (radiative corrections or normal- 
ization uncertainties ) or something wrong in the analysis of the 
previous Rosenbluth separations data? 

• Is there a missing correction (radiative corrections or proton 
spin precession determination) or something wrong in the anal- 
ysis of the previous recoil polarization data? 

2. Which form factors are the correct ones to use? 
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• Is there something fundamentally wrong in the physics of one 
or both techniques? 

3. What about all the conclusions based on the old and new theoretical 
models and calculations concerning the proton form factors ? 

In order to provide a check on the consistency of the world's cross section mea- 
surements, several global analysis have been done. Walker [IB*] lo^j combined cross 
sections from different experiments and performed a global extraction of the elastic 
form factors in the range of 0.1 < Q 2 < 10.0 GeV 2 . The result of the global fit 
supported the ansatz of form factor scaling as shown in Figure 12.181 

An empirical fit to the world data of the proton form factors was made by Bosted 
[HI for the region 0.0 < Q 2 < 30.0 GeV 2 . The data of both §^ and j n the region 
Q 2 < 7.0 GeV 2 are from the global analysis of Walker [T3], while for Q 2 > 9.0 GeV 2 , 
form factor scaling was assumed to extract Gmp- A good fit to the data was achieved 
when the form factors were described as an inverse polynomial in Q: 

Gep ^ = 1.0 + 0.62Q + 0.68Q 2 + 2.80Q 3 + 0.83Q 4 ' ^ 



G Mp (Q 2 ) _ i.o 



(2.3) 



u p 1.0 + 0.35Q + 2.44Q 2 + 0.50Q 3 + 1.04Q 4 + 0.34Q 5 ' 

An extensive examination of the form factors extractions from cross section mea- 
surements was done by Arrington EH] ■ A similar global analysis to that of Walker 
was performed. The fit included two extra data sets from Stein jl^] and Rock [23] 
in addition to the data sets used in Walker's analysis. Also, the final published cross 
sections in the work of Sill [19j and Andivahis [15 j were used which were not available 
for Walker's analysis. Some recent measurements of elastic scattering at Jefferson 
Lab by Dutta [HI], Niculescu [201, and Christy [17] were added. Also, the results 
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of Borkowski ^H], Murphy jllj . and Simon were added to constrain the low Q 2 
behavior. In addition, all the high Q 2 data up to 30 GeV 2 were included in the fit. 

For each dataset, an overall normalization uncertainty was determined. The nor- 
malization uncertainties were taken either from the original published work or from 
Walker's global analysis. Independent normalization uncertainties were assigned to 
data taken by different detectors in the same experiment. In the work of Bartel |4Uj . 
data were taken using three different spectrometers, so the data were divided into 
three sets with different normalization uncertainty factor assigned to each data set. 
Higher order terms such as the Schwinger term and additional corrections for 
vacuum polarization contributions from muon and quark loops were added to the 
radiative corrections applied to the work of Janssens J7j, Bartel |4TH 15%] HI]. Alber- 
cht jSH], Litt [US], Goitein [37], Berger [TH], and Price [TH]. Finally, the small-angle 
(< 20°) data from Walker ^3] were excluded since an error was identified in that 
data. The form factors are parameterized as: 



Ge p {Q 



2 . G Mp (Q 2 ) 



1.0 



Hp 1.0 + p 2 Q 2 + p 4 Q 4 + . . . + P2N Q 2N 
where the parameters of the fit are listed in Table 12.21 



(2.4) 



Parameter 


Gep 


(-*Mp 


P2 


3.226 


3.19 


P4 


1.508 


1.355 


P6 


-0.3773 


0.151 


P8 


0.611 


-1.14 x 10- 2 


PlO 


-0.1853 


5.33 x 10" 4 


Pl2 


1.596 x 10" 2 


-9.0 x 10" 6 



Table 2.2: Form factors parameters, equation (|2.4p . as extracted by Arrington fit to 



a. 



It was hypothesized that the discrepancy between the Rosenbluth and the po- 
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larization data was coming from a common systematic error in the cross section 
measurements and a (5-8)% e-dependent systematic error in the cross sections could 
resolve the discrepancy (see section |H3 for more detail). Therefore, a combined anal- 
ysis was done by Arrington [22] where an ^-dependent correction of 6% was applied 
to all cross sections: 

a c = aJl.O - 0.06(e - 1.0)) , (2.5) 

and then the recoil polarization data were included in the fit. Here a c and a are 
the corrected and uncorrected cross sections respectively. The form factors were 
parameterized using the same form as in equation (|2.4j) . The parameters of the fit 
are listed in Table 1231 



Parameter 


Gep 




P2 


2.94 


3.00 


P4 


3.04 


1.39 


P6 


-2.255 


0.122 


P8 


2.002 


-8.34 x 10~ 3 


PlO 


-0.5338 


4.25 x 10~ 4 


Pl2 


4.875 x 10" 2 


-7.79 x 10~ 6 



Table 2.3: Form factors parameters as extracted from Arrington combined fit to cross 
sections and recoil polarization data. 



The ratio of the electric to magnetic form factor of the proton as extracted from 
the Arrington's global fit to the world's elastic cross section data supports the results 
of previous Rosenbluth extractions. The result indicated a good consistency between 
all different data sets and ruled out any possibility for a single or two bad data sets 
or incorrect normalization in the combined Rosenbluth analysis. 

An empirical fit to the proton form factors was done by Brash jU] where most 
of the higher-Q 2 elastic e-p cross sections data were reanalyzed by using the ratio 
between the proton form factors r = ^ Gfip as provided by the recoil polarization 



'Mp 
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data constraint: 



r(Q 2 ) 



G 



Mp 



1 -0.130(Q 2 -0.04) 







: if Q 2 < 0.04 GeV 2 , 

: if 0.04 < Q 2 < 7.7 GeV 2 , 

: if Q 2 > 7.7 GeV 2 , 



(2.6) 



that way, r(Q 2 ) fixes the ratio between the intercept (a = rG 2 Mp ) and the slope 
(b = G 2 Ep ) from the linear fit of the reduced cross section <jr to e. The form factors are 
associated with the parameters of the linear fit, <7r = a(l + be), where a 



bR 



and R = -±f. A new parameterization of the proton magnetic form factor Gm p was 
obtained and the electric form factor was then calculated using the recoil polarization 
constrained ratio: 



Gm p (Q 2 



1 



1.0 + 0.116Q + 2.874Q 2 + 0.241Q 3 + 1.006Q 4 + 0.345Q 5 



(2.7) 



G Ep (Q 2 )=r(Q 2 )^^ 1 - (2-8) 
Hp 

The magnetic form factor of the proton as parameterized by Brash jU] was shown 
previously in Figure 11.41 in section 11.41 In Figures 12.181 12.191 and 12.201 the fits of 
Arrington, Bosted, and recoil polarization are shown along with the world data on 
proton form factors in addition to Walker's global analysis. 
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Figure 2.18: Global analysis of the proton form factors ratio [i p Ge p /Gm p by Walker. 
In addition, the fits of Arrington, Bosted, and recoil polarization are also shown. 
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Figure 2.19: Proton form factors ratio h p Ge p /Gm p by Rosenbluth separations. In 
addition, the fits of Arrington and Bosted are also shown. 
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Figure 2.20: Proton form factors ratio h p Ge p /Gm p by Rosenbluth separations and 
recoil polarization. In addition, the fits of Arrington, Bosted, and recoil polarization 
are also shown. 



Chapter 3 Experimental Setup 



3.1 Overview 

Due to the inconsistency in the results of the ratio of the electric to magnetic form 
factors of the proton, ^rf 2 , as extracted from the Rosenbluth and recoil polariza- 
tion techniques, and due to the fact that the reported uncertainties in in the 
various Rosenbluth determinations are much larger than those quoted for the polar- 
ization transfer measurements, a high-precision measurement of fl ^ Sp using the L-T 
separation technique in the Q 2 > 1.0 GeV 2 region is important to: 

• Provide a comparison between the two techniques in the region where 
both can extract the value of M ^P fip with high precision. 

• Achieve uncertainties comparable to or better than the uncertainties 
quoted by the recoil polarization measurements. 

• Provide a check on any possibility of additional and unaccounted 
for systematic uncertainties in the L-T or recoil polarization mea- 
surements. 

Experiment E01-001 sought to achieve these goals. It ran in May 2002 and was 
carried out in Hall A of the Thomas Jefferson National Accelerator Facility (formerly 
known as Continuous Electron Beam Accelerator Facility, or CEBAF) which is lo- 
cated in Newport News Virginia in the U.S.A. An incident electron beam of energies 
in the range of 1.912< E <4.702 GeV was directed on a 4-cm-long unpolarized liq- 
uid hydrogen target. High precision measurements of the elastic e-p cross sections 

53 
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were made to allow for an L-T separation of the proton electric and magnetic form 
factors. Protons were detected simultaneously using the two identical high resolution 
spectrometers (HRS) or what is known by the left and right arm spectrometers of 
Hall A. The left arm spectrometer was used to measure three Q 2 points of 2.64, 3.20, 
and 4.10 GeV 2 . Simultaneously, measurements at Q 2 = 0.5 GeV 2 were carried out 
using the right arm spectrometer which served as a monitor of beam charge, current, 
and target density fluctuations. 

A total of 12 points (5 e points for Q 2 = 2.64 GeV 2 , 4 e points for Q 2 = 3.20 GeV 2 , 
and 3 e points for Q 2 = 4.10 GeV 2 ) were measured covering an angular range of 12.52° 
< 6l < 38.26° for the left arm, while the right arm was at Q 2 = 0.5 GeV 2 , and used to 
simultaneously measure 5 e points covering an angular range of 58.29° < Or < 64.98°. 
Here, 0^ and Or are the nominal angle of the struck proton with respect to the beam 
electron for the left and right spectrometer, respectively. Figure 13.11 and Table 13.11 
show and list the nominal kinematics covered in the E0 1-001 experiment and their 
settings. Small offsets were determined and applied to the energy and scattering 
angles. See section I4~4l for details. The final kinematics used in the analysis are listed 
in Table IP 

The final momentum of the scattered protons, P me asured, measured using the high 
resolution spectrometer, was compared to the final momentum of the scattered pro- 
tons, PcaicuiatediQp) j calculated from two-body kinematics using the measured scatter- 
ing angle of the protons V (see equation ()5.8j0 : 

2E beam (M 2 + E beam M p ) cos Op 

Pcalculate^) ~ M 2 + £ 2 _ + ^ beam M p - E 2 eam COS 2 p ' ^ 

and the difference in momentum AP was then constructed: 

AP Pmeasured PcalculatediOp) , (3.2) 




Figure 3.1: Plot of e vs Q 2 showing the kinematics covered. 

The AP spectrum is made of several contributions. The main contributions come 
from the elastic peak which is due to the elastic ep — > ep scattering and the radia- 
tive tail. In addition there are backgrounds due to quasi-elastic scattering from the 
aluminum target windows and high energy protons generated from photoreactions 
{lP ~* 7r °P an d 7£> — > 7p) that contribute to the AP spectrum. For each of the 
kinematics covered, data were taken with the dummy target to subtract away the 
endcaps contribution from the spectrum. The photoproduction of tt° or 7^ — > 7i°p 
events were simulated using a calculated bremsstrahlung spectrum and ^ oc s~ 7 and 
then subtracted away from the spectrum as well. 

The net number of elastic events from data is then compared to the number of 
elastic events in the e-p peak as simulated using the Monte Carlo simulation program 
SIMC |dTH 1ST] , under the same conditions (cuts), for a given narrow window cut 
on the AP spectrum. The ratio of the number of events from the data to that of 
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Setting 


Ebeam 


£l 


Ql 


Ql 


Pl 


£r 


0R 


Pr 




(GeV) 




(GeV 2 ) 


(°) 


(GeV/c) 




(°) 


(GeV/c) 





1.912 


0.117 


2.64 


12.631 


+ 2.149 


0.914 


58.288 


+ 0.756 


a 


2.262 


0.356 


2.64 


22.166 


+ 2.149 


0.939 


60.075 


+ 0.756 


i 


2.842 


0.597 


2.64 


29.462 


+ 2.149 


0.962 


62.029 


+ 0.756 


1 


3.772 


0.782 


2.64 


35.174 


+ 2.149 


0.979 


63.876 


+ 0.756 


I 


4.702 


0.865 


2.64 


38.261 


+ 2.149 


0.986 


64.978 


+ 0.756 


b 


2.262 


0.131 


3.20 


12.525 


+ 2.471 


0.939 


60.075 


+ 0.756 


J 


2.842 


0.443 


3.20 


23.395 


+ 2.471 


0.962 


62.029 


+ 0.756 


P 


3.772 


0.696 


3.20 


30.501 


+ 2.471 


0.979 


63.876 


+ 0.756 


m 


4.702 


0.813 


3.20 


34.139 


+ 2.471 


0.986 


64.978 


+ 0.756 


k 


2.842 


0.160 


4.10 


12.682 


+ 2.979 


0.962 


62.029 


+ 0.756 


r 


3.772 


0.528 


4.10 


23.665 


+ 2.979 


0.979 


63.876 


+ 0.756 


n 


4.702 


0.709 


4.10 


28.380 


+ 2.979 


0.986 


64.978 


+ 0.756 


coinl 


2.262 


0.356 


2.64 


22.166 


+ 2.149 


0.356 


71.481 


- 0.855 


coiril 


2.842 


0.597 


2.64 


29.462 


+ 2.149 


0.597 


47.439 


- 1.435 


coin?) 


3.362 


0.398 


4.10 


20.257 


+ 2.185 


0.398 


61.184 


- 1.177 



Table 3.1: The nominal kinematics setting for the E01-001 experiment. For each 
setting, E beam is the nominal electron beam energy, e L (er) is the virtual photon po- 
larization parameter for the left (right) arm spectrometer, Q\ is the four-momentum 
transfer squared for the left arm spectrometer, 6^ (Or) is the spectrometer nominal 
scattering angle for the left (right) arm spectrometer, and Pl {Pr) is the central mo- 
mentum for the left (right) arm spectrometer. Note that the four-momentum transfer 
squared for the right arm Q\ (not listed) was at 0.50 GeV 2 for all the settings ex- 
cept for the coincidence kinematics settings coinl, coiril, and coin3 where Q 2 R = Q\. 
Small offsets were determined and applied to the energy and scattering angles. See 
section 14.41 for details. The final kinematics used in the analysis are listed in table 
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the simulation normalized to input e-p cross section is then determined for that AP 
window cut. 

In addition, three coincidence kinematics were taken in order to allow for a sepa- 
ration of the elastic events from background events. Such separations helped us test 
our calculations of the lineshapes of the AP spectrum. The left arm was used to 
detect protons while electrons scattered in coincidence with protons were detected 
using the right arm spectrometer. The coincidence data were also used to provide 
a check on the scattering kinematics and to measure the proton detection efficiency 
and absorption. 

3.2 Why Detect Protons? 

The majority of the previous elastic e-p cross section measurements were made by 
detecting electrons rather than protons. However, proton detection has several ad- 
vantages over electron detection: 

1. Protons at moderately large angles correspond to electrons scattered at small 
angles. Detecting protons allows us to go to lower values of electron scattering 
angles (down to ~ T) than would normally be possible. 

2. Detecting protons reduces the variation of the cross section with the scattering 
angle. The cross section for the forward angle electrons varies rapidly with the 
scattering angle, while such variation for the corresponding protons is much 
smaller. The reverse is true for the backwards angle electrons where the varia- 
tion of the cross section with the scattering angle is greater for the corresponding 
forward angle protons, however, it is still a smaller effect than that for the for- 
ward angle electrons. 

3. Detecting protons reduces the variation of the cross section with the beam energy. 
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4- Detecting backwards angle electrons results in a reduced cross section within the 
angular acceptance of the spectrometer. On the other hand, the corresponding 
protons fall within a narrow angular window resulting in higher counting rates 
in less running time. That is, electrons are cross section limited at large angles 
(small e), while protons cross section is 10-20 times larger for small e. 

5. The proton momentum is constant for all e values for a given Q 2 . 

6. The linear e dependence of the radiative corrections is smaller for the protons 
than for the electrons. 

Because of these advantages, detecting protons greatly reduces the e-dependent 
systematic corrections and associated uncertainties applied to the measured cross 
sections as compared to detecting electrons. Figure E21 shows the momentum of the 
proton and electron as a function of the virtual photon polarization parameter e for 
Q 2 = 2.64 GeV 2 . The momentum is constant for the proton, while it varies by a factor 
of ~ 20 for the electron. The fact that the proton momentum is the same for all e 
values at a given Q 2 point means that there is no e dependence due to any momentum- 
dependent corrections due to detector efficiency, particle identification, and multiple 
scattering. It should be mentioned that any momentum- dependent correction will 
introduce an uncertainty in the reduced cross sections at a given Q 2 point which in 
turn introduces an uncertainty in both Ge p and Gm p but not in the ratio. 

Figure shows the cross section for both the proton, J^, and electron, J^, 
as a function of e for Q 2 = 2.64 GeV 2 . The cross section, and thus rate for fixed 
conditions, is nearly constant for the proton and that will reduce the effect of rate- 
dependent uncertainties dramatically. Also, the low e cross sections are no longer 
rate limited for the proton as it is the case for the electron which will require more 
running time for comparable statistics. 
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Figure 3.2: Plot of the proton and electron momentum as a function of e at Q 2 = 
2.64 GeV 2 . 




Figure 3.3: Plot of the proton and electron cross section (a x 10 35 ) (cm 2 ) as a 
function of s at Q 2 = 2.64 GeV 2 . 
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Figures 13.41 and 13.51 show the sensitivity of the cross section to the beam energy 
(percentage change in the cross section per one percent change in the beam energy) 
and scattering angle (percentage change in the cross section per one degree change in 
the scattering angle) for both particles as a function of e for Q 2 = 2.64 GeV 2 . The 
variation of the cross section with beam energy and especially with scattering angle 
as a function of e is less for the proton than that of the electron. 

Figure 13.41 shows the sensitivity of the electron and proton cross section to a 
1% change in the beam energy. For the electron, the maximum variation of the cross 
section is ~ 7% and for the proton it is 4%. If a 0.5% measurement in the cross section 
is desired, then the beam energy must be known with a precision of 7.14xl0 -4 for 
the electron and 1.25 xlO -3 for the proton. Such high precision is achievable for both 
particles knowing that energy measurements at JLAB can performed with precision 
as high as ^ = 2xl0~ 4 . 

Looking at Figure l3~B1 it can be seen that the biggest variation of the cross section 
of the electron for a one degree change in the scattering angle is ~ 55%. If we want 
to achieve a 0.5% measurement in the cross section, we need to know the scattering 
angle to within ~ 0.158 mrad which is very difficult to achieve. On the other hand, for 
the proton, over the same e range, the variation is a factor of 3 less and the scattering 
angle need only to be known to within ~ 0.476 mrad which is much more readily 
achievable. 

Finally, Figure l3~o1 shows the radiative correction factor (internal corrections only) 
as a function of e for Q 2 = 2.64 GeV 2 based on calculations done by Afanasev et al |49j . 
While the magnitude of the corrections is similar and both show an approximately 
linear dependence on e the dependence is much smaller for protons, ~ —8%, than for 
the electron ~ 17%. Note that the radiative corrections have an e dependence that is 
comparable to the slope brought about by the form factors. It is therefore extremely 
important that the radiative corrections (section 15 .6|) are correctly handled. 
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Figure 3.4: Plot of the proton and electron cross section sensitivity to the electron 
beam energy offset (%/%) as a function of e at Q 2 = 2.64 GeV 2 . 
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Figure 3.5: Plot of the proton and electron cross section sensitivity to the scattering 
ane le offset {%) (%/degree) as a function of e at Q 2 = 2.64 GeV 2 . 
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Figure 3.6: Plot of the proton and electron radiative correction factor (internal cor- 
rections only) function of e at Q 2 = 2.64 GeV 2 . 

3.3 The Continuous Electron Beam Accelerator 
Facility (CEBAF) 

During the run of the E01-001 experiment, the Continuous Electron Beam Acceler- 
ator Facility of the Thomas Jefferson National Accelerator Facility provided an 
unpolarized electron beam in the range of 1.912< E <4.702 GeV with beam currents 
up to 70 /j,A. Figure I3.7I shows the layout of the Jefferson Lab accelerator. 

Using a state-of-the-art strained GaAs photocathode gun system with maximum 
current of (J ~ 200 //A) and polarization above 70%, continuous-wave (CW) beams 
of high current and polarization are delivered to both Hall A and Hall C (J ~ 100 
//A). Meanwhile, a high polarization and low current beam is delivered to Hall B (/ ~ 
100 nA). 

First, electrons from the photocathode gun is accelerated to 50 MeV, and then 
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Figure 3.7: Schematic layout of the Thomas Jefferson National Accelerator Facility. 

injected into the north-linac. The north-linac consists of 20 Radio Frequency (RF) 
cryomodules. Each cryomodule has eight accelerating superconducting niobium cavi- 
ties kept at a temperature of 2 K using liquid helium coolant form the Central Helium 
Liquefier. By the time the electrons reach the end of the north-linac, they could have 
been accelerated up to 600 MeV by the 160 cavities. At the end of the north-linac, 
180° bending arcs (east arc) with a radius of 80 meters join the north-linac to the 
identical and antiparallel superconducting south-linac forming a recirculating beam- 
line. The beam through each arc is focused and steered using quadrupole and dipole 
magnets located inside each arc. The east arc has a total of 5 arcs on top of each 
other each with different bending field. The beam is steered through the east arc and 
passed on to the south-linac where it gets accelerated again and gain up to 600 MeV. 
At the end of the south-linac, the beam can be sent through the west arc for another 
pass or it can be sent to the Beam Switch Yard with a microstructure that consists 
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of short (1.67 ps) bursts of beam coming at 1497 MHz. Each hall receives one third 
of these bursts, giving a pulse train of 499 MHz in each hall. 

If another beam pass is desired for higher energy, the beam can be sent through the 
west recirculating linac for additional acceleration in the linacs, up to 5 passes through 
the accelerator. At the west recirculating linac, there are 4 different arcs on top of 
each other each with different bending field. It must be mentioned that the energy 
of the extracted beam is always a multiple of the combined linac energies Ei inacs = 
Enorth iinac + E south Hnac , plus the initial injector energy E injector or {E injector + nE Unacs ) 
where n is the number of passes. 

3.4 Hall A Beam Energy Measurements 

Precise measurements of the beam energy are required in order to extract the form 
factors of the protons from the elastic e-p cross sections. There are two different 
measurements that can be performed to determine the energy of incident electron 
beam with precision as high as = 2x 10 4 . These measurements are known as the 
arc and ep measurements. 

3.4.1 The Arc Beam Energy Measurement 

In the arc measurement jHSl EH EHj; the momentum p of the incident electrons is 
determined by knowing the net bend angle 9 of the electrons and the integral of the 
magnetic field along the arc section (electron path) of the beam line: 




where c is the speed of light. Equation ()3.3|) is derived based on the fact that when 
an electron moves in a circular motion with velocity v in region of magnetic field B, 
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— * 2 

where the magnetic force on the electron ^v~xB provides the central force v^-. Here, 
e is the charge of the electron, m is the mass of the electron, and r is the radius 
of circular path which is related to the net bend angle and linear electron path I as 
I = r6. 

In the arc measurements, simultaneous measurements of the bend angle and the 
magnetic field integral of the 8 bending dipoles (based on the measurement of the 9 th 
dipole as a reference) in the arc section of the beamline are made. The bend angle is 
determined by measuring the beam position and profile at the entrance and exist of 
the arc using four wire scanners or superharps. During the bend angle measurement, 
the quadruples are turned off (dispersive mode). The nominal bend angle of the beam 
in the arc section of the beamline is <ft — 34.3°. 

A superharp consists of three wires, two vertical wires for the horizontal beam 
profile measurement and one horizontal wire for the vertical beam profile. The signal 
from the wires gets picked up by an analog-to-digital converter (ADC). In addition, 
a position encoder measures the position of the ladder as the wires pass through 
the beam. The signal from the ADC and the position of the ladder determines the 
position and profile of the beam. Figures and 1331 show the superharp system and 
the arc section of the beamline. 



Preamp 




ADC 



Computer 
Readout 



Figure 3.8: Schematic of the superharp system. 
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Figure 3.9: Schematic Layout of the Arc Section of the Beamline. 
3.4.2 The ep Beam Energy Measurement 

The second method of measuring the beam energy is the ep measurement [HU EH1 
155] . A stand-alone device along the beamline located 17 m upstream of the target 
is used for the ep measurement. Based on the two-body kinematics in the elastic 
1 if(e, e'p) reaction, an incident electron is scattered by a CH 2 film target enclosed by 
an aluminum cover. The scattered electron and recoiled proton are detected using 
two identical arms, each of which has a detector package as shown in Figure 13.101 
The angles of the scattered electron 9 e and recoiled proton 9 P are measured, and the 
incident electron energy E Q is then determined from the two-body kinematic equation: 



E = ^ \ ^ + , (3.4) 



cos# e + (sin0 e /tan# p ) — 1\ n( m e 
1 - cos9 e ) + ^~E 

where M p is the mass of the proton, m e is the mass of the electron, and E is the 

final energy of the scattered electron. The energy measurements by the arc and 

ep methods show an excellent agreement with each other within an experimental 

uncertainty of < 3 x 10~ 4 . During the run of the E01-001 experiment, three arc and 
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two ep measurements were made to determine the beam energy. 




Figure 3.10: Schematic Layout of the ep energy measuring system. 

Table l3~2l lists all of the beam energy measurements for the E01-001 experiment. 
The Tiefenback energy measurement is similar to the arc measurement but is deter- 
mined using Hall A arc beam position monitors (BPMs) instead of the superharps to 
calculate the beam energy. The Tiefenback value for the beam energy, which is con- 
sistent with the arc measurements, is used in the analysis of the E0 1-001 experiment. 

Random (point-to-point) and scale uncertainties of 0.01% and 0.05% have been 
reported on the non-invasive arc measurements (the Tiefenback energy), respectively, 
and a 0.02% random uncertainty in the full invasive arc measurement or the ep 
measurement [U3]. Since the Tiefenback results were consistent with the full arc and 
ep measurements where we took them, we can assume that the absolute uncertainty 
in the Tiefenback is closer to the 0.02%. Therefore, an overall scale and random 
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Pass 


(Tiefenback) 


AE 


Arc 


ep 


Arc /Tiefenback 


ep / Arc 




(MeV) 


(MeV) 


(MeV) 


(MeV) 






2 


1912.94 


0.69 










2 


2260.00 


0.81 


2260.20 


2260.83 


1.000088 


1.000279 


3 


2844.71 


1.03 










4 


3772.80 


1.36 


3773.10 


3775.23 


1.000080 


1.000565 


5 


4702.52 


1.70 










5 


5554.60 


2.00 


5555.17 




1.000103 





Table 3.2: Tiefenback, Tiefenback quoted uncertainty AE, arc, and ep beam energy 
measurements of the E0 1-001 experiment. Tiefenback energy was used for the analysis 
with final uncertainty of 0.03% offset combined with 0.02% point-to-point uncertainty. 



uncertainties of 0.03% and 0.02% are estimated on Tiefenback energy, respectively. 

The classification of the uncertainty into a slope, random, and scale will be dis- 
cussed in detail in section I4~T1 The sensitivity of the cross sections to a 0.02% beam 
energy offset has been studied. Figures 13.111 and 13.121 show the relative difference 
between the nominal cross sections, determined at the nominal scattering angle and 
energy, and cross sections with a 0.02% shift in the beam energy for the right and 
left arms and at all kinematics. A 0.02% energy fluctuation changes the cross section 
by 0.01-0.03% for the right arm and 0.04-0.08% for the left arm. Such change in the 
cross sections was applied as a random uncertainty to each e point. 

Similarly, the sensitivity of the cross sections to an overall 0.03% beam energy 
offset has been studied. Figures 13.131 and 13.141 show the relative difference between 
the nominal cross sections and cross sections with a 0.03% beam energy offset for the 
right and left arms and at all kinematics. In order to determine the slope uncertainty 
in the cross sections due to a 0.03% beam energy offset, a fit of the relative difference 
in the cross sections to a straight line was performed at each Q 2 value. The average 
slope was used as the overall slope uncertainty. On the other hand, the deviation of 
the relative difference in the cross sections from the data and fit was then used as a 



69 

measure of the random uncertainty in each e point, while the average of the relative 
difference in the cross sections was used as an overall scale uncertainty. The right arm 
results show an average scale, random, and slope uncertainties of 0.034%, 0.01%, and 
0.29%, respectively. The left arm results show an average scale, random, and slope 
uncertainties of 0.13%, 0.02%, and 0.073%, respectively It should be mentioned that 
the same procedure of estimating the scale, random, and slope uncertainties in the 
sensitivity of the cross sections to a 0.10 mrad and 0.18 mrad angle offset will be used 
in section [Ol 
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Figure 3.11: The relative difference between the right arm nominal cross sections and 
cross sections with a 0.02% beam energy offset (red squares) as a function of e at 
all 5 incident energies. The black squares are the nominal cross sections relative to 
themselves. 
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Figure 3.12: The relative difference between the left arm nominal cross sections and 
cross sections with a 0.02% beam energy offset (red squares) as a function of e. The 
black squares are the nominal cross sections relative to themselves. 




Figure 3.13: The relative difference between the right arm nominal cross sections and 
cross sections with a 0.03% beam energy offset (red squares) as a function of e at all 
5 incident energies. The solid blue line is a linear fit to the data. 
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Figure 3.14: The relative difference between the left arm nominal cross sections and 
cross sections with a 0.03% beam energy offset (red squares) as a function of e. The 
solid blue line is a linear fit to the data. 

3.5 Beam Position Measurements 

Knowing the position and direction of the beam on the target is crucial since any 
beam offset could be translated into an uncertainty in the effective target length used 
in the analysis. During the E01-001 experiment, the beam first was rastered using a 
rectangular raster producing a ~ 2 mm x 2 mm spot size at the target to prevent 
any damage to the target by overheating it or local density fluctuations. 

In order to determine the beam position and direction at the target, two beam 
position monitors (BPMs), BPMA and BPMB, [HH E3 E3 located 7.516 and 2.378 
m upstream from the target were used. The BPM is a cavity with 4-wire antenna 
in one plane with frequency tuned to match the RF frequency of the beam (1497 
MHz). The absolute position of the beam is determined by the BPMs by calibrating 
the BPMs with respect to wire scanners (superharps) located adjacent to each BPM 
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at 7.345 and 2.214 m upstream of the target. 

The standard difference-over-sum method is used to determine the relative posi- 
tion of the beam to the wires. The position of the beam averaged over 0.3 sec, as 
read from the BPMs, is written into a data stream system on a event-by-event basis. 
The actual beam position and angle at the target can be reconstructed as: 

%beam I bpma^bpmb •^bpmb^bpma 

(3.5) 

a 3 \ / 



Vbeam ^ ( ybpma^bpmb Vbpmb^bpma 

(3.6) 



@beam — > (3-7) 

a 3 



, _ a 2 

'beam 



(3.8) 



where (L\ Xjjpmb Xfjpmaj &2 Vbpmb 1/bpmai and GI3 Zbpmb ^bprna- Also, X\,pmai 

Vbpma, Xbpmb, and ybpmb are the x and y coordinates of the beam as determined by 
beam position monitor a and b (BPMA and BPMB). Finally, %, ma and Zb pm b are the 
locations of the superharps for both BPMA and BPMB and have a value of -2.214 
and -7.345 m, respectively. The values of z bpma and Zb pm b are negative because the 
superharps are located upstream with respect to the target. 

Using the above equations, the position of the beam on the target was recon- 
structed. Figure 13.151 top (bottom) shows the x(y) coordinate of the beam at the 
target for all the runs in order of increasing e and Q 2 . For example, the point Q 2 
= 2.64 GeV 2 has 5 e points and a different color is given for each e value starting 
with black for the lowest e point and then red for the next and higher e point and 
so on. See Table 13.11 for kinematics description. The beam was well focused on the 
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target with an average (x,y) position of (-0.20,-0.10) mm, and an average beam drift 
of ±0.30 mm. The position uncertainties of 0.30 mm in two BPMs a ~ 5 m apart 
yield a ~ 0.07 mrad angle uncertainty in the beam angle. 
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Figure 3.15: The top(bottom) plot shows the x(y) coordinate of the beam at the 
target. 



3.6 Beam Current Measurements 

In order to measure the current used and therefore the total charge accumulated 
during each run in each kinematics, two identical beam current monitors (BCMs) 
located ~ 25 m upstream from the target are used. The two BCMs are calibrated at 
several beam currents relative to a parametric current transformer or Unser monitor. 
The Unser monitor is located halfway between the two BCMs and is calibrated by 
passing a precisely known current through it. Hall A beam current monitors jM] are 
designed to provide a non-invasive continuous measurements for the beam current. 
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The Hall A BCM is a stainless steel cylindrical cavity 15.24 cm in length and 15.48 
cm in diameter. The cylindrical axis of the cavity coincides with the beam direction. 
When the electron beam enters the cavity, it excites the transverse magnetic mode 
of the cavity. The resonant frequency of the cavity is tuned to the frequency of the 
electron beam (1497 MHz) and hence an output signal from one of the two antennas 
inside the cavity is produced. The output signal is then amplified and split into 
two components, one of which is sent to an AC Voltmeter which measures the beam 
current averaged over 1.0 sec period. The second signal is then converted into an 
analog DC voltage level by using an RMS-to-DC converter. The analog DC voltage 
level is then converted to a frequency signal using a voltage-to-frequency converter. 
Scalers gated by the start and end of each run receive the frequency signal and then 
provide a measurement of the accumulated charge during the runs. Figure l3~. 1 61 shows 
the schematic layout of the BCM readout block. 
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Figure 3.16: Schematic layout of the Hall A Beam Current Monitors readout block. 
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3.7 Hall A Target System 

The target system of Hall A [BUIEE] consists of the scattering vacuum chamber and the 
cryogenic targets. The 4-cm long liquid hydrogen target was used in this experiment. 

3.7.1 The Scattering Chamber 

The scattering vacuum chamber of Hall A consists of three main sections or rings. The 
first section, the base ring, is fixed on the pivot of the hall and made of stainless-steel. 
It contains a vacuum pump-out port and several viewing and electrical ports. The 
middle section, middle ring, is located at the beam height and is made of aluminum. 
It has a 104 cm inner diameter and 5 cm thick aluminum wall. Also a 15.2 cm vertical 
cutout on each of the beam side over the whole angular range of 12.5°< 9 <167.5° 
has been introduced. A 0.38 mm thin aluminum foil covers the vertical cutout on 
both side of the beam. The middle section has beam entrance and exit ports which 
are vacuum coupled to the electron beamline to prevent the beam from interacting 
with any material except the target. The third section, the upper ring, houses the 
cryogenic targets system. 

3.7.2 The Cryogenic Target 

The cryogenic target system of Hall A is mounted on a ladder inside the scattering 
vacuum chamber. The ladder contains sub-systems for cooling, gas handling, tem- 
perature and pressure monitoring, and target control and motion. The ladder also 
contains a selection of solid targets such as dummy target for background measure- 
ments and BeO, 12 C, and optics targets for beam viewing and calibration. The desired 
target can be selected from the control room (counting house) by moving the ladder 
vertically up and down until the target is aligned with the beam. 

The cryogenic target has three independent loops. Two loops were configured to 
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hold liquid hydrogen (LH 2 ) or liquid deuterium (LD 2 ), and one for a gaseous helium. 
Fans are used to circulate the liquid or the gas through each loop. Each of the LH 2 
and LD 2 loops has two aluminum cylindrical target cells of either 4 or 15 cm in length 
and 6.35 cm in diameter. The sidewalls of the cells are 0.178 mm thick with entrance 
(upstream) and exit (downstream) windows of 0.071 and 0.102 mm thick, respectively. 

During the E01-001 experiment, the 4-cm LH 2 target was used and operated at a 
constant temperature of 19 K and pressure of 25 psi with a density of 0.0723 g/cm 3 . 
The temperature of the target was stabilized using a high-power heater to compensate 
for the effect of any beam intensity variation on the target's temperature. The target 
was cooled using a target coolant (liquid helium) at 15 K supplied by the End Station 
Refrigerator (ESR). Because of the small spot size of the beam, the beam was rastered 
using a rectangular raster producing a ~ 2 mm x 2 mm spot size at the target to 
prevent any damage to the target by overheating it. 

3.8 Hall A High Resolution Spectrometers 

Protons were detected simultaneously using the two identical high resolution spec- 
trometers (HRS), also called the left and right arm spectrometers of Hall A. The left 
arm spectrometer was used to measure three Q 2 points of 2.64, 3.20, and 4.10 GeV 2 . 
Simultaneously, measurements at Q 2 = 0.5 GeV 2 were carried out using the right arm 
spectrometer which served as a luminosity monitor to remove any uncertainties due 
to beam charge, current, and target density fluctuations. 

The two spectrometers are identical in their design jM]. They can provide a 
maximum central momentum of ~ 4 GeV/c with momentum resolution better than 
2x 10~ 4 and a horizontal angular resolution better than 2 mrad. Figure 13". 1 71 shows the 
basic layout of the high resolution spectrometer, as well as the detector hut which will 
be discussed later on in more detail. Table l3~3~l lists some of the characteristics of the 
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Figure 3.17: Schematic layout of the Hall A high resolution spectrometer (HRS) and 
the detector hut. Dimension are in meters. 



Configuration 


QQDQ 


Momentum Range p (GeV/c) 


0.30 - 4.0 


Bend Angle (°) 


45 


Optical Length (m) 


23.4 


Momentum Acceptance 8p/p (%) 


± 4.50 


Dispersion (D) (cm/%) 


12.4 


Radial Linear Magnification (M) 


2.5 


D/M 


5.0 


Momentum Resolution (FWHM) 5p/p 


lxl0~ 4 


Angular Acceptance (Horizontal) (mrad) 


± 28 


Angular Acceptance (Vertical) (mrad) 


± 60 


Solid Angle AQ (msr) 


~ 6.7 


Angular Resolution (FWHM) Horizontal (mrad) 


0.6 


Angular Resolution (FWHM) Vertical 9 (mrad) 


2.0 


Transverse Length Acceptance (cm) 


± 5.0 


Transverse Position Resolution (FWHM) (mm) 


1.50 


Spectrometer Angle Determination Accuracy (mrad) 


0.10 



Table 3.3: Characteristics of the Hall A high resolution spectrometer. 
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Hall A high resolution spectrometer. A detailed description of the spectrometer design 
can be found in [fJDIHH] and references therein. Each high resolution spectrometer uses 
a QQDQ configuration of super- conducting magnets to focus charged particles onto 
their focal planes. Each spectrometer consists of two super-conducting quadrupoles 
followed by a 6.6 m long indexed dipole magnet for bending with focusing entrance 
and exit windows. Due to the trapezoidal cross sectional shape of the dipole, the 
bending field inside the dipole is not radially uniform. The first quadrupole is used 
to focus in the vertical plane, while the second and the third quadrupoles are used to 
focus in the horizontal plane. The net effect of the QQDQ configuration is to provide 
a vertical bending with nominal bending angle of 45°. 

3.9 Detector Package 

The detector package of each HRS is located in a large steel and concrete detector 
hut after the magnet system. The purpose of these detectors is to select and identify 
the charged particles coming through the spectrometers. During the E01-001 experi- 
ment, the two spectrometers (arms) included similar detector packages with a slight 
difference between the two arms. A detailed description of the detector package on 
each arm can be found in |M1 ITU] and references therein. The following detectors 
were used on each arm: 

• A set of two Vertical Drift Chambers (VDC's) for tracking purposes. 

• Two scintillator planes, Si and S2, for trigger activation and time-of- 
flight determination. In addition, an extra plane So was added to the 
left arm for efficiency determination. 

• Aerogel Cerenkov detectors A\ for the left arm and A 2 for the right 
arm for particle identification (p and 7r + separation). 
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• A gas Cerenkov detector for particle identification. The right arm 
gas Cerenkov detector was used during the coincidence runs for e~ 
and 7r separation. 

In addition, the right and left arm detector packages included some extra detectors 
that were not used during the E01-001. The right arm had a mirror aerogel detector 
and calorimeter. The left arm included a gas Cerenkov, pion rejector (two layers), 
and focal plane polarimeter (FPP) (rear straw chambers). Although these extra 
detectors were not used during the EO 1-001 experiment, they do contribute to proton 
absorption (absorption in the material of each detector) if they were located before 
the scintillator plane 5*2. The mirror aerogel detector and the gas Cerenkov detector 
are the most relevant in this case and their contribution to proton absorption should 
be taken into an account. On the other hand, both layers of the pion rejector and the 
two rear straw chambers of the FPP were located behind the 5*2 scintillator and were 
not relevant for the proton absorption determination. Figures 13.181 and 13.191 show 
a schematic layout of the detector package used with each arm during the EO 1-001 
experiment. 

3.9.1 Vertical Drift Chambers 

Vertical drift chambers (VDCs) jHH E] are used for tracking (position and slope 
of the particle trajectory) of the scattered particles. Each spectrometer of Hall A 
has two vertical drift chambers. The first VDC is located at the focal plane of the 
spectrometer, and the second VDC is positioned parallel to and 23 cm away from 
the first VDC. Both VDCs intersect the spectrometer central ray at an angle of 45° 
which is the approximate nominal angle at which the particle trajectory crosses the 
wire planes of each VDC. Each VDC has two wire planes known as the U and V 
planes. Each wire plane contains 368 parallel 20/xm-diameter gold-coated tungsten 
wires sandwiched in between gold-coated mylar planes. The wires of the U plane are 
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FPP Rear Straw Chambers 




Figure 3.18: Schematic layout of the left arm detector package used during the E01- 
001 experiment. 




Figure 3.19: Schematic layout of the right arm detector package used during the 
E01-001 experiment. 
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perpendicular to that of the V plane and they are oriented with an angle of 45° (-45°) 
with respect to the dispersive (transverse) direction. We refer to the U and V planes 
as Ui and Vi planes in the first VDC and U2 and V2 planes in the second. Figures 
13.201 and 13.211 show a pair of vertical drift chambers and a side view of a nominal 
45° particle trajectory crosses the planes of the VDCs. To operate the VDCs, a gas 



Figure 3.20: Schematic layout of the Hall A vertical drift chambers. 

mixture of argon (62%) and ethane (38%) (C2H 6 ) fills the area between the mylar 
planes. The mylar planes are kept at negative potential of 4KV, while the tungsten 
wires are grounded. When a charged particle crosses the VDC planes, it ionizes the 
atoms in the gas mixture. The released electrons are accelerated (drift) due to the 
potential difference between the mylar planes and the tungsten wires taking the path 
of least time (also known as the geodetic path). Once the electrons are close to the 
sense wires where the electric field is the strongest, the drifting electrons create an 
electron avalanche. This avalanche of electrons hits the sense wire and then generates 
a signal which gets amplified, discriminated, and then sent to a multihit time-to- 
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Figure 3.21: Side view of the Hall A vertical drift chambers planes. The red arrow is 
the nominal 45° particle trajectory crosses the planes of the VDCs. 

digital converter (TDC) to measure the time for the least path defined above. The 
TDC is started by the signal from the sense wire and stopped by the event trigger 
supervisor. 

In order to get the tracking information or the drift distance of the charged particle 
from each fired wire, we use the drift velocity of the ionized electrons in the C2H6 
gas mixture (known to be 50 /zm/ns), and the output of the TDC signal (time). 
Having calculated the drift distance of the charged particle from each fired wire, 
the trajectory (track) of the charged particle can be reconstructed. A 45° nominal 
trajectory (track) typically triggers five wires. Figure 13.221 shows the reconstruction 
of the charged particle track in the VDC planes. 

3.9.2 Scintillators and Triggers 

For most of the data taking, proton singles events were detected. In addition, co- 
incidence runs at Q 2 = 2.64 and 4.10 GeV 2 were taken. Electrons in coincidence 
with protons were detected using the right arm spectrometer while the protons were 
detected using the left arm. Triggering of such events was made possible by the use 



83 



Cathode Plane 




Figure 3.22: The avalanche process in the VDC and the reconstruction of the charged 
particle track. 

of the Hall A scintillators and trigger system |M1 172*] . 

For each spectrometer in Hall A, there are two scintillators planes labeled Si and 
S2. The two planes are parallel to each other and both are perpendicular to the 
nominal central ray of the spectrometer. The two scintillators planes are ~ 2 meters 
apart. Each scintillator plane has 6 identical overlapping scintillators paddles made 
of thin plastic (0.5 cm thick BICON 408) to minimize hadron absorption. A photo- 
multiplier tube (PMT) is attached to each end of each scintillator paddle to collect 
the photons produced by particles passing through scintillator. We refer to these 
PMTs as the left and right PMT of the scintillator paddle. Figure 13.231 shows the 
scintillator configuration. 

The active area of Si is ~ 170 cm x 36 cm (30 cm x 36 cm for each paddle) and 
~ 220 cm x 60 cm (37 cm x 60 cm for each paddle) for S2. The time resolution 
for each scintillator plane is ~ 0.30 ns. For the E01-001 experiment, an additional 
scintillator counter , So, was added to the left arm trigger system for a more accurate 
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measurement of the efficiency. So was installed before the Si scintillator plane and 
it is a 1 cm thick scintillator paddle with an active area of ~ 190 cm x 40 cm. The 
So paddle has 2 PMTs labeled top and bottom. The analog signals from the PMTs 




Figure 3.23: Scintillator plane configuration. Each scintillator plane is perpendicular 
to the nominal central ray of the spectrometer. 

are first sent to a discriminator (LeCroy model 4413/200) providing both digitized 
and analog outputs. The analog signals are sent to ADCs and the digitized signal are 
split into three signals. One signal is sent to TDCs, the second is sent to a scalers 
gated by the start and the end of each run, and the third is sent to a logical unit 
called the AND unit to make a coincidence between the pairs of PMTs viewing the 
same paddle. For each of the spectrometers, the 12 outputs of the logical AND unit 
are fed into the Memory Lookup unit (MLU) which takes in a combination of logical 
signals at its input and gives out a combination of logical signals at its output. 

There are five basic triggers (event types) generated from the scintillators timing 
information. They are classified as main physics triggers and loose (auxiliary) physics 
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triggers. The conventional triggers types used for the left arm are the main singles 
trigger T 3 and the loose singles trigger T 4 . For the right arm, we refer to the main 
singles trigger by Ti and the loose singles trigger by T2. Finally, there is the coinci- 
dence trigger T5 which is the coincidence of Ti and T3 implying that the two events 
Ti and T 3 were produced at the target simultaneously. The loose triggers are usually 
used to determine the efficiencies of the main triggers and the scintillators planes in 
general. 

The output of each ADC is proportional to the number of photons produced inside 
the scintillator and in turn represents the ionization energy loss of the particle, or 
equivalently, the energy deposition in the scintillator. For low momenta, the heavier 
the particle, the more energy it deposits in the scintillator. The TDC signal provides 
timing information for the different type triggers used. By determining the time-of- 
flight of the particle between the scintillator planes form the TDC signals, and by 
knowing the distance between the two scintillator planes, the velocity of the particle 
v or (3 = v/c can be determined. Usually (3 is used for particle identification (PID) 
since different charged particles of a known momentum can be separated by knowing 
their velocity and hence (3. Figure 1.1241 shows the ADC and TDC spectrum from 
one of the PMTs from the second scintillator plane S2 for the low-momentum setting 
spectrometer (right arm). A (3 spectrum where different particles are clearly separated 
is also shown. 

The left and right arm MLUs were programmed to generate a logical signal which 
defines and assigns the trigger (event type) Ti(T 3 ) to an event (particle) if the event 
satisfies the following: 

• The main Ti(T 3 ) event type is generated if: 

1. A scintillator paddle is said to have fired if and only if a signal from both 
the right and left PMT in that paddle is generated. The logic AND is then 
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Figure 3.24: Typical ADC and TDC spectra from a scintillator PMT (upper two 
plots) from the low-momentum setting spectrometer. Different particles have been 
identified using the ADC spectrum based on the amount of the energy deposited in 
the scintillator. Pions show first (time wise) in the TDC spectrum since they are 
lighter than the protons and therefore have smaller time-of-flight (larger (3). The 
calculated /3 spectrum based on the TDC timing information and distance between 
the two scintillators planes is also shown (lower plot). The (3 spectrum is used for 
particle identification. 
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formed between the right and left PMT's signal to generate that paddle's 
signal. 

2. Each scintillator plane will have 6 signals (one signal from each paddle in 
each plane). 

3. The logic OR of the 6 signals in that plane generate the main signal for 
that plane (one main signal for each scintillator plane). 

4. The right arm main trigger (Ti) is defined as: 

(Si main signal) AND (S2 main signal). 

5. A third scintillator So was added to the left arm only. Therefore, the left 
arm main trigger (T3) is defined as: 



• The loose T 2 (T 4 ) event type is generated if: 

1. T 2 is the same as Ti defined above except the AND in point number 4 is 
replaced with an OR: 

(Si main signal) OR (S 2 main signal). 

2. T 4 is the same as T 3 defined above except the AND between Si and S 2 in 
point number 5 is replaced with an OR: 



• The coincidence trigger T5 is generated as the coincidence between Ti and T3. 
3.9.3 Gas Cerenkov 

Coincidence kinematics were taken at Q 2 = 2.64 and 4.10 GeV 2 . Electrons were 
detected using the right arm and protons using the left arm. See Table 13.11 for 
more details. In order to separate electrons from background particles such as the 



(Sq main signal) AND (Si main signal) AND (S 2 main signal). 



(Sq main signal) AND ( (Si main signal) OR (S 2 main signal) 



8. 



S 



negatively charged pions, a Cerenkov detector was used in the right arm for particle 
identification. 

The gas Cerenkov detector used in Hall A (HH EH] was installed between the 
scintillator Si and S2 planes. It is a rectangular chamber filled with CO2 gas at 
atmospheric pressure. When a high energy charged particle travels through a material 
of index of refraction n with velocity v larger than the velocity of light in that material 
(c/n) where c is the speed of light, electromagnetic radiation (light) is emitted. This 
is known as the Cerenkov effect. For the Hall A Cerenkov detector, 10 spherical 
mirrors at the back wall focus the Cerenkov radiation on 10 PMT photocathodes. 

A charged particle can emit Cerenkov radiation only when its velocity is larger 
than a threshold velocity set by the index of refraction of that material: 



For the C0 2 gas used, the index of refraction is n = 1.00041. Equation (|3.10j) 
then gives a momentum threshold of p t h ~ 17 MeV/c for electrons and p t h ~ 4.8 
GeV/c for pions which is larger than the momentum range of the Hall A HRS which 
is 0.30—4.0 GeV/c. So, only electrons give a signal in the analog-to-digital (ADC) of 
the Cerenkov detector. Figure 13*. 251 shows the distribution of the sum of the 10 ADCs 
signals of the gas Cerenkov PMTs. 




(3.9) 



and therefore a threshold momentum: 



mvth 



(3.10) 



Pth = 
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Figure 3.25: Gas Cerenkov ADCs sum signal. The peak at GAS ADCSUM = zero 
represents the ir~ which did not fire the gas Cerenkov. GAS ADCSUM > 200 repre- 
sents the e~ signal which fired the gas Cerenkov. The red solid line is the boundary 
between the vr" region (GAS ADCSUM < 200) and e~ region (GAS ADCSUM > 
200). 
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3.9.4 Aerogel Cerenkov 

Two Aerogel Cerenkov detectors Ai and A 2 [OH EH EE] were used for particle identi- 
fication where protons and positively charged pions could be separated based on the 
Cerenkov radiation effect discussed above. The A 2 aerogel detector was used with 
the right arm, while the Aj aerogel was used with the left arm. In addition, a mirror 
aerogel, Am, was installed on the right arm between the Si and S 2 scintillators but 
was never used. The Am Aerogel contributes to the proton absorption in the right 
arm. 

The two aerogel detectors Ai and A 2 are designed in the same way. The Ai aerogel 
has a 9 cm aerogel radiator with index of refraction ua 1 = 1.015 while A 2 has a 5 cm 
aerogel radiator with index of refraction tia 2 = 1.055. There are 24 PMTs installed 
in Ai and 26 PMTs in A 2 . 

For the Ai aerogel (n^ = 1.015), equation (J3.10j) gives a momentum threshold of 
Pth ~ 5.4 GeV/c for protons and p t h ~ 0.8 GeV/c for positive pions, while, for the 
A 2 aerogel (ua 1 = 1.055), equation ()3.1U|) gives a momentum threshold of p t h ~ 2.8 
GeV/c for protons and p t h ~ 0.40 GeV/c for positive pions. For the A x aerogel (left 
arm), the momentum threshold for the protons is larger even than the momentum 
range of the Hall A HRS (0.30—4.0 GeV/c) and that makes it impossible for the 
protons to fire the aerogel Cerenkov detector. On the other hand, the momentum 
threshold for the positive pions is within the momentum range of the Hall A HRS and 
hence pions above 0.8 GeV/c should produce a signal in the analog-to-digital (ADC) 
of the Ai aerogel Cerenkov detector. For the right arm where A 2 aerogel was used, 
only positive pions could trigger the analog-to-digital and fire the aerogel Cerenkov 
detector since the central momentum for the right arm spectrometer was always set 
to 0.756 GeV/c which is less than the threshold momentum needed for the protons 
to trigger the aerogel Cerenkov detector. Figure 13.261 shows the distribution of the 
sum of the ADCs signals of the A 2 (Ai) aerogel Cerenkov PMTs. 



91 



C 1n 5 

w 

10 4 



10 



10' 
10 











All A 2 ADCSUM Sp< 


jctrum 
















Protons 




















Deuterons 

A„ ADCSUM = 125( 


) 




























































^iP^M rinflnU n 





s 

10 



10 

10 : 

10 

1 



1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

A 2 ADCSUM 













A t ADCSUM > 


350 
















P = 


A 1 ADCSUM < 


350 


















A, ., 


ADCSUM = 350 




















































II n r 



200 



400 



600 800 1000 1200 1400 1600 1800 2000 



A, ADCSUM 



Figure 3.26: A 2 (A!) aerogel Cerenkov ADCs sum signal top(bottom) in blue. The 
peak at A 2 (A!) ADCSUM = zero represents the protons which did not fire the aerogel 
Cerenkov. The contribution of protons, 7r + , and deuterons to A 2 ADCSUM (selected 
by cuts on the time-of- flight and amount of energy deposited in the Si scintillator) is 
shown. The contribution of protons, tt + , and deuterons to A x ADCSUM is not shown 
since they cannot be separated using any other detector and it will be discussed in 
the next chapter in more detail. The A 2 (Ai) ADCSUM > 1250(350) represents the 
cut we use to separate the ir + signal which fired the aerogel Cerenkov. Notice the 
leakage of protons into the ir + area and the small contamination of the protons area 
with 7r + and deuterons as can be seen form A 2 ADCSUM spectrum. 
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3.10 Data Acquisition System 

The data acquisition system in Hall A (DAQ) was used during the run of experiment 
EO 1-001. The data acquisition system is controlled by CEBAF On-line Data Acquisi- 
tion system (CODA) [HUEE] which is developed by the JLAB data acquisition group 
and designed for nuclear physics experiments. CODA is a toolkit that is composed 
of a set of software and hardware packages from which DAQ can be built to manage 
the acquisition, monitoring, and storage of data. 

Data acquisition in Hall A takes the following steps: 

• The data is read out from Read-Out Controllers (ROCs). The ROCs are CPUs 
in Fastbus and VME crates in the hall and in the electronics room. These crates 
contain the ADCs, TDCs, and scalers that contain the event information. 

• The Trigger Supervisor controls the state of the run, and generates the triggers 
that cause the ROCs to be read out. 

• The Event Builder subsystem (EB) is the part of CODA that reads in the data 
fragments from the ROCs and puts the data together into an event, incorporat- 
ing all of the necessary CODA header information needed to describe and label 
the event and the data fragments. 

• CODA manages the data acquisition system, and takes care of handling the 
data from the events. 

• After the event is built by the EB, it is placed into a buffer, after which it can 
be tested and recorded using the event recorder or rejected if desired. 

• The accepted events (data) are written to a local disk and then transferred to 
the Mass Storage System. 
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• Data from various control systems and scalers are injected into the data stream 
every few seconds using the event transfer. 

• In addition to running the data acquisition, CODA also includes a graphical 
user interface (RunControl) which allows the user to start and stop runs, as 
well as define run parameters. 

• Events are classified as physics events which come from the spectrometers (de- 
tectors) and beamline information (Beam position monitors, beam loss moni- 
tors, and beam raster readback values recorded for each event), or EPICS events 
[77] such as the readout of hardware (spectrometers magnets settings, angles, 
and target controls). 



Chapter 4 Data Analysis I: Efficiencies 
and Corrections 



4.1 Analysis Introduction 

In this chapter a description of the event reconstruction procedure and the corrections 
and efficiencies applied to the measured beam charge, Q, will be presented. These 
corrections and efficiencies include the computer and electronics livetimes , CLT, and 
ELT, respectively, VDCs tracking efficiency, €vdc, VDCs hardware cuts efficiency 
€vdch, scintillators efficiency (product of the two scintillators efficiency or e Sl x es 2 ), 
particle identification (PID) efficiency, epm-, proton absorption correction, C Absorption, 
target boiling correction, Ctb, and target length correction, Ctl- In addition, a 
brief description of the spectrometer optics calibration and spectrometer mispointing 
measurements will be given. Having applied these corrections and efficiencies to the 
beam charge, we refer to the corrected charge as the effective charge or Q e ff and it 
is defined as: 



Qeff — ( Q x ELT X C LT X 6vDC x €VDCH x 6s 1 X €s 2 X €piD X C Absorption 



where ps is the prescale factor. The prescale factor n for the trigger type T, (i = 
1, • • • ,5) means that the Trigger Supervisor will read out every nth event of type Tj. 

The advantages of detecting protons as against electrons were discussed in section 
13.21 Because of these advantages, the E0I-00I experiment was able to greatly reduce 





(4.1) 
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any e-dependent systematic corrections and associated uncertainties applied to the 
measured cross sections. In addition, measurements at Q 2 = 0.5 GeV 2 using the right 
arm spectrometer which served as a luminosity monitor will check the uncertainties 
due to beam charge, current, and target density fluctuations. 

The systematic uncertainties in these efficiencies and corrections will be discussed. 
These systematic uncertainties are a measure of how accurately we know the various 
efficiencies and corrections. The systematic uncertainties are classified into three 
types, based on how they contribute to the quantities we want to extract: 

• Scale Uncertainty: Sometimes referred to as normalization uncertainty as it 
has the same effect on all e points at a given Q 2 . An example would be the 
uncertainty in the target length. Such uncertainty would affect both Ge p and 
Gm p the same way but not the ratio. 

• Random Uncertainty: Sometimes referred to as point-to-point uncertainty. It 
has an uncorrelated effect on each e point at a given Q 2 . Two examples are 
the statistical uncertainty and the uncorrelated shifts in the beam energy and 
scattering angle for different kinematics. Such uncertainties affect the extraction 
of Ge p , Gm p , and the ratio. 

• Slope Uncertainty: Refers to a correlated uncertainty in a correction that varies 
linearly with e. An example would be the effect of a fixed scattering angle or 
beam energy offset for all kinematics which also leads to a scale uncertainty. 
Such an uncertainty will modify the slope of the reduced cross section verses e, 
and thus Ge p and off-, but will not spoil the expected linearity. 

It must be mentioned that in some cases it is unclear whether the uncertainty is 
scale, random, or slope. In these cases, the worst case of the three uncertainties is 
assumed. Table 14.11 indicates whether the different types of uncertainties effect the 
individual form factors, their ratio, and the linearity of the L-T plots. Note that the 
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slope uncertainty is the change in slope between e = 0.0 and e = 1.0. So the value 
for the slope uncertainty is the change over the actual e range in data divided by 
Ae. Since Ae ~ 0.70 for the left arm and Ae ~ 0.07 for the right arm, the slope 
uncertainty is much larger for the right arm. 



Uncertainty Type 


Gep 


Gm p 


I^pGepI Gup 


Linearity 


Scale 


Yes 


Yes 


No 


No 


Random 


Yes 


Yes 


Yes 


Yes 


Slope 


Yes 


Yes 


Yes 


No 



Table 4.1: The effect of the scale, random, and scale uncertainties on the individual 
form factors, their ratio, and the linearity of the L-T plots. 



4.2 Event Reconstruction 

The raw data files imported from the DAQ system were replayed using the standard 
Hall A event processing software ESPACE (Event Scanning Program for Hall A Col- 
laboration Experiments) [7H]. The multipurpose ESPACE reads the raw events and 
then decodes, filters, calibrates, and reconstructs physical variables of interest for the 
analysis on an event-by-event basis. These physical variables can be detector hits, 
tracks and particle identification (PID) signals, or the coordinates of the reaction 
vertex in the target. The outputs of ESPACE are histograms and ntuples (multi- 
dimensional matrices of the physical quantities of interest) saved in HBOOK formated 
files [75] • An HBOOK file is generated for each production run. These HBOOKs are 
read using the Physics Analysis Workstation (PAW) software jHD] where cuts can be 
applied on these physical quantities. 
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4.3 Optics 

There are different coordinate systems used in Jefferson Lab Hall A. These coordinate 
systems are described in detail in [THJ |H] and can be classified as: the hall coordinate 
system, the target coordinate system (spectrometer reconstructed coordinate system), 
the detector coordinate system, the transport coordinate system, and finally the focal 
plane coordinate system. The hall and target coordinate systems are commonly used 
to define and calculate the physical quantities of interest for data analysis (see Figures 
14.11 and l4.2|) . while the detector, transport, and focal plane coordinate systems are 
used for Monte Carlo simulations and optics calibrations. 



X 



Beamline 



Y 



to beam dump 



Figure 4.1: The hall coordinate system viewed from above: Z is the beam direction, Y 
is vertically pointing up, and X is to the left of the beam direction and perpendicular 
to both Z and Y. 

For event reconstruction, ESPACE starts by reconstructing the event at the focal 
plane of the spectrometer (passing through the center of the Ui plane of the first 
VDC). That includes determination of the electrons drift times in the wire chambers 
and the drift distances (see section IH.9.1|) . reconstruction of the particle trajectories, 
and calculation of (3 = v/c and positions and angles of the track in the focal plane 
coordinate system. The detector coordinate system, the transport coordinate system, 
and the focal plane coordinate system all share the same origin defined as the inter- 
section of central wire 184 of the VDC U\ plane and the central wire projection of 
the VDC V 1 plane. 



98 



scattered Z (HRS 
electron central ray) 




X 

(vertical down) 



Figure 4.2: The target coordinate system: Z is the spectrometer central ray direction, 
X is pointing down (dispersive direction), and Y is perpendicular to both Z and X 
(transverse direction), and 9 (usually called (p tg and 9 tg where tg is short for target) 
are the in-plane and out-of-plane scattering angles, respectively, as measured with 
respect to the spectrometer central ray. 

Having reconstructed the focal plane variables, the reconstruction of the target 
variables or reconstruction of the trajectories at the target can be done by using the 
optics database or transformation matrix between the focal plane and target variables 
|HT] . These reconstructed target variables are the in-plane and out-of-plane scattering 
angles or <p tg , and 9 or 9 tg where tg is short for target, the y-coordinate of the 
extended target length y tg , and the deviation from the central momentum 5 = 2 ^ ££ - 

In order to study the optical properties of the spectrometer, sieve slits collimators 
[HUES positioned 1.184±0.005 and 1.176±0.005 m from the target on the left and 
right arm spectrometers, respectively, are used. The sieve slit is a 5 mm thick stainless 
steel sheet with a pattern of 49 holes (7x7), spaced 12.5 mm apart horizontally and 
25 mm apart vertically. Two of the holes, one in the center and one displaced two 
rows vertically and one horizontally, are 4 mm in diameter, while the rest are 2 mm 
in diameter. Figures 14.31 and 14.41 show the reconstruction of the sieve slit image 
due to electrons scattering from a thin 12 C target as seen by the right and left arm 
spectrometers, respectively. The fact that the optics database used by ESPACE 
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during the E01-001 analysis produces the correct position of the holes in the sieve 
slit for both spectrometers indicates that the optics database is well calibrated and 
no further calibration is needed. The calibration procedure is documented in detail 
in [EH- 
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Figure 4.3: The reconstruction of the sieve slit image into the right arm spectrometer 
due to electrons scattered from a thin 12 C target. 



During the analysis of the EO 1-001 experiment, a 1.6 msr solid angle software 
cut was applied. This cut corresponds to an out-of-plane angle of -40. 0< 6 tg <40.0 
mrad and in-plane angle of -10. 0< fa g <10.0 mrad. See section IS~2l and Table I5~2l 
for a full description of the cuts applied during the analysis. We estimate that we 
know the edge of the 9t g and fa g cuts from looking at the adjacent hole separation 
in the vertical and horizontal directions to within 0.20 mrad. This translates into 
a 2.0% scale uncertainty in the in-plane angle (0.2 mrad/10 mrad), and 0.5% scale 
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uncertainty in the out-of-plane angle (0.2 mrad/40 rnrad). The sum in quadrature 
of the two angle scale uncertainties was used to determine the final estimated scale 
uncertainty in the 1.6 msr solid angle cut of 2.06%. Because the 1.6 msr solid angle 
cut is identical for all kinematics, the uncertainty in the solid angle only contributes 
to the scale uncertainty. 
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Figure 4.4: The reconstruction of the sieve slit image into the left arm spectrometer 
due to electrons scattered from a thin 12 C target. 

4.4 Spectrometer Mispointing 

Experiment E01-001 requires a precise knowledge of the e-p reduced cross sections. 
In order to achieve that, an accurate knowledge of the scattering angle is required. 
Knowing the scattering angle accurately requires accurate knowledge of the spec- 



101 



trometer optics and offsets, target position, and beam position. Figure 14.21 shows the 
coordinate system of the spectrometer. Ideally, the target center coincides with the 
hall center. Due to translational movements of the spectrometer around the hall cen- 
ter, the central ray of the spectrometer can miss the hall center in both the horizontal 
and vertical directions. The horizontal offset or the horizontal distance between the 
hall center and the central ray of the spectrometer is referred to as the spectrometer 
mispointing. 

There are two different and reliable methods by which the horizontal and vertical 
offsets can be measured and hence the spectrometer angle: the survey method and 
the carbon-pointing method. The survey method is the most precise. The two spec- 
trometers used during the E0 1-001 experiment were surveyed at several kinematics 
settings by the JLAB survey group. 

For the right arm spectrometer, the spectrometer angle was surveyed and deter- 
mined at all 5 e points, and the survey angles were used in the analysis. For the left 
arm, there were several spectrometer settings where a survey was not performed, and 
so the carbon-pointing method was used. In the carbon pointing method, the spec- 
trometer mispointing, Ah, and spectrometer angle, 9 S , are determined using carbon- 
pointing runs where electrons are scattered by a thin carbon foil. By knowing the 
spectrometer central angle, o , (defined as the spectrometer angle with no mispoint- 
ing), as determined from the hall floor marks, target position as reconstructed by the 
spectrometer, y tg , and the target offset along the beam direction, z //> as measured 
by the target survey group, Ah can be determined as: 



and hence the spectrometer scattering angle (angle setting) can be calculated as: 



Ah = ±y tg + z of f sm(9 ) , 



(4.2) 




(4.3) 
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where the plus(minus) sign in front of y tg in equation (|4.2|) is used with the right (left) 
arm and L is the distance between the hall center and the floor marks where the 
angles are scripted and has a value of 8.458 m. Note that ^ = A0 O in equation (|4.2|) 
above represents the correction to the central scattering angle of the spectrometer, 
9 S = 9 Q ± A9 a . The spectrometer is said to be mispointed downstream (upstream) if 
Ah is positive(negative). 

We can test the nominal kinematics by looking at the reconstructed kinematics 
for elastic scattering. The reconstruction of the AP spectrum is discussed in section 
15.21 The elastic peak position should be near AP = 0.0 MeV, but there are small 
corrections due to energy loss and radiative corrections, which are modeled in the 
Monte Carlo simulation program SIMC. See section l5.5.1l for details. The elastic peak 
position in the AP spectrum from data, 5p RIGHT (LEFT), is compared to that in 
simulations Sp SIMC. This comparison was done for the right and left arms and at all 
kinematics. For the left arm, an overall angular offset of 0.19 mrad was applied to the 
pointing angles to best center the elastic peak position from data to that of simulations 
at each kinematics. Note that an offset of 0.28 mrad is needed to match the carbon 
angles to the survey ones. The two offsets are consistent taking into account the 
0.18 mrad scale uncertainty assigned in the scattering angle (discussed below). For 
the right arm, the survey angles are used and yielded a good AP peak position. 
Therefore, no additional offset was needed. Figures l4~51 and l4~o1 show the difference in 
the elastic peak position from data and that of simulations after applying the angular 
offset and for both arms. The error bars assume random uncertainties of 0.10 mrad 
for the angle (see discussion below) and 0.02% for the beam energy (see section EO)) . 
With these uncertainties, the values of 5p from data are in good agreement with those 
of simulations. In fact, the difference is zero to better than the error bars, indicating 
that we might be slightly overestimating the random uncertainties, but we will stick 
with these uncertainties to be conservative. 
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Figure 4.5: The right arm difference in the elastic peak position from data (5p 
RIGHT) and that of simulations (5p SIMC) after applying the angular offset. The 
points are sorted according to the kinematics of the left arm. Kinematics 1-5 cor- 
respond to Q 2 = 2.64 GeV 2 , while kinematics 6-9 (10-12) correspond to 3.20 (4.10) 
GeV 2 . For each Q 2 value, the points are sorted by e (low to high). Note that for the 
right arm only the first five measurements are truly independent. The later points 
are repeats at the same kinematics. 
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Figure 4.6: The left arm difference in the elastic peak position from data (dp LEFT) 
and that of simulations (Sp SIMC) after applying the angular offset. Kinematics 1-5 
correspond to Q 2 = 2.64 GeV 2 , while kinematics 6-9 (10-12) correspond to 3.20 (4.10) 
GeV 2 . For each Q 2 value, the points are sorted by e (low to high). 
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A 0.10 mrad random uncertainty in the angular offset is assigned. This is based 
on the following contributions combined in quadrature: 0.07 mrad due to drifts in the 
beam angle as determined from the BPMs, ~ 0.05 mrad uncertainty that comes from 
the ~ 0.1 MeV uncertainties in determining the peak position, and the ~ 0.05 mrad 
uncertainty from the pointing determination (run to run scatter and uncertainty in 
determining the target position). Table 14.21 shows the Tiefenback energy and the 
spectrometer settings used for the left and right arms. 



Kinematics 


Tiefenback Energy 


6 L 


Or 




(MeV) 


(°) 


(°) 





1912.94 


12.6311 


58.3090 


a 


2260.00 


22.1592 


60.0700 


i 


2844.71 


29.4590 


62.0380 


Q 


3772.80 


35.1512 


63.8710 


I 


4702.52 


38.2512 


64.9810 


b 


2260.00 


12.5226 


60.0700 


3 


2844.71 


23.3896 


62.0380 


V 


3772.80 


30.4802 


63.8710 


m 


4702.52 


34.1225 


64.9810 


k 


2844.71 


12.6807 


62.0380 


r 


3772.80 


23.6586 


63.8710 


n 


4702.52 


28.3735 


64.9810 



Table 4.2: The Tiefenback energy and spectrometer settings used in the actual analy- 
sis of the E0 1-001 kinematics. The right arm spectrometer settings Or are the survey 
angles, while the left arm spectrometer settings 6l are the angles as determined by 
carbon-pointing runs measurements corrected by 0.19 mrad to center the elastic peak. 



A 0.18 mrad scale uncertainty in the angular offset is assigned. The scale uncer- 
tainty is determined based on the following contributions combined in quadrature: 
0.13 mrad uncertainty due to uncertainty in modeling the energy loss and radiative 
effects on the elastic peak (if we shift the peak centers by 0.2 MeV due to error in 
SIMC energy loss, or small errors in the smearing we use, the angle offset varies by 
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0.13 mrad), 0.07 mrad uncertainty for possible beam angle offset, and 0.10 mrad 
uncertainty coming from the fact that a small shift in the beam energy can give a 
solution with the wrong angle offset (if we shift the beam energy by 0.03%, we still 
get decent elastic peak positions by shifting the angle by 0.1 mrad. The \ 2 is worse, 
but not totally unreasonable). 
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Figure 4.7: The relative difference between the right arm nominal cross sections 
and cross sections with a 0.10 mrad offset in the scattering angle (red squares) as 
a function of e at all 5 incident energies. The black squares are the nominal cross 
sections relative to themselves. 

In estimating the scale, random, and slope uncertainties in the sensitivity of the 
cross sections to a 0.10 mrad and 0.18 mrad angle offset, the same procedure described 
in section 13.4.21 is used. The sensitivity of the cross sections to a 0.10 mrad offset in 
the scattering angle has been studied. Figures I4~7I and l4~%l show the relative difference 
between the nominal cross sections, determined at the nominal scattering angle and 
energy, and the cross sections with a 0.10 mrad change in the scattering angle for the 
right and left arms at all kinematics. A 0.10 mrad angle fluctuation changes the cross 
section by (0.10-0.12)% for the right arm and (0.02-0.10)% for the left arm. Such 
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Figure 4.8: The relative difference between the left arm nominal cross sections and 
cross sections with a 0.10 mrad offset in the scattering angle (red squares) as a function 
of e. 

change in the cross sections was applied as a random uncertainty to each e point. 

Similarly, the sensitivity of the cross sections to a 0.18 mrad offset in the scattering 
angle has been studied. Figures T4.9I and 14. 101 show the relative difference between the 
nominal cross sections and the cross sections with a 0.18 mrad angle offset for the 
right and left arms and at all kinematics. The right arm results show an average 
scale, random, and slope uncertainties of 0.20%, 0.02%, and 0.67%, respectively. 
While, the left arm results show an average scale, random, and slope uncertainties of 
0.13%, 0.01%, and 0.18%, respectively. 



4.5 VDC Tracking Efficiency 

The VDC tracking inefficiency is defined as the fraction of good events where we do 
not reconstruct a track in the VDCs. It is determined by taking the fraction of the 
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Figure 4.9: The relative difference between the right arm nominal cross sections and 
cross sections with a 0.18 mrad offset in the scattering angle (red squares) as a function 
of e at all 5 incident energies. The solid blue line is a linear fit to the data. 
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Figure 4.10: The relative difference between the left arm nominal cross sections and 
cross sections with a 0.18 mrad offset in the scattering angle (red squares) as a function 
of e. The solid blue line is a linear fit to the data. 
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zero- and multiple-track events that passed through both planes of the two VDCs 
and caused a trigger in the scintillators planes: (N zero + N mu [ tip i e ) / 'N to tai where N zero 
and N mu iti P i e are the number of events with zero- and multiple-track, respectively, 
and Ntotai is the total number of good events. That in turn defines the VDC tracking 
efficiency €vdc as the fraction of one track events N one /N to tai- in the analysis of E01- 
001, we require only one-track events. Therefore, we reject multi-track events and 
correct for lost events. A good-track event is defined as a one-track event that fired 
both VDCs planes and was in the the fiducial area, i.e, hit the central paddles of the 
Si and 62 scintillator planes. That requires a hit in paddle 3 or 4 of each scintillator 
plane. In addition, tracks with zero or many hits (hits per plane >30) are excluded. 
In addition on both arms, particle identification (PID) cuts using A\ and A2 aerogels 
were applied to exclude ir + . See section I4~TI for a more detailed description of the PID 
cuts. Figure E. Ill shows the distribution of the number of reconstructed tracks that 
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Figure 4.11: Track multiplicity in the left arm spectrometer for run 1597, kinematics 
i. 
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satisfy the good-track event condition defined above for the left arm spectrometer. 
Although events with one track are desired, events with zero- or multiple-track are 
also present. Table POl shows the number of zero-, one-, multiple-track events, and 
the VDC tracking efficiency for the right and left arm spectrometers as determined 
for a single run. The efficiency is calculated and applied to each run separately. The 
zero- and multiple-track inefficiencies for the left arm are < 0.21% and < 0.38%, 
respectively, and < 0.017% and < 1.26% for the right arm, respectively. 

In order to determine the uncertainty in the VDC tracking efficiency, the average 
zero- and multiple-track inefficiencies for all the runs in all the kinematics were de- 
termined by applying several cuts. These cuts include: cutl = fiducial cut, cut2 = 
fiducial plus no zero-hit cut, cut3 = fiducial plus no multi-hit cut, and cut4 = fiducial 
plus no zero-hit plus no multi-hit cut. Figures 14.121 14.131 14.141 and 14.151 show the 
average zero- and multiple-track inefficiencies as determined under these cuts for the 
right and left arms, respectively. For the zero-track inefficiency, Figures r4.12l and l4.14l 
show that cutl and cut3 produce similar results as do cut2 and cut4. Clearly, exclu- 
sion of the no multi-hit cut does not alter the results and suggests the importance of 
the fiducial and no zero-hit cut applied only. Therefore, in the analysis of the E01- 
001, cut4 was used to determine the zero- and multiple-track inefficiencies. Clearly 
the zero-track inefficiency for the left arm is larger than that for the right arm and 
has a range of (0.05-0.20)% with Q 2 dependence but not any e dependence. There 
is an estimated random or point-to-point uncertainty (fluctuation) of < 0.01% in the 
average zero-track inefficiency for both arms. Therefore, the random uncertainty will 
be set to 0.0%. The average zero-track inefficiency on the right arm does not have 
any significant e dependence or scale offset. Therefore, a scale and slope uncertainties 
of 0.0% are assigned. On the other hand, a scale and slope uncertainties of 0.1% and 
0.0% are assigned for the left arm. 

As for the multiple-track inefficiency, Figures 14.131 and 14.151 show that the differ- 
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ent cuts seem to produce relatively close results. The multiple-track inefficiency for 
the right arm has a range of (0.60-1.30)% mainly rate dependence with an average 
inefficiency of ~ 1.0% and shows an e dependence of 0.05%. If a 10% measurement is 
assumed, this will yield a 0.10% scale uncertainty. A slope uncertainty ofO.05%/0.07 
~ 0.70% is assigned when we consider the Ae range of 0.07. On the other hand, the 
multiple-track inefficiency for left arm has a range of (0.15-0.30)% which varies with 
Q 2 . A scale uncertainty of ~ 0.10% will be assigned as well. The multiple-track inef- 
ficiency shows insignificant e dependence and a slope uncertainty of 0.0% is assigned. 
In order to determine the random uncertainty in the multiple-track inefficiency for 
both arms, the average multiple-track inefficiency from all the cuts (cutl-cut4) were 
plotted relative to cut4 and in order of increasing e. This was done for the average 
zero-track inefficiency as well. Figure 14.161 and 14.171 show the results. Again, the 
multiple-track inefficiency for the left arm is larger than that for the right arm and 
has some Q 2 dependence as well. The average multiple-track inefficiency for the right 
arm has on the average a random fluctuation of 0.01% (set to 0.0%), while the left 
arm average multiple-track inefficiency shows a random fluctuation of 0.02%. 



Number of Tracks 


Right Arm 


Left Arm 


Zero 


10 


63 


One 


198563 


154028 


Two 


2231 


240 


Three 


16 


1 


Four 








Five 








Total 


200820 


154332 


Zero-track fraction (%) 


0.00498 


0.0408 


Multiple-track fraction (%) 


1.119 


0.156 


tVDC 


0.988 


0.998 



Table 4.3: The number of zero-, one-, multiple-track events, and the VDC tracking 
efficiency for the right and left arm spectrometers for run number 1597, kinematics i. 
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Figure 4.12: The average zero-track inefficiency for the right arm with several cuts 
applied plotted as a function of kinematics. The points are sorted according to the 
kinematics of the left arm. Kinematics 1-5 correspond to Q 2 = 2.64 GeV 2 , while 
kinematics 6-9 (10-12) correspond to 3.20 (4.10) GeV 2 . For each Q 2 value, the points 
are sorted by e (low to high). See Table I3~T1 for detail. Note that cut4 was used in 
the analysis. 
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Figure 4.13: The average multiple-track inefficiency for the right arm with several 
cuts applied plotted as a function of kinematics. The points are sorted according to 
the kinematics of the left arm. Kinematics 1-5 correspond to Q 2 = 2.64 GeV 2 , while 
kinematics 6-9 (10-12) correspond to 3.20 (4.10) GeV 2 . For each Q 2 value, the points 
are sorted by e (low to high). See Table I3~T1 for detail. Note that cut4 was used in 
the analysis. 
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Figure 4.14: The average zero-track inefficiency for the left arm with several cuts 
applied plotted as a function of kinematics. Kinematics 1-5 correspond to Q 2 = 2.64 
GeV 2 , while kinematics 6-9 (10-12) correspond to 3.20 (4.10) GeV 2 . For each Q 2 
value, the points are sorted by e (low to high). See Table l3~T1 for detail. Note that 
cut4 was used in the analysis. 
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Figure 4.15: The average multiple-track inefficiency for the left arm with several cuts 
applied plotted as a function of kinematics. Kinematics 1-5 correspond to Q 2 = 2.64 
GeV 2 , while kinematics 6-9 (10-12) correspond to 3.20 (4.10) GeV 2 . For each Q 2 
value, the points are sorted by e (low to high). See Table I3~T1 for detail. Note that 
cut4 was used in the analysis. 
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Figure 4.16: Top: The right arm average zero-track inefficiency with several cuts ap- 
plied (cutl = solid black circles, cut2 = solid red squares, cut3 = solid green triangles) 
relative to cut4 (open blue squares) plotted as a function of kinematics. Note that 
the open blue squares are the average zero-track inefficiency from cut4 relative to 
itself. The points are sorted according to the kinematics of the left arm. Kinematics 
1-5 correspond to Q 2 = 2.64 GeV 2 , while kinematics 6-9 (10-12) correspond to 3.20 
(4.10) GeV 2 . For each Q 2 value, the points are sorted by e (low to high). See Table 
13. II for detail. Bottom: The right arm average multiple-track inefficiency with several 
cuts applied (cutl = solid black circles, cut2 = solid red squares, cut3 = solid green 
triangles) relative to cut4 (open blue squares) plotted as a function of kinematics. 
The points are sorted according to the kinematics of the left arm as discussed in top 
caption above. 
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Figure 4.17: Top: The left arm average zero-track inefficiency with several cuts ap- 
plied (cutl = solid black circles, cut2 = solid red squares, cut3 = solid green triangles) 
relative to cut4 (open blue squares) plotted as a function of kinematics. Note that 
the open blue squares are the average zero-track inefficiency from cut4 relative to 
itself. Kinematics 1-5 correspond to Q 2 = 2.64 GeV 2 , while kinematics 6-9 (10-12) 
correspond to 3.20 (4.10) GeV 2 . For each Q 2 value, the points are sorted by e (low to 
high). See Table I3~D for detail. Bottom: The left arm average multiple-track ineffi- 
ciency with several cuts applied (cutl = solid black circles, cut2 = solid red squares, 
cut3 = solid green triangles) relative to cut4 (open blue squares) plotted as a function 
of kinematics. The points are sorted according to the kinematics of the left arm as 
discussed in top caption above. 
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4.5.1 VDC Hardware Cuts Inefficiency 

In the previous section we discussed the VDC tracking efficiency. Even after selecting 
single-track event, not every reconstructed track corresponds to the true particle 
trajectory . A typical single-track event making an angle of 45° with the VDC surface 
has a multiplicity, number of hit wires in the cluster that were fired by the event, of 
4-6 [751 r7Tj . Some of the good single-track events do not have enough hits in the 
VDC cluster or have noise hits not associated with the true particles. This can lead 
to tracks that do not reproduce the true particle trajectory. Such tracks lead to a 
long tails in the distribution of the reconstructed physical quantities. 

In order to eliminate these long tails, we apply a VDC hardware cuts. These cuts 
require that the number of clusters for each VDC plane = 1, minimum number of 
hits per cluster in each VDC plane = 3, maximum number of hits per cluster in each 
VDC plane = 6. However, removing these long tails comes at the expense of loosing 
good events which must be accounted for. We account for such loss by multiplying 
the beam charge Q with the efficiency of the VDC hardware cuts 6vdch- 

Figure 14*. 1 81 shows the inefficiency of the VDC hardware cuts applied for the right 
arm as a function of e for all 5 incidents energies. This inefficiency represents the 
fraction of good events that were lost after applying the VDC hardware cuts to 
remove the long tails seen in the distribution of the reconstructed physical quantities. 
In determining the inefficiency, we applied a set of physics/data quality cuts to get rid 
of junk events (see section EIH for a full list of cuts applied) and the results have some 
cut dependence. The results show an average inefficiency of 11.35% giving eyDCH = 
0.8865. We do not expect this efficiency to depend on e and the results are consistent 
with no slope. Therefore a slope uncertainty of 0.0% is assigned. The results also 
suggest a scale and random uncertainties of 0.5% and 0.1%, respectively, due to the 
dependence on the quality cuts applied. 

Figure HUH] shows the inefficiency of the VDCs hardware cuts applied as a function 
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Figure 4.18: The right arm VDCs hardware cuts inefficiency as a function of e at all 
5 incident energies (solid red circles). The solid green line is the average inefficiency 
of 11.35% used in the analysis. 
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Figure 4.19: The left arm VDCs hardware cuts inefficiency as a function of e for all 
kinematics. The Q 2 = 2.64 GeV 2 kinematics are shown as solid black circles, Q 2 = 
3.20 GeV 2 as solid red circles, and Q 2 = 4.10 GeV 2 as solid blue circles. For each Q 2 
value, the points are sorted by e (low to high). The solid green line is the inefficiency 
we use in the analysis as given by 7.0%+0.25%(£-0.50). See text for details. 
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of £ for all kinematics for the left arm. Each of the three Q 2 values gives a small slope. 
If we average the three slopes we get an average slope of (0.60 ± 0.14)%. However, 
most of the slope comes from the lowest e point. If we remove the lowest point for 
each Q 2 , the result is consistent with zero slope. We do not know of any reason why 
the efficiency of the VDCs hardware cuts should depend on e. However, we do know 
of reasons why our estimate of the efficiency might have an e dependence and that 
is because of the overlap of the VDCs cuts with some of our physics/data quality 
cuts. So we take our best set of cuts, and apply a correction that has a large enough 
uncertainty to be consistent with no e dependence. Therefore, we apply a small slope 
with a large uncertainty of (0.25 ± 0.25)% to try and cover both possibilities. The 
green line shown in Figure 14.191 is the inefficiency or correction we apply to each 
kinematics and is given by: 7.0% + 0.25%(e-0.5). The results also suggest a scale 
and random uncertainties of 0.5% and 0.1%, respectively. 

4.6 Scintillator Efficiency 

The scintillator planes and trigger system are described in detail in section l3.9.2l Too 
little or no energy deposited by the charged particles in the scintillator paddles, the 
inefficient transmission of light emitted by the charged particle inside the scintillators 
paddles to the PMTs, or inefficiencies of the PMTs used, will result in scintillator 
inefficiency and hence trigger inefficiency. Trigger inefficiency is estimated by taking 
the fraction of good events that were not counted by the main physics trigger Ti(T 3 ) 
in the right (left) arm but still caused a trigger as a T2{T±) event type in the right (left) 
arm. 

In calculating the scintillator efficiency, single-track events that fell inside the 
scintillator boundaries were kept by projecting the track to the scintillator plane and 
excluding events that missed the detectors. This is done by applying cuts on the VDCs 
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planes multiplicity and focal plane x and y axes, xfp and yfp, and angles, xpfp and 
ypfp. PID cuts were used to separate 7r + background from protons. Photomultiplier 
hits in the left and right side of each scintillators plane were required and software 
solid angle cut of 1.6 msr was applied to assure that these protons within the desired 
spectrometer's angular acceptance. 

Each scintillator plane efficiency esi t s2 was calculated using: 

esiS2 = N ^ +N ^ (4 4) 

where iVj (i — 1, • - • , 5) is the number of events of trigger type i corrected for prescal- 
ing factor and electronic and computer deadtimes that fell inside the scintillators 
boundaries as defined by Table 14.41 It must be clear that event type T 2 (4) requires a 
hit in S\ or S2 but not both. In determining the efficiency of Si, N 2 (&) in equation 
(14.4)1 is the number of events that fired S 2 scintillator plane but not necessarily S\. 
Similarly, if the efficiency of S2 is desired, N 2 u) will be the number of events that 
fired Si scintillator plane but not necessarily S 2 - 



Arm 


Si plane boundary 


S2 plane boundary 


Right 


-1.05 < xfp + 1.381xpfp < 0.90 
-0.18 < yfp+1.381ypfp < 0.18 


-1.30 < xfp + 3.3Uxpfp < 1.00 
-0.32 < yfp + 3.3Uypfp < 0.32 


Left 


-1.05 < xfp + 1.287xpfp < 0.90 
-0.18 < yfp + 1.287ypfp < 0.18 


-1.30 < xfp + 3.Ulxpfp < 1.00 
-0.32 < yfp+3.Ulypfp < 0.32 



Table 4.4: Track projected from the focal plane to the location of S\ and 5*2 scintil- 
lators and required to be inside of the scintillator plane. The focal plane variables 
used are the x and y axis coordinates, xfp and yfp (measured in meter), and angles, 
xpfp and ypfp. 



Figure E.20I shows how the efficiency for the two scintillators planes and es 2 , 
as defined by equation ()4.4|) . for the right arm spectrometer is calculated. Similar 
plots (not shown) were generated for the left arm spectrometer and for each run. The 
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Figure 4.20: Events that fell within the scintillators boundaries as defined by Table 
I4.4I and needed for determination of Si and S% scintillators efficiencies for run 1597, 
kinematic i. Top left: Events of type (7\ + T 5 ) that fell within the right arm Si 
scintillator boundary as projected from the focal plane (xfp+1.381xpfp). Top right: 
Events of type (T x + T 2 + T 5 ) that fell within the right arm Si scintillator boundary as 
projected from the focal plane (xfp+1.381xpfp). The Si scintillator efficiency is ratio 
of events in the two plots. Bottom left: Events of type (7\ + T5) that fell within the 
right arm S2 scintillator boundary as projected from the focal plane (xfp+3.314xpfp). 
Bottom right: Events of type (Ti+T 2 +T 5 ) that fell within the right arm S2 scintillator 
boundary as projected from the focal plane (xfp+3.314xpfp). The S2 scintillator 
efficiency is ratio of events in the two plots. Similar plots can be generated for the 
left arm scintillators with T^) — > T 3 ^. 
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final and total scintillators efficiency for any run is the product of the two scintillators 
efficiencies or es 1 es 2 . Table 1431 lists the results for eg i; es 2 , and £s±£s 2 f° r both the 
right and left arm spectrometers for a selected run. The scintillators efficiencies were 
typically > 99.5% and > 99.6% for Si and S2, respectively, and were calculated and 
applied to each run separately. 



Arm 


£51 


tS2 




Right 
Left 


0.996 
0.999 


0.997 
0.999 


0.994 
0.998 



Table 4.5: The Si and S 2 efficiencies, eg 1 and es 2 , and total scintillators efficiency, 
C5 1 es 2 j f° r f ne right and left arm spectrometers as determined for run 1597, kinematics 
i. 

Figures 14.211 and 14.221 show the efficiency of the two scintillators planes, and 
€s 2 , for all the LH 2 and dummy data runs for both arms. For the both arms, we 
estimate the random and scale uncertainties to be 0.05% and 0.10%, respectively. 
There is not any significant e dependence observed in the results, therefore a 0.0% 
slope uncertainty is assigned. 

4.7 Particle Identification Efficiency 

Particle identification (PID) cuts are needed in order to obtain a clean proton sam- 
ple and the efficiency of these cuts must be accurately determined as well as any 
misidentification of other particles as protons. Two spectrometers were used, the low 
Q 2 spectrometer (right arm) and the high Q 2 spectrometer (left arm), and PID cuts 
and their efficiencies are needed for both. 
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Figure 4.21: The right arm Si and S2 scintillators efficiency for both LH 2 (red) and 
dummy (black) for all of the elastic kinematics. The efficiency is plotted as a function 
of run number. The run number is sorted according to the kinematics of the left arm 
and by increasing Q 2 . For each Q 2 , the kinematics are sorted by e (low to high). See 
Table EH] for details. 
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Figure 4.22: The left arm Si and S2 scintillators efficiency for both LH 2 (red) and 
dummy (black) for all of the elastic kinematics. The efficiency is plotted as a function 
of run number. The run number is sorted by increasing Q 2 and for each Q 2 by e (low 
to high). Kinematics a and b have lower Si scintillator efficiency for their LH 2 and 
dummy runs compared to the rest of the kinematics because one of the six paddles 
in the Si scintillator plane was tilted at the time of the run forming a gap with the 
adjacent paddle. Because the gap was far from the elastic acceptance, we increased 
the efficiency of the Si scintillator by 0.15% for all the runs in kinematics a and b. 
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4.7.1 The efficiency of the right arm (3 cut 

In order to determine the efficiency of the (3 cut, ep, that we are applying, a two- 
dimensional plot of the sum of the ADC signals for the A 2 aerogel, A 2 ADCSUM, 
and the energy deposited in the first right arm scintillator, S\ dE/dx, was generated. 
Figure 14.231 shows such a plot. In Figure 14.231 deuterons, protons, and pions have 
been identified and separated from each other by applying two-dimensional cuts on 
A 2 ADCSUM vs SidE/dx. A plot of the right arm (3 for the three particle types 
was then generated using two-dimensional cuts shown in Figure 14.231 Figure 14.241 




Figure 4.23: Two-dimensional plot of A 2 ADCSUM vs S\ dE/dx for run 1730 from 
kinematics n. Deuterons (red boundary), protons (blue boundary), and pions (7r + ) 
(green boundary) have been identified and separated from each other using the tight 
two-dimensional cuts shown. The cut for deuterons is offset to reduce proton con- 
tamination since the efficiency of this cut is not important. 
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shows the contribution of each particle to the full ft spectrum. Pions (which usually 
measure 0.85< ft <1.20) and deuterons (usually 0.02< ft <0.45) leak into the protons 
area with 0.45< ft <0.85. Although the tails of the distributions (tails of the pions, 
deuterons, and protons distributions) with ft < zero and ft > 1.0 are unphysical and 
are due to finite timing resolution of the PMTs, we account for the events in these 
tails when we calculate the pion and deuteron contaminations as well as the proton 
loss. 
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Figure 4.24: The contribution of pions (green), deuterons (red), and protons (blue) 
to the full right arm ft spectrum (solid black line). 

Based on Figure I4.241 the deuteron and pion contamination and proton loss were 
determined to be: 

• Total number of deuterons = 1 185 

• Number of deuterons into protons area = 9 
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• Number of deuterons into pions area = 

• Total number of pions = 3051 

• Number of pions into protons area = 132 

• Number of pions into deuterons 14 

• Total number of protons = 64587 

• Number of protons into pions area = 62 

• Number of protons into deuterons area = 271 

• Number of protons into protons 64089 

• Deuteron Contamination = (Number of deuterons into protons area/Number 
of protons into protons area)xl00% = (9/64089) x 100% = 0.014% 

• Pion Contamination = (Number of pions into protons area/Number of protons 
into protons area) x 100% = (132/64089) x 100% = 0.206% 

• Proton Efficiency (e^) = (Number of protons into protons area/Total number 
of protons) = (64089/64587) = 0.99243 

This analysis was done at all 5 incident energies and for multiple runs. Table 14.61 
summarizes the results. The deuteron contamination is clearly negligible while the 
pion contamination is small and is negligible after applying the aerogel A 2 ADCSUM 
cut. 

Figure 14.251 shows the deuteron contamination, pion contamination, proton loss, 
and the efficiency of the f3 cut, e@, for a multiple runs from all the kinematics. The 
results are plotted in order of increasing e. The efficiency of the (3 cut is typically 
> 99.2% with an average random fluctuation of 0.05% and without any noticeable e 
dependence. 
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Figure 4.25: The right arm deuteron contamination (top left), pion contamination 
(top right), proton loss (bottom left), and efficiency of the (3 cut, e@, (bottom right) for 
sample runs from all kinematics. The points are sorted according to the kinematics of 
the left arm, with three runs shown for each kinematics setting. Runs 1-15 correspond 
to Q 2 = 2.64 GeV 2 , while 16-27 (28-36) correspond to 3.20 (4.10) GeV 2 . For each Q 2 
value, the points are sorted by e (low to high). See Table I3~T1 for details. 
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Run 


Incident Energy 


Maximum Deuteron 


Maximum Pion 


Proton 


Number 


(MeV) 


Contamination 


Contamination 


Efficiency 






(%) 


(%) 


M 


1252 


2260.00 


0.0100 


0.1886 


0.9943 


1653 


2844.71 


0.0086 


0.1827 


0.9932 


1730 


4702.52 


0.0140 


0.2060 


0.9924 


1772 


1912.94 


0.0000 


0.2290 


0.9934 


1823 


3772.80 


0.0120 


0.1937 


0.9925 



Table 4.6: The incident energy, deuteron contamination, pion contamination, and 
proton efficiency for the (3 cut. 



Having determined the proton efficiency from the (3 cut, ep, we would like to 
determine the uncertainty in ep. Note that ep was determined assuming a pure proton 
sample. However, this is not the case since any pion and deuteron contamination will 
change the value of ep. In order to determine the size of possible pions/deuterons 
contamination, the j3 spectrum is used and the number of deuterons and pions which 
can be misidentified as protons determined. This was done for multiple runs at all 5 
incident energies. First, the maximum number of deuterons (maximum bin content) 
in the deuterons area is determined. Then the maximum number of protons which 
leaked into the deuterons area is determined. This is an upper limit estimate and 
based on the content of the bin directly below the bin that contained the maximum 
number of deuterons. This procedure is illustrated below for run number 1250 as in 
Figure 14.241 and for deuterons: 

• Total number of deuterons in all (3 spectrum = 1326. 

• Maximum deuterons bin content in deuterons area = 506. 



Maximum protons into deuterons area (upper limit estimate) = 65. 
Ratio of protons to deuterons in deuterons area = a = 12.85%. 



130 

• Number of deuterons misidentified as protons = N Deuterons = a (total number 
of deuterons) ~ 170. 

• Total number of protons in all (3 spectrum = N proton = 161019. 

The same procedure used to determine the uncertainty in is used to determine 
the uncertainty in €a 2 - Similarly the number of pions (in the pions area) which are 
misidentified as protons are determined. The ratio of protons to pions is estimated 
to be 7 = 0.69%, giving the number of pions misidentified as protons = N n + = 
7 (total number of pions) ~ 47. The scale uncertainty is determined as (iV£> etlteroris 
+ N n +) x 100%/ Np ro ton = 0.135%. The scale uncertainty in the efficiency of the (3 
cut as determined at all the incident energies was fairly constant and an overall scale 
uncertainty of 0.15% was assigned. 

4.7.2 The efficiency of the right arm Aerogel (A2 ADC- 
SUM) cut 

We looked at the A 2 ADCSUM for each of the five runs listed in Table 14.61 For 
each run, we overlayed the contribution from each of the three particles to the A 2 
ADCSUM spectrum using tight cuts on f3 and Si dE/dx as shown in Table I4~7I Figure 
14.261 shows the full A 2 ADCSUM spectrum and its constituents. In addition, the PID 
boundary cut used in the analysis is also shown. 



Particle 


(3 range 


Si dE/dx range 


Deuterons 
Protons 
Pions 


0.20< (3 <0.40 
0.60< (3 <0.70 
0.80< (3 <1.20 


2500< St dE/dx <3500 
1000< St dE/dx <1800 
400< St dE/dx <900 



Table 4.7: The tight (3 and S\ dE/dx cuts range. 



131 



10 



10 



10 



10 



10 



Full A 2 ADCSUM Spectrum 

Protons 

Pions (tc + ) 

DeUterons 



A 2 ADCSUM PID Cut = 1250 




1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

A 2 ADCSUM 



Figure 4.26: The full A 2 ADCSUM spectrum (blue). Also shown the deuterons 
(cyan), protons (red), and pions (tt + ) (green) contribution to the full A 2 ADCSUM 
signal for run 1730 from kinematics n. The magenta line at A 2 ADCSUM = 1250 is 
the PID cut used. 
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It can be seen that some of the protons leak into the pion area, A 2 ADCSUM > 
1250. Similarly, some pions make it into the proton region, A 2 ADCSUM < 1250. 
Only events with A 2 ADCSUM < 1250 are kept knowing that these good protons 
are definitely contaminated with some pions and deuterons while events with A 2 
ADCSUM > 1250 will be rejected with some good protons lost. 

The contamination of pions and the protons efficiency were determined for the 
data shown in Figure l4".26l as the following: 

• Total Protons = 66503 

• Protons greater than 1250 = 322 

• Protons less than 1250 = 66181 

• Total Pions = 3208 

• Pions greater than 1250 = 3176 

• Pions less than 1250 = 32 

• Pion Contamination = (Number of pions less than 1250/Number of protons less 
than 1250) x 100% = (32/66181) x 100% = 0.05% 

• Proton Efficiency (e^ 2 ) = (Number of protons less than 1250/Total number of 
protons) = (66181/66503) = 0.99516 

Similar analyses were done on runs taken at the other four incident energies. Table 
14.81 summarizes the results. The final proton efficiency for each of the 5 e points is 
defined as the product of the proton efficiency as determined from the (3 cut and the 
proton efficiency as determined from the A 2 ADCSUM cut or e^e^. As can be seen 
from Table 14.91 the final proton efficiency as determined for all of the 5 e points is 
fairly constant as expected since all data are taken at a fixed proton momentum. The 
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particle identification efficiency, epjn, for the right arm is taken as the product of the 
average value of ep and €a 2 or epir> =< £p >< £a 2 > — 0.9886. Table [Ql summarizes 
the results. It can be seen from Table l4~£l that the pion contamination is < 0.1% and 
can be reduced by applying the f3 cut. Such contamination can be reduced further 
by kinematics cuts and endcaps subtraction and it is taken to be negligible in the 
analysis. 



Run 


Incident Energy 


Pion 


Proton 


Number 


(MeV) 


Contamination 


Efficiency 






(%) 


( e A 2 ) 


1252 


2260.00 


0.0410 


0.9960 


1653 


2844.71 


0.0353 


0.9950 


1730 


4702.52 


0.0500 


0.9952 


1772 


1912.94 


0.0404 


0.9960 


1823 


3772.80 


0.0397 


0.9950 



Table 4.8: The pion contamination and proton efficiency as determined using the A 2 
ADCSUM cut. 

Figure 14.271 shows the deuteron contamination, pion contamination, proton loss, 
and the efficiency of the A2 ADCSUM cut, 6a 2 , for multiple runs from all the kinemat- 
ics. Again, the results are plotted in order of increasing e and the efficiency of the A2 
ADCSUM cut is typically > 99.4% with an average random fluctuation of 0.05% and 
without any significant e dependence. Similarly, and based on the same procedure 
used earlier with the f3 cut efficiency, the scale offset uncertainty was determined by 
looking at the A2 ADCSUM spectrum and determining the number of deuterons and 
pions which can be misidentified as protons. The scale uncertainty in the efficiency of 
the A2 ADCSUM cut as determined at all the incident energies was fairly close and 
an overall scale uncertainty of 0.25% was assigned. Therefore, a scale uncertainty of 
~ 0.30% (adding in quadrature the scale uncertainty from ep and ca 2 ) on the right 
arm particle identification efficiency epm is assigned. 
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Figure 4.27: The right arm deuteron contamination (top left), pion contamination 
(top right), proton loss (bottom left), and efficiency of the A 2 ADCSUM cut, e^ 2 , 
(bottom right) for sample runs from all kinematics. The points are sorted according 
to the kinematics of the left arm, with three runs shown for each kinematics setting. 
Runs 1-15 correspond to Q 2 = 2.64 GeV 2 , while 16-27 (28-36) correspond to 3.20 
(4.10) GeV 2 . For each Q 2 value, the points are sorted by e (low to high). See Table 
13. II for details. 
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Run 


Energy 


Proton Efficiency 


Proton Efficiency 


Final Proton 


Number 


(MeV) 






Efficiency ep €a 2 


1252 


2260.00 


0.9943 


0.9960 


0.9903 


1653 


2844.71 


0.9932 


0.9950 


0.9882 


1730 


4702.52 


0.9924 


0.9952 


0.9876 


1772 


1912.94 


0.9934 


0.9960 


0.9894 


1823 


3772.80 


0.9925 


0.9950 


0.9875 






< e p > = 0.9932 


< e A2 > = 0.9954 


e PID = 0.9886 



Table 4.9: The final proton efficiency as determined from the efficiency of the (3 and 
A 2 ADCSUM cuts. 

4.7.3 The efficiency of the left arm A x Aerogel (A x ADCSUM) 



cut 

Understanding the PID cuts and efficiencies for the left arm is not as straightforward 
as it is for the right arm. In the right arm, we had a good separation of particles 
using /3, dE/dx, and the A 2 ADCSUM cuts. For the left arm, using the A\ aerogel 
to separate the protons from pions works well, however, the f3 cut cannot be used 
because at such high Q 2 protons and pions have almost the same value for (3. 

The coincidence runs with protons in the left spectrometer and electrons in the 
right spectrometer were used to generate a pure proton spectrum. First, all the 
runs in each of the three coincidence kinematics were added. A plot of logio(Ai) for 
each coincidence kinematics was generated. For simplicity, I will refer to logio(^i) 
by A\ ADCSUM throughout this section. The total number of events under the 
Ai ADCSUM spectrum as well as the number below and above A\ ADCSUM = 
2.544, which is the PID cut used for A\ ADCSUM and is represented by the blue 
line in Figure I4.28| was determined for pure proton sample. A pure proton sample 
was generated by applying cuts on the elastic peak AP and/or y-coordinate of the 
extended target length y tg . The results were insensitive to the exact cuts applied. 
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This allowed for determination of the fraction of protons (proton inefficiency) that 
leaked into A\ ADCSUM > 2.544 which was fairly constant regardless of the various 
cuts applied. The two Q 2 = 2.64 GeV 2 spectra yield consistent proton inefficiencies, 
(1.037±0.029)% and (1.030±0.013)%, while the Q 2 = 4.1 GeV 2 has a higher proton 
inefficiency of (1.940±0.054)%. Figure ED shows the A 1 ADCSUM spectrum for 
kinematics coinl (see Table l3~Tj) . We do not have any coincidence kinematics with 
Q 2 = 3.2 GeV 2 to extract clean proton sample. Instead, we take the A\ ADCSUM 
spectrum for each of the kinematics at this Q 2 value, and get clean proton sample by 
applying kinematics cuts on the elastic peak AP, -15. 5< AP <30.5 MeV, and y tg , 
-0.0044< y tg <-0.001 m, to eliminate the endcaps contribution. The same procedure 
applied above was used and gave a value of (1.550±0.049)% for the proton inefficiency. 



I 10 5 
10 4 
10 3 
10 2 
10 



1 1.5 2 2.5 3 3.5 4 4.5 

log 10 (Aj ADCSUM) 

Figure 4.28: The Ai ADCSUM spectrum for a pure proton sample from coincidence 
kinematics coinl. The blue line is the PID cut used. 

We also need to determine the 7r + contamination. To do that, we need the number 
of 7r + and fraction of tt + that are identified as protons (ir + efficiency). We compared 
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Figure 4.29: The Ai ADCSUM spectra from singles kinematics n (top black), coin- 
cidence kinematics coinl (top red), scaled kinematics n (top blue). The difference 
between the A\ ADCSUM spectrum from the scaled kinematics n and the coinl 
spectrum (blue minus red) is the pions spectrum (bottom magenta). 
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the A\ ADCSUM spectrum from trigger type T 3 events (LH 2 singles events, which 
include both protons and 7r + ) with that of the trigger type T 5 events (LH 2 coincidence 
events, which have only protons). The difference in the spectra should represent 
the A\ ADCSUM spectrum for the background 7r + . When we compare the two A\ 
ADCSUM spectra, we scale the A\ ADCSUM spectrum for the singles events so that 
the two peaks match. Figure I4".29l illustrates this procedure where the A\ ADCSUM 
spectrum with trigger type T 3 events is compared with that of trigger type T 5 events. 
The peak of the A\ ADCSUM spectrum from singles was scaled to match that of 
coincidence. The difference between the A\ ADCSUM spectrum from the scaled 
singles and the coincidence spectrum should represent the A\ ADCSUM spectrum 
for pions. The total number of pions under the A\ ADCSUM spectrum and the 
number of pions below and above A\ ADCSUM = 2.544 are determined. The pion 
efficiency is defined as (number of pions below A\ ADCSUM = 2.544)/(total number 
of pions). 

Having determined the pion efficiency and proton inefficiency, the pion contamina- 
tion is then calculated for each kinematics based on the A\ ADCSUM spectrum and 
without applying any kinematics cuts or subtracting off the yield from the aluminum 
endcaps using the following procedure: 

• Total number of protons ~ number of events < 2.544. 

• Total number of pions ~ number of events > 2.544 - (total number of protons) x 
(proton inefficiency). 

• Number of pions < 2.544 = (pion efficiency) x (total number of pions). 

• Pion contamination = (number of pions < 2.544)/(total number of protons). 

For the proton inefficiency we considered the average weight of (1.031±0.012)% 
for Q 2 = 2.64 GeV 2 , while for Q 2 = 4.10 GeV 2 , a value of 1.88% was used. For Q 2 = 
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3.20 GeV 2 , the average value of the proton inefficiency for Q 2 = 2.64 and 4.10 GeV 2 
or 1.45% was finally used. Figure I4*.30l shows the results for the proton inefficiency as 
determined by applying cuts on the elastic peak AP and y tg for multiple runs selected 
at the three Q 2 values. For Q 2 = 2.64 GeV 2 , the proton inefficiency seems to have, on 
the average, a random fluctuation and scale offset uncertainties of < 0.1% each, and 
no significant e dependence. For Q 2 = 3.20 GeV 2 , the proton inefficiency shows, on 
the average, a random fluctuation and scale offset uncertainties of < 0.1% and (0.10- 
0.20)% respectively, and does not show any significant e dependence based on the data 
centered around the 1.5% proton inefficiency. Finally, for Q 2 = 4.10 GeV 2 , the proton 
inefficiency shows a random fluctuation and scale offset uncertainties of < 0.1% and 
0.10% respectively, and does not show any significant e dependence. Therefore, and 
mainly based on the results of the Q 2 = 3.20 GeV 2 point, a scale, random, and slope 
uncertainties of 0.20%, 0.10%, and 0.0%, respectively, will be assigned for the A 1 P1D 
efficiency at all three Q 2 points. 

To estimate the pion contamination, we also need the pion efficiency. We take 
the worst case estimate to give a conservative limit on the pion contamination. The 
pion efficiency for Q 2 = 2.64 and 4.10 GeV 2 was found to be (4.230±1.32)% or 5.5% 
as a worst case (defined as the maximum value) and (8.935±1.93)% or 11% as a 
worst case, respectively. For Q 2 = 3.20 GeV 2 , pion efficiency value of 11% was used 
as a worst case. Table 14.101 summarizes the results for the pion efficiency, proton 
inefficiency, pion contamination at the three Q 2 values. The proton efficiency ea x , 
= 1.0 — proton inefficiency, or the efficiency of the A\ ADCSUM cut is also listed. The 
pion contamination is always < 1.0%. Applying kinematics cuts and subtracting the 
target endcaps reduces the pions contaminations as most of the pions are generated 
by bremsstrahlung scattering in the endcaps. Thus, pion contamination in the final 
results, after the dummy subtraction, is always < 0.1%. 



140 





o 

Q 


= 2. 


64 Ge\ 


o 

r 








Q 2 = 3. 


20 Ge\ 


r 
















Q 2 -4. 


10GeV 2 


— 








































— 




















































— 


































i 


> 




< 












— 




















































— 
































< 


» 


















— 




































< 


» 














— 




















































































• 




































> 








i 




p 










































« 












































































































< 


































4 


i 












»--- 


< 


i 






























1 


i 




< 


i 


< 












































i 


► 




i 


► 



















































































































































1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 

Run Number 



Figure 4.30: The left arm proton inefficiency for sample runs. Runs 1-6 correspond 
to the two coincidence kinematics with Q 2 = 2.64 GeV 2 , 7-18 correspond to singles 
kinematics with Q 2 = 3.20 GeV 2 , and 19-21 correspond to coincidence kinematics 
with Q 2 = 4.10 GeV 2 . For each Q 2 value, the points are sorted by e (low to high). 
See Table I3~T1 for details. Note that for Q 2 = 3.20 GeV 2 , we do not have a direct 
measurements, and so the scatter of the results shown in the plot is due to uncertainty 
in the procedure used rather than an actual variation in the proton inefficiency. 
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Table 4.10: The pion efficiency as determined using the A\ ADCSUM cut. In addition, 
the proton inefficiency, pion contamination, and proton efficiency are also listed. 



4.8 Proton Absorption 

Once the protons are struck, they have to travel through material in the target, spec- 
trometer, and the detector stack before they can cause a trigger in the scintillators. 
While traveling, some protons undergo nuclear interaction in these materials and 
hence get absorbed or scattered. Such loss due to the nuclear interaction is called 
nuclear absorption. In order to determine the proton absorption for each arm, we 
have to account for all the target materials the protons have to pass through on their 
way to the scintillators. Also the thickness (Ml ESI EH EH1 ESI an d density of each 
absorber jHTl EH EI3 EI] , and the mean free path between nuclear collision A (effective 
absorption length) are needed in order to calculate the absorption. 

To determine A, the mean free path between nuclear collision (total interaction 
length), Ay, and the mean free path between inelastic interactions (inelastic interac- 
tion length), A/, must be known [H71 IHH] . In the analysis of the E01-001 experiment, A 



142 

was determined using two different definitions. The first definition, A is estimated as 
2\ T \ I /(\ T + A/) assuming that half of the elastic and inelastic scattering contribute 
to the absorption. The second definition, A is estimated as the average of two lengths 
or (Ar + A/)/2 since we cannot determine exactly the full contribution from the elastic 
scattering to the absorption. I will refer to A estimated by the first (second) definition 
as X act (X avg ) where act and avg stand for actual and average. 

The ratio of X/X is calculated for each absorber first using X act and then X avg . 
Here X is the product of the absorber's density and thickness. The ratios are then 
added, i.e, ^™ =1 (^«/Ai) where i runs over all absorbers n. The proton absorption is 
given by: 

proton absorption = 1.0 — e^^^ 1 ^/^ , (4.5) 

where e~^=^ Xi / x ^ defines the proton transmission. The final proton absorption used 
is the average value of the calculated proton absorption from the two definitions or 
5.19% for the right arm and 4.91% for the left arm giving a proton absorption cor- 
rection of C Absorption = 0.948 and 0.951 for the right arm and left arm, respectively. 
Tables 14.111 and 14.121 list the absorbers used and their properties for both arms. Un- 
certainty in the final proton absorption used rise due to our lack of precise knowledge 
of the actual thickness of each absorber used and the use of different cross sections in 
determining the mean free path lengths discussed above. Therefore, a scale offset of 
1% is assigned. Due to the fact that the scattered proton will have to travel different 
path length inside the LH 2 target at each kinematics, calculations of the proton ab- 
sorption inside the LH 2 target at different scattering angles for both arms show an e 
dependence uncertainty of 0.03% and 0.10% for the left and right arm, respectively. 
Random uncertainty of 0.0% is assigned. 
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Table 4.11: The absorbers used and their properties for the right arm spectrometer. 
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Table 4.12: The absorbers used and their properties for the left arm spectrometer. 
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4.9 Target Length Correction 

Figure E.31I shows the geometry of the 4 cm LH2 cell used in the experiment. The 
cell wall is made of 0.14 mm Al. The cell has two endcaps. The upstream endcap 
(not shown) is the beam entrance window and it is made of Al 7075 T6 with 0.142 
mm thickness, while the downstream endcap window is uniformly machined Al in the 
shape of a hemisphere and has a thickness of 0.15 mm and radius R = 20.33 mm. The 
length of the central axis of the cell is 40.18 mm (black dashed line) |8*9*1 1901 I^T| . 

The 4 cm dummy target used is made of Al 6061 T6 with density 2.85 g/cm 3 . The 
thickness is 0.2052 g/cm 2 for the upstream foil and 0.2062 g/cm 2 for the downstream 
one. The distance between the two foils is (40±0.13) mm [SHI HUB EU EH] • 



DownStream Endcap 




R = 20.33 mm 



Beam 



Beam 



(40.18-R) mm 



40.18 mm 



R = Radius 

X = Beam Offset 



Figure 4.31: Geometry of the 4 cm LH2 cell. The electron beam (blue arrow) is going 
from the right to the left side with an offset X (green line) from the central axis of 
the cell (black dashed line). 
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If the beam (the blue arrow going from left to right) is displaced a distance X 
(green line) from the central axis, then the length that the beam will travel in the 
target is L = Li + (40.18 - R) mm where Li < R and is shown as a magenta arrow. 
This length is effectively the new target length. A survey showed that the target was 
displaced down (2.0±0.5) mm relative to the center of the beam. Furthermore, the 
beam was rastered to produce an ~ 2 mm x 2 mm spot. Therefore, the average beam 
path length in the target (new target length) is then L = <Lx> + (40.18 - R) where 
<Lx> is given by: 

[* ma * V-R 2 - X 2 dX 

< Lt > = Xmm y , (4.6) 

[* max dX 

and X min and X maa . are the minimum and maximum beam displacement. 

The target length correction needed is defined as SL = (40.18-L)xl00%/40.18. 
Table l4~T3*l lists the results of the calculations. The nominal value of SL, SL = 0.266%, 
gives a target length correction to Q e ff of Ctl = 0.9973, which changes by ~ 0.12% 
(average difference from the nominal value) for a 0.50 mm beam shift. Therefore, a 
scale uncertainty in the target length of 0.12% is assigned for both arms. 

The beam position measurements are discussed in detail in section 13.51 Figure 
13.151 shows the x and y coordinates of the beam at the target. The beam was well 
focused on the target with an average beam drift of 0.30 mm. This drift in the beam 
position translates into an uncertainty in the effective target length and introduces 
a correction to Q e //- Table |4~T41 lists the results. The target length correction SL 
changes by ~ 0.074% (average difference from the nominal value) for a 0.30 mm beam 
drift. Therefore, a random uncertainty in the target length of 0.074% is assigned for 
both arms. 
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0.27651 



Table 4.13: The target length correction SL as determined for several beam offset 
ranges. 
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1.3 


3.3 
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40.0411 


0.345694 


Average 








0.26987 



Table 4.14: The target length correction 5L as determined for several beam drift 
ranges. 
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4.10 Computer and Electronics Deadtime 

The computer deadtime is the portion of the time that the DAQ system is unable to 
record an event because it is busy recording another event. Scalers are used to count 
the total number of triggers generated for each event or trigger type. By knowing the 
number of events recored by the DAQ system, N 1 DA q, and the total number of events 
fed to the DAQ, N^ otal , as recorded by scalers for a particular event type {% where 
% = 1, • • • , 5), one can account for the missing events and determine the computer 
livetime CLT = pSiN D AQ/N Tota i. The computer deadtime is then defined as t\ C dt = 
1.0 - CLT where psi is the prescale factor for that event type. 

The electronics deadtime is the portion of the time during which the triggers are 
missed because the hardware is busy when a second event comes in. When an event 
causes a trigger, the logic gate in the trigger is activated and the TDC signal stays 
high for a fixed time r, which is the gate time width of the logic signal and has a 
value of (100±20) nsec. If a second event tries to trigger and activate the gate within 
that time window r, the second event will be ignored. If the event rate is R, then 
the probability of having n counts in a time t is given by the Poisson distribution 
P{n) = (Rt) n e~ Rt /n\, and the probability distribution for the time between the events 
is P(t) = Re~ Rt . For small electronics deadtime, the fraction of measured events or 
the electronics livetime ELT is equal to the probability that the time between events 
will be greater than r or: 



and the electronics deadtime is then defined as t/edt = 1-0 - ELT. 

The raw event rates were < lKHz for the left arm and < 30KHz for the right 
arm, respectively, but the right arm was heavily prescaled. The computer deadtime 
is in the range of (9—20)% for the left arm and (2—16)% for the right arm with an 




(4.7) 
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average of ~ 10% correction on both arms. The computer deadtime is well measured 
and known to within 1.0%. Therefore a 1.0% uncertainty in a 10% correction yields a 
scale uncertainty of 0.10% for both arms. Similarly, we assign a slope uncertainty of 
0.10% for the left arm, and ~ 1.0% (0.10%/0.07) for the right arm when we consider 
the Ae range of 0.07. A random uncertainty in the computer deadtime of 0.0% is 
assigned for both arms. 

The electronic deadtime is in the range of (0.001—0.008)% for the left arm with 
a random uncertainty, determined based on the 20% uncertainty in r, of 0.002% (set 
to 0.0%), and (0.03—0.30)% for the right arm with a maximum random uncertainty 
of 0.06%. If we take the average random uncertainty of 0.04% for the right arm and 
consider the Ae range of 0.07, we estimate the slope uncertainty to be 0.04%/0.07 ~ 
0.50%. The slope uncertainty is 0.0% for the left arm. A 0.0% scale uncertainty in 
the electronic deadtime is assigned for both arms. 

4.11 Beam Current Monitor Correction 

The beam current measurement was described in detail in section 13.61 The basic 
technique of calibrating the beam current monitors is well documented jB3J E21 EH]- 
The idea is to normalize the BCMs cavities to the OL02 cavity (Unser monitor) at 
the injector. During a calibration run, the beam was interrupted using a Faraday cup 
inserted after the OL02 cavity so that the zeros of the cavities and the Unser could 
be determined. The beam current was first set to 80 yuA for one minute and then the 
Faraday cup was inserted for another minute. The beam current was then stepped 
down to 70 /j,A and the Faraday cup was inserted again for another minute. This 
procedure was done down to a beam current of 10 /zA and then reversed up to beam 
current of 80 /iA. 

While the beam was stepped down through various currents, the EPICS informa- 
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tion such as the values of the upstream and downstream BCM voltage, Unser current, 
injector Faraday cup current, and the OL02 current were saved to a file. The EPICS 
BCM constants were determined from these values as: 

constant = — — ° L ° 2 — ^ (4.8) 

< V cavity > - zero offset 

where < Ioloi > an d < V cavity > are the average value for the OL02 current and 
cavity voltage, respectively. Table l4~T5l summarizes the results. The zero offsets for 
the cavities and Unser were determined from the beam off periods when the Faraday 
cup was inserted and the calibration constants and the zero offsets for converting 
the voltage-to-frequency (V-to-F) scalers to charge were also determined to give an 
average current of: 

S ,^ aler - zero offset 

Average current = e . (4.9) 

constant 

The values for V-to-F zero offsets and constants are listed in Table 14.161 Based on 
the analysis of all of the calibration runs, the gain and offset drifts were stable to 
within 0.1% and 0.01%, respectively. The effect of the offset drift is negligible, so a 
random uncertainty in the charge measurement of 0.1% is assigned for both arms. It 
is reported that jHlj the BCMs with the Unser are calibrated at the same time as 
the charge monitors every 2-3 months, and the results are stable to within ±0.5%. 
Therefore, a scale uncertainty in the charge measurement of 0.5% is assigned for both 
arms. 



Constant Type 


Constant Value 


EPICS Upstream cavity 
EPICS Downstream cavity 


77.44 ± 0.04 
79.08 ± 0.04 



Table 4.15: The EPICS BCM calibration constants. 



151 



Constant Type 


Constant Value 


Constant Zero Offsets 


V-to-F Ulx 


1331.2 ± 0.7 


152.78 


V-to-F U3x 


4094.8 ± 2.1 


163.88 


V-to-F UlOx 


12451.7 ± 9.6 


360.37 


V-to-F Dlx 


1353.4 ± 0.7 


34.63 


V-to-F D3x 


4190.0 ± 2.2 


110.01 


V-to-F DlOx 


13187.7 ± 10.2 


307.91 



Table 4.16: The V-to-F calibrations constants and offsets. 



4.12 Target Boiling Correction 

The LH 2 target density po can decrease due to boiling caused by the electron beam 
current. In order to study the boiling effect in the LH 2 , runs were taken at several 
beam current values using LH 2 and carbon targets. We define the normalized yield Y 
as the total number of events N detected in a chosen acceptance region by applying 
kinematics and acceptance cuts and then normalized to the effective beam charge 
Qeff or: 

N 

Y = — . (4.10) 

We// 

The normalized yield as a function of current showed a linearly decreasing re- 
lationship. The deviation of the normalized yield from 1.0 over the current range 
represents the instability of the target density due to target boiling effect. The target 
density p was parameterized as a function of current / as: 

p(I)=p (l.0-BI) , (4.11) 



where po = p(I = 0) and B is the slope giving a target boiling correction Ctb = 
p(I) j po as listed in Table 14.171 for the different kinematics. 
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Kinematics 




Ctb ± SCtb 


6CTB(relative) % 


a and run< 1269 


30 


0.9925 ± 0.0045 


± 0.60 


a and run> 1269 


50 


0.9875 ± 0.0075 


± 0.30 


b 


50 


0.9875 ± 0.0075 


± 0.30 


i-r 


70 


0.9825 ± 0.0105 


± 0.00 



Table 4.17: The target boiling correction Ctb an d its uncertainty 5Ctb as determined 
for different kinematics. In determining 5Ctb — ISB, a SB value of 1.5%/100/xA is 
used. See text for details. Note that 5CTB(relative) is the uncertainty relative to 
data at 70yuA. 

In the boiling studies, we analyze the data in two different methods. First, we 
plot the events/charge from scalers corrected for current-independent rate (rate of 
cosmic ray triggers) as a function of current. Figure 14.321 shows the results. For 
carbon, the slope is consistent with zero (0.32 ± 0.32)%/100/xA, as it should be. 
For LH 2 , the slope is (1.38 ± 0.15)%/100/xA. However, there is a slight nonlinearity 
at lower current, possibly due to nonlinearity in the BCM calibration at very low 
current or uncertainty in the correction for cosmic trigger ray, so this may be a slight 
overestimate. This may also be a slight underestimate, since the thickness of the 
endcaps does not depend on current. 

In the second method, we plot the yield as defined by equation ()4.10|) as a function 
of current where N is the number of the T 3 type tracked events. For carbon, the check 
is not of good quality due to poor statistics. For LH 2 , the results yield a slope of 
(1.12 ± 0.65)%/100/xA, but the statistics are low and so the uncertainty is large. In 
addition, the results are fairly cut-dependent, yielding slopes up to 3%, but with large 
uncertainties. 

Clearly, neither test is ideal and both give lower corrections than observed by previ- 
ous experiments. The scalers analysis is sensitive to the endcaps contribution and the 
current-independent rate. The tracked analysis is cleaner, but has low statistics. So 
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we take a correction of 2.5%/100/zA, but apply a large uncertainty (2.5±1.5)%/100/xA 
or (1.75±1.05) %/70/zA. Therefore a scale uncertainty of 1.05% in the target boiling 
correction is assigned. One thing to mention is that this effect cancels out if we take 
the ratio of the left to the right arm cross sections, and even for the left arm alone, 
it is mostly a normalization uncertainty since most of the kinematics were run at the 
same current and at a fixed Q 2 the proton corrections do not vary much with e. 
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Figure 4.32: Target boiling studies. The left arm normalized yield for LH 2 target 
(top) and carbon target (bottom) as a function of current. 
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We check this correction by looking at kinematics a and b, where part of the data 
is taken at 30/xA and part at 50/iA. The left arm has a relatively low rate (making it 
more sensitive to cosmic ray rates) and has a larger contribution from the aluminum 
endcaps, and so we focus on the right arm to check the correction. Also, since we 
mainly want this correction for the left arm, it is good to have an independent check 
of the boiling. For the right arm data, we see a clear difference between the yield at 
30 and 50/iA if we do not apply the boiling correction as can be seen in Figure d. 331 
The yield differs by roughly (0.6 ± 0.1)% for a 20/iA difference, yielding a boiling 
correction of (3.0 ± 0.5)%/100/iA. This is consistent with the (2.5±1.5)%/100/iA, so 
we stick with the more conservative number and use the consistency at 30 and 50/iA 
as a test of the correction. Figure 14.341 shows the yield after applying the boiling 
correction. 

Since most of the data are taken at 70/xA, we want the uncertainty relative to the 
data at 70//A. So we set the random uncertainty for all kinematics to 0.0% except 
that of kinematics a and b. We estimate the random uncertainty for kinematics a and 
b to be 0.45% or 1.5%(70/iA - 40/iA)/100/iA assuming an average current of 40/iA 
and 0.30% or 1.5% (20/iA)/ 100/iA, respectively. Furthermore, the beam was rastered 
to produce an ~ 2 mm x 2 mm spot with uncertainty in the raster size of < ±0.2 
mm which translates into a ~ 0.1% random uncertainty in the target boiling. 

4.13 Summary of Systematic Uncertainties 

Tables 14.181 and 14.191 summarize the systematic uncertainties for the right and left 
arm, respectively. The systematic uncertainties are discussed in the sections listed in 
these tables. The systematic uncertainty in each listed source is broken down into 
three types: scale, random, and slope. The contribution of each uncertainty type 
from all sources is then added in quadrature to form the total uncertainty in an for 
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Figure 4.33: The right arm LH2 yield before applying the target boiling correction 
for all the runs from kinematics a (1250-1285) and b (1293-1316) normalized to the 
yield of run 1250. The yield is shown in red(green) for runs taken at I — 30(50)/iA. 
Runs 1272-1285 from kinematics a are taken at I — 50/iA and their yield is shown 
in green as well. The solid black line is the average normalized yield for that current 
setting. The dashed black lines are the average uncertainty in the normalized yield. 
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Figure 4.34: The right arm LH2 yield after applying the target boiling correction for 
all the runs from kinematics a (1250-1285) and b (1293-1316) normalized to the yield 
of run 1250. The yield is shown in red(green) for runs taken at I = 30(50)yuA. Runs 
1272-1285 from kinematics a are taken at / = 50/iA and their yield is shown in green 
as well. The solid black line is the average normalized yield for that current setting. 
The dashed black lines are the average uncertainty in the normalized yield. 
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that uncertainty type. Only the range on total 5R a ndom is given in these tables, while 
the actual value of total 5 Random in &r at each e point is listed in Tables 15.51 and 15.61 
It must be mentioned that the scale uncertainty in the pion subtraction for the left 
arm is 0.20% for all kinematics except those of Q 2 = 4.10 GeV 2 which is 0.40%. This 
results in a Sgiope — 0.539% for all kinematics with Q 2 = 2.64 and 3.20 GeV 2 , and 
Ssiope = 0.641% for all kinematics with Q 2 = 4.10 GeV 2 . In addition, 5sio P e for Q 2 
= 4.10 GeV 2 is then scaled up by 25% to account for the Ae range difference among 
the three Q 2 points which results in a Ssiope = 0.801%. See Table l4~T9l for details. 
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Table 4.18: Summay of the systematic uncertainties for the right arm. Numbers in 
bold show the biggest contribution to that uncertainty type. Sources marked with 
(*) contribute to the uncertainty in the luminosity monitor. For (fa), the value is 
0% for all kinematics except a and b which is 0.45% and 0.30%, respectively. For 
(|6), the value is applied on kinematics-by- kinematics basis as shown in Figure I4~7I 
and it ranges from 0.097%-0.125%. For (fc), the value is applied on kinematics-by- 
kinematics basis as shown in Figure 13.111 and it ranges from 0.015%-0.029%. Since 
As ~ 0.07 for the right arm, the slope uncertainty is much larger for the right arm. 
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Table 4.19: Summay of the systematic uncertainties for the left arm. Numbers in 
bold show the biggest contribution to that uncertainty type. Sources marked with 
(*) contribute to the uncertainty in the luminosity monitor. For (fa), the value is 0% 
for all kinematics except a and b which is 0.45% and 0.30%, respectively. For (ffo), 
the value is 0.15% for all kinematics except those of Q 2 = 4.10 GeV 2 which is 0.30%. 
The random uncertainty in the 0.10 mrad angle fluctuations and 0.02% beam energy 
fluctuations is applied on kinematics-by- kinematics basis as shown in Figure 14.81 and 
Figure EH21 respectively. 



Chapter 5 Data Analysis II: Reduced 
Cross Sections Extraction 



5.1 Overview 



Correct Dummy Data 
For Effective Charge 
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Monte Carlo 
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Figure 5.1: E01-001 experiment analysis flow chart. 

In this chapter I will present the procedure used to extract the e-p reduced cross 
section or. Figure 15.11 shows a flow chart of the procedure used. First, ESPACE 
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generates HBOOKs for a list of LH 2 and dummy target runs in a particular kinematics 
which are read by Recon ( section 15.2)) . Recon applies a set of cuts and corrections 
on an event-by-event basis and generates a new HBOOK file for each run on the list. 
In addition, Recon generates the AP spectrum for each LH2 and dummy run on the 
list. For each kinematics the AP spectra from all of the LH 2 and dummy runs are 
added separately after applying a new set of cuts to form the final AP spectrum for 
both LH 2 and dummy fsection l5.3p . The final AP spectrum for both the LH 2 and the 
dummy is then normalized to its total Q e ff ( section l5.4j) . The final corrected dummy 
AP spectrum is subtracted from the final corrected LH 2 AP spectrum. Monte Carlo 
simulations of the 7£> — > 7i°p and 7p — > jp backgrounds are performed ( section l5.5.2|) . 
and then normalized and subtracted from the final corrected LH 2 AP spectrum to 
produce the residual elastic e-p spectrum. The net number of the elastic e-p events 
is then compared to the number of elastic events in the e-p peak as simulated using 
the Monte Carlo simulation program SIMC (section 15.5. 1|) in a window in the AP 
spectrum. The reduced cross section or is taken as the value of the input e-p cross 
section used in the simulation normalized by the ratio of number of elastic events in 
the data to that in the simulation (section l5.7)l . Finally, a linear fit of <jr to e gives 
rG 2 Mp (Q 2 ) as the intercept and G 2 Ep (Q 2 ) as the slope. 

5.2 Event Reconstruction by Recon 

The HBOOK file for each physics run generated by ESPACE (see section 14. 2j) was 
read using the in-house event analyzer Recon. Recon reads in the HBOOK file for a 
given run as Column- Wise Ntuples, applies a set of cuts on an event-by-event basis 
on all the physical variables saved in these Ntuples, and calculates the necessary 
physics quantities. A new HBOOK file for each run is then generated for the final 
data analysis. Table 15.21 lists the cuts applied on these physical variables. Events 
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that satisfy these cuts are accepted for the final data analysis. 

First, a single-track event type T!(T 3 ) for the right(left) arm is selected. The 
multiplicity of a cluster of hit wires in each plane is defined as the number of hit 
wires in the cluster. A single-track event making an angle of 45° with the VDC 
surface has typically a multiplicity of 5. We require the multiplicity of each of the 
four planes for both VDCs to be greater than two. The Q% cut is applied to reject 
events that scrape off the exit pipe of the Q 3 quadrupole. The Q 3 cut is defined as: 



where xfp, xpfp, yfp, and ypfp are the focal plane positions (xfp, yfp) and angles 
(xpfp, ypfp) variables. The Q% was applied to the left arm events only since there is 
no evidence that events were scraping off the exit pipe of the Q% quadrupole in the 
right arm spectrometer. It must be mentioned that the elastic events are far from the 
Qs edges and are not affected by the Q3 cut. Figure IS~2l shows how the Q3 cut was 
applied. 

The hourglass cut (all good trajectories that make it through the HRS form an 
hourglass pattern ~ 69 cm before the first VDC) is then applied to remove events 
which come from multiple scattering inside the spectrometer and located outside the 
acceptance. The hourglass cut is defined as: 



where C\ has a value of 0.045(0.017) for the right(left) arm. Figure I5~31 shows how 
the hourglass cut is applied. 

Finally cuts on the target reconstructed variables such as the momentum accep- 
tance 6, the y-coordinate of the extended target length y tg , the out-of-plane angle 
6 tg , and the in-plane angle <p tg are applied. Figures 15.41 and 15.51 show the target 



Q 3 = a/ (xfp - 2.6Axpfp) 2 + (yfp - 2.64ypfp) 2 < 0.29 m , 



(5.1) 



yfp — O.Q9ypfp + 0.005 — C\ xfp — 0.69xpfp , 



(5.2) 
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Figure 5.2: The Q3 cut, solid red circle with radius of 0.29 m, applied to the left arm 
events (black) for a single run from kinematics k to reject events that were scraping 
off the exit pipe of the Q3 quadrupole. Only events within the area of the solid red 
circle are accepted. 
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Figure 5.3: The hourglass cut applied to the right and left arms events from run 1282. 
The right arm events before the cut (top left) and after the cut (top right). The left 
arm events before the cut (bottom left) and after the cut (bottom right). 
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reconstructed variables after the cuts have been applied. 
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Figure 5.4: The right arm target reconstructed variables with cuts from run 1282. 
The right arm momentum acceptance 5 (top left), in-plane angle 4> tg (top right), out- 
of-plane angle 6 tg (bottom left), and y-coordinate of the extended target length y tg 
(bottom right). 

It should be mentioned that the momentum acceptance, 5, showed an out-of-plane 
angle dependence on both the right and left arms, and a raster, ry, dependence on 
the left arm only. The momentum acceptance 5 had first to be adjusted to correct 
for the 6 tg dependence for the left arm as: 

6 1 = S + 0.025# t9 + 0.250# t 2 , (5.3) 
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Figure 5.5: The left arm target reconstructed variables with cuts from run 1282. 
The left arm momentum acceptance 5 (top left), in-plane angle 4> tg (top right), out- 
of-plane angle 9 tg (bottom left), and y-coordinate of the extended target length y tg 
(bottom right). 



while for the right arm: 



S 1 = S + O.O27420 ts 
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(5.4) 



where 8q and 8\ are the momentum acceptance before and after the adjustment, 
respectively. The raster correction was then applied as: 



5 2 = Si + constant x ry 



(5.5) 



where ry is the raster's ADC signal (raster number 2 or rast2adc): 

i (rast2adc- 1560) 



ry = 0.002- 



raster width 



(5.6) 



constant is the correction applied, and the raster width is the width of the ADC 
signal whose average width is 1560 and is energy dependent. Again, 8% and 82 are the 
momentum acceptance before and after the raster correction, respectively. Figure l5~B1 
shows the AP spectrum before and after the raster correction. Table 15.11 lists the 
values for the constant and raster width. 



Tiefenback Energy 
(MeV) 


constant 


raster width 


1912.94 


0.42 


-1400 


2260.00 


0.38 


+ 1100 


2844.71 


0.00 


-1000 


3772.80 


0.00 


+ 1200 


4702.52 


0.30 


-600 



Table 5.1: The Tiefenback energy, constant (correction applied) and raster width 
needed to correct the momentum acceptance 8 for the raster dependence. 



In order to reconstruct AP, equation ()3.2|) . the following procedure is used. The 
final momentum of the scattered protons P me asured (measured using the high resolution 
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Figure 5.6: Two dimensional histogram of AP (x-axis) vs raster ry (y-axis). Top 
black: the AP spectrum before the correction shows a raster ry dependence as the 
events centered around AP = MeV are slightly tilted counter clockwise. Bottom 
red: the AP spectrum after the correction where the events are tilted clockwise (see 
equation ([5.5)1 ) and become mainly centered around AP = MeV. 
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spectrometer) is reconstructed as: 

Pmeasured P. L(-R) (1 + S 2 ), (5.7) 

where Pl(r) is the central momentum setting for the left (right) arm spectrometer as 
listed by Table Kill The scattering angle of the proton is then reconstructed: 

C0S e p = c ^m±M^ , (5.8) 

where the plus(minus) sign is used with the right(left) arm respectively, and &r(l) is 
the right (left) arm spectrometer scattering angle as listed by Table I4~2l Finally, the 
momentum of the scattered protons P ca icuiated{6 P ) , equation (jH.ljl . is then obtained 



Cut Type 
Applied 


Left Arm 
Cut Value 


Right Arm 
Cut Value 


Event Type 
Number of Tracks 
VDC Multiplicity 
Qs (m) 
Hourglass (m) 
5 Momentum (%) 


3 
1 

Multiplicity>2 
<0.29 
<0.02sin# L +0.01 
-5.0<<5<5.0 


1 
1 

Multiplicity>2 

<O.O2sin0iH-O.O2 
-5.0<5<5.0 


PID Aerogel ADCSUM 
PID (3 

Number of Clusters/VDC Plane 
Minimum Number of Hits/Cluster 
Maximum Number of Hits/Cluster 

Vtg ( m ) 

Out-of-Plane Angle (9 tg ) (mrad) 
In-Plane Angle (<f>tg) (mrad) 


A!ADCSUM<350 

1 

3 
6 

-0.05<?/ tg <0.05 
-40.0<^ tg <40.0 
-lO.O<0t 9 <lO.O 


A 2 ADCSUM<1250 
0.45</3<0.85 
1 
3 
6 

-0.05<^<0.05 
-40.0<^ a <40.0 
-1O.O<0 49 <1O.O 



Table 5.2: Cuts applied in the analysis of the E01-001 data. Cuts in the upper half of 
table are applied in Recon and cuts in the bottom half are applied later as discussed 
in section 15.31 See text for more details. 
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using the measured scattering angle of the proton as given by equation (|5.8|) . 

5.3 Adding Histograms 

For the data analysis, the HBOOK file generated by Recon for each physics run in 
a given kinematics is read by a kumac file (code) written using PAW the Physics 
Analysis Workstation software. The PAW kumac starts by reading the run number 
and then loading its HBOOK file. A new set of cuts such as the PID cuts A2 ADCSUM 
and (3 cuts for the right arm and A x ADCSUM cut for the left arm (to select good 
protons) and tracking cuts (number of clusters for each VDC's plane = 1, minimum 
number of hits per cluster in each VDC's plane = 3, and maximum number of hits 
per cluster in each VDC's plane = 6) are applied to the LH2 (dummy) AP ntuple of 
that run. See section 14.71 for the range of PID cuts used. The AP ntuple for LH2 
or dummy is then projected onto a histogram where a distribution of the number of 
counts per one MeV bin can be plotted. This procedure was done for all the LH2 and 
dummy runs in that kinematics. The LH 2 and dummy AP histograms from all the 
runs are added to form the final Q e ff uncorrected LH 2 and dummy AP histogram. 

The LH 2 AP spectrum is made of several contributions. The main contribution 
is from elastic ep — *> ep scattering which has a peak and a radiative tail. In addition 
there are backgrounds due to quasi-elastic and inelastic scattering in the aluminum 
target windows and walls and protons generated the photoreactions ^p — > ir p and 

To account for the quasi-elastic contribution to the LH2 AP spectrum, events 
scattered from the aluminum target windows were measured using an empty dummy 
target. The windows of the dummy target were thicker than the windows of the LH 2 
target and the determination of the effective thickness of the dummy target, referred 
to as "dummy thickness" will be discussed is section 15.71 Figures 15.71 and 15.81 show a 
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Figure 5.7: Q e ff uncorrected right arm AP spectrum for LH 2 (red) and dummy (blue) 
targets from kinematics b. 




Figure 5.8: Q e ff uncorrected left arm AP spectrum for LH 2 (red) and dummy (blue) 
targets from kinematics b. 
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final Qeff uncorrected AP spectra for the LH 2 and dummy targets for both the right 
and left arm spectrometers. The resolution of the elastic peak is dominated by the 
angular resolution. The right arm elastic peak has a broader width than the left arm 
elastic peak. This is due to the fact that at low momentum, AP is more sensitive to 
the angular resolution of the spectrometer, causing a broadening of the AP peak. 

5.4 The Effective Charge Normalization 

After adding all the LH 2 and all of the dummy histograms from the runs in a given 
kinematics, the final Q e ff uncorrected AP histogram must be divided by Q e ff- The 
effective charge as defined in section |4~T1 was calculated for each LH 2 run using equa- 
tion (tO . 

The dummy target is an empty target with thicker aluminum windows than the 
actual entrance and exit windows of the LH 2 target. In determining Q e ff, equation 
(14.1)1 . we must correct for the effective thickness of the dummy target. The reason 
we use an effective thickness for the dummy target is the bremsstrahlung initiating 
reactions, primarily in the downstream endcap. The bremsstrahlung are different 
for LH 2 and dummy because of the different thicknesses and also the presence of 
the hydrogen. The dummy target spectrum is used to determine the shape of the 
endcaps contribution and the procedure of determining the effective dummy thickness 
is discussed in section l5~71 The final Q e ff is the sum of the effective charge from all 
the LH 2 (dummy) runs, i.e., 

Ei Qlff where Qi // is the effective charge for run 
number i and N run is the total number of LH 2 (dummy) runs. 

The default and initial value of 4.11 for the effective dummy thickness was used 
to calculate Q e ff for the dummy data. The procedure of extracting the final effective 
dummy thickness used to extract the reduced cross section <jr for each kinematics 
is discussed in section 15.7.21 This initial value seems to do a good job producing 
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the high AP side of the LH 2 spectrum which is mainly dominated by the quasi- 
elastic contribution from the endcaps. Figures 15.91 and I5.1UI show normalized LH 2 
and dummy AP spectra for the right and left arms, respectively. 

! 1 ° 7 
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10 5 
10 4 
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10 2 

-60 -50 -40 -30 -20 -10 10 20 30 40 50 60 
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Figure 5.9: Q e ff corrected right arm AP spectrum for LH 2 (red) and dummy (blue) 
targets from kinematics b. 

5.5 Monte- Carlo Simulations 

5.5.1 Elastic e-p Simulations (SIMC) 

The Monte Carlo simulation program SIMC is used to simulate elastic scattering for 
all kinematics and for both arms. The elastic e-p simulations are a crucial compo- 
nent of the analysis as they are used to extract the reduced cross section. SIMC was 
adapted from the (e,e'p) simulation program named SIMULATE that was written 
for SLAC experiment NE18 joPl Id"!] , and was converted to simulate the Jefferson 
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Figure 5.10: Q e ff corrected left arm AP spectrum for LH 2 (red) and dummy (blue) 
targets from kinematics b. 

Lab Hall C and A spectrometers. The two main components of SIMC are the event 
generator, which includes the cross-section weighting and radiative corrections, and 
the spectrometer models. First, SIMC randomly generates the energy and position 
of the incident electron at the target to match the energy and spatial spread of the 
actual beam. In doing so, SIMC takes into account the target length and the beam 
raster. The beam energy is then corrected for ionization losses in the target. More- 
over, SIMC randomly generates the momenta and angles of the scattered electron 
and proton vectors with a flat distribution over limits exceeding the actual experi- 
mental spectrometer acceptance. In our case, we generate electrons and protons in 
coincidence and use the Bosted cross section to weight each generated event. Having 
generated a basic event at the scattering vertex, SIMC allows for any or all of these 
events to emit real or virtual photons where the corresponding event vectors are ad- 
justed and radiatively corrected The scattered electron and proton vectors are 



174 

transported through the target where ionization energy losses and multiple scattering 
in the target material, cells, and chamber are applied. Finally, only the scattered 
protons are transported through the spectrometers. 

Transporting the protons through the spectrometer was done using the spectrom- 
eter optics models built in the Monte Carlo simulation program COSY [95J. COSY 
generates both the forward and backward sets of matrix elements to simulate the 
optical resolution of the magnetic systems in the spectrometer. Note that all the 
material the protons have to travel through in the spectrometer and the detector 
stack are included, and the multiple scattering due to all the material is calculated. 
However, we do not account for proton absorption in the simulations since we correct 
the data for that. The forward matrix elements transport the particle vectors from 
the entrance window of the spectrometer to its focal plane going through every major 
aperture in the spectrometer. SIMC assures that these particles have gone through 
each of these apertures by checking the acceptance of each aperture using a set of 
apertures cuts. These apertures include the front, middle, and back aperture of each 
magnetic element (QQDQ configuration), the aperture of the vacuum can after the 
Q3 quadrupole, the aperture of the rectangular collimator, and the aperture of each 
detector that was used as a fiducial cut. In addition, VDCs smearing is applied to 
the particle positions at the two VDCs to match those of the actual VDCs. The 
backward matrix elements then reconstruct or transport back the particle vectors to 
the target vertex where they are corrected for energy loss and multiple scattering and 
the elastic AP spectrum for a given kinematics is then reconstructed. 

When comparing the shape of the elastic peak from data to that of simulations, 
it was obvious that the width of the elastic peak from data was broader than that 
of simulations. In addition, there was a clear mismatch between the shape of elastic 
tail from data and that of simulations. To resolve these issues, we compare the AP 
spectrum from coincidence data to that of simulations since the coincidence data 
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represent a pure protons sample without any background. Since the peak width is 
limited by the angular resolution of the spectrometer, we smear the angular resolution 
of the spectrometer in the simulations using a symmetric gaussian function to better 
fit the observed peak widths. We also apply an additional non-gaussian smearing to 
better match the tails. Figures 15. Ill and 15 . 1 21 show the normalized LH 2 AP spectrum 
along with the elastic e-p simulations for the right and left arms, respectively. 




-40 -30 -20 -10 10 20 30 40 

Right AP (MeV) 



Figure 5.11: Q e ff corrected right arm LH 2 AP spectrum (red) and elastic e-p simu- 
lation (blue) for kinematics b. 

5.5.2 7J9 — > tt°p and ^p — > ^p Simulations 

When the electron beam passes through the target (protons in our case), electrons 
lose energy by radiating real photons which is known by bremsstrahlung. These real 
photons which impinge upon the target have a maximum energy at the beam energy. 
To model the AP spectrum for these high energy 7r° protons, we first calculate the 
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Figure 5.12: Q e ff corrected left arm LH 2 AP spectrum (red) and elastic e-p simulation 
(blue) for kinematics b. 

bremsstrahlung cross section and reconstruct the E 1 spectrum for these photons jHH| , 
and then randomly generate photons according to the E 1 spectrum. The next step 
is to uniformly and randomly generate protons over the acceptance and this is done 
according to the event generation procedure discussed before. The generated event 
is then weighted by an s~ 7 cross section dependence, as predicted by the high energy 
approximation and based on the constituent counting rules or s~ n , forming the shape 
of the AP spectrum used for the ^p — > 7t°p contribution. Note that the generated 
proton is transported through the spectrometer using the spectrometer optics models 
mentioned before. Finally, the shape of the 7p — > ir°p AP spectrum from simulations 
is normalized to that of data as discussed in section 15.7.31 

Previous experiments [HE] indicated that the ratio of the jp —>■ 7p to the jp — ► n p 
cross sections (relative cross section) was 1-5%. We took the ratio from the original 
reference jHE] arid plotted it as a function of beam energy. Note that the beam energy 
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range quoted was very similar to that of the EO 1-001. There was a clear dependence 
on energy, with the ratio well fitted by (0.92E — 1.2)%. Therefore, we rescaled the 
7£> — > 7r°p AP spectrum by (0.92E — 1.2)% to generate the Compton AP background. 
Figures 15. 13l and l5. 14l show the simulated AP spectrum for the jp — > n°p and 7p — > ^p 
backgrounds for the left arm, respectively. 
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Figure 5.13: Left arm AP simulation for 7p — > ir°p background (7r°AP) for kinematics 
b. 

5.6 Radiative Corrections 

In this section I will start by describing the radiative corrections procedures for in- 
clusive p(e,e')p scattering, and then discuss details of implementation for E01-001. 
The Lorentz invariant transition scattering amplitude equation (jl.5j) and hence the 
elastic e-p scattering differential cross section equation (|1.10J) were derived to lowest 
order in a including only the amplitude due to the exchange of a single virtual pho- 
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Figure 5.14: Left arm AP simulation for 7p — > 'jp background (Compton AP) for 
kinematics b. 

ton between the incident electron and struck proton as shown in Figure 11.11 This 
is known as the Born approximation or single-photon exchange. However, in reality 
higher order processes in a and beyond the Born approximation also affect the cross 
section such as radiation of additional real and virtual photons. These processes 
clearly must be included. The incident electrons radiate in the presence of a nuclear 
field due to changes in their velocities caused by Coulomb interactions. We refer to 
the radiation resulting from such deceleration of the electron by bremsstrahlung or 
"Braking Radiation" . Corrections for these higher order processes are referred to as 
radiative corrections. Radiative corrections are classified into internal and external 
radiative corrections. 

• Internal Radiative Corrections: When the incoming and outgoing electrons 
interact with the Coulomb field of the nucleus involved primary in the scattering 
process, this results in emission and reabsorption of virtual photons and emission 
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Figure 5.15: Feynman diagrams for the elastic e-p scattering including both the first- 
order (Born) and higher-order QCD radiative corrections (a)-(h). 

of real soft photons. Figure lB*. 15l shows the Born term and higher-order Feynman 
diagrams used in the calculation of the internal radiative corrections. Note that 
these corrections are corrections to the principle scattering vertex itself. These 
internal radiative corrections are also classified into elastic radiative corrections 
such as vertex corrections (for both the electron and proton) Figure 15.15( a) 
and (g), vacuum polarization or loop diagrams Figure 15. 15( b). and two-photon- 
exchange (TPE) Figure 15.15( e) and (f). The inelastic radiative corrections 
involve the emission of real photons such as internal bremsstrahlung (for either 
electron or proton) Figure 13. 15( c). (d), and (h). Higher-order processes due to 
three or more photon exchange (Coulomb corrections) and emission of multiple 
real photons are also possible. 

Corrections that involve only lepton-photon or lepton-lepton vertices (electron 
vertex and electron bremsstrahlung) are exactly calculable in QED. Corrections 
that involve photon-proton vertices (proton vertex, TPE, proton bremsstrahlung, 
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and Coulomb corrections) are not exactly calculable in QED and must be either 
estimated or measured. 

• External Radiative Corrections: When the incoming and outgoing electrons 
interact with the Coulomb field of a nucleus other than the one involved in the 
primary scattering vertex, this also results in emission of real photons. Such 
radiative processes are caused by bremsstrahlung in material that the electron 
passes through before and after the scattering. 

As we see, the radiative processes are real physical processes which depend on the 
target material, geometry, and kinematics of the scattering. They affect the data by 
modifying the cross section of the process and the kinematics (momentum, energy, 
and angle) of the incident electron or final state of the electron and proton. Emis- 
sion of these real photons causes the detected particle's momentum to be different 
from the actual momentum at the scattering vertex and therefore distorts the exper- 
imental spectra. On the other hand, amplitudes involving the emission of additional 
virtual particles such as TPE processes affect only the magnitude of the measured 
cross section. The procedure for doing such radiative corrections was first derived by 
Schwinger (33] and later modified by Mo and Tsai 38J. 

The 5-function that multiplies equation (Jl.llj) is used to assure elastic scattering in 
the case of the electron, that is, at a given energy E and angle 8 e , the elastic differential 
cross section is a 5-function in E' at E' = E e i astic = E — where E' is given by 
equation (|1.14j) . See section fOl for more details. Bremsstrahlung radiative processes 
cause the elastic cross section to change from the simple 5-function (modified by finite 
resolution, small angle scattering, and energy spread of the beam) to an asymmetric 
peak with an extended tail at lower energies as shown in Figures lo~7l and IB~51 in the case 
of the proton. Note that this radiative tail extends down to values of AP where other 
inelastic processes such as pion production yields high energy electron, or for the case 
of proton detection, high energy protons generated from photoreaction 7^ — > 7r°p and 
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IV ~ > IP also occur. The procedure of extracting the elastic e-p scattering reduced 
cross section or will be discussed in section 15.71 and sub-sections therein. 

To convert the measured cross section to the single-photon-exchange cross section, 
the Born cross section of 0(a), we use the radiative corrections procedure of Mo and 
Tsai [HE] as modified by Walker |13| 15*^] and Ent jMJ EH] and implemented in the 
elastic e-p simulation code SIMC jHOllS]- Therefore, we write: 



da 

dn 



da 



Measured (1 &corr) 



Born 



(5.9) 



where 5 corr = (6 int + S int + 5 ext ) and 5 int and 5 ext represent the internal and external 
radiative corrections, respectively. The additional term Si nt results from the improve- 
ments made to the internal radiative corrections by Walker. For corrections higher 
than 0(a 3 ), we exponentiate 5 corr and we write: 



da 

dn 



Measured 



n &corr 



da 

dn 



Born 



(5.10) 



5.6.1 Internal Radiative Corrections 5i nt 

The internal radiative corrections 6i nt are determined based on the modification made 
by Walker [TJ IMl and Ent jM] to the work of Mo and Tsai [Ml- They include 
the processes of: vacuum polarization, TPE, electron and proton vertex corrections, 
and internal bremsstrahlung. Note that the approximation used in equation (|5.1()j) 
for higher order corrections is mainly valid only for the infrared divergent terms 
and the error caused by neglecting the nondivergent terms is estimated to be < 1% 
|97| 19*51 ESI HHDU mi] ■ ^ ne contributions from vacuum polarization and electron vertex 
corrections are calculated exactly. However, the contributions from TPE are limited 
only to the infrared divergent contributions (nondivergent terms have been neglected). 
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Therefore, the contributions to 5i nt can be written as: 

$int = (~Z 7T^ n ( q) ~\~ ^ int.br ems , (5.11) 

7T V 9 / 

where 5i nt .bremss. is the internal bremsstrahlung contribution and will be discussed 
next. 

For the internal bremsstrahlung, the calculations are done under the assumption 
used by Tsai [32] that AE(1 + jj-) <C E' where M p is the mass of the proton, E is the 
incident energy of the electron, E' is the final energy of the electron, and AE is the E' 
cutoff of the elastic peak or AE = E e i ast — E cuto ff. The resulting expression is rather 
lengthy and complicated. In general it has a complicated dependence on the incident 
energy of the electron, final energy of the electron, final energy of the proton, and 
AE. To achieve uncertainty better than 1% in the radiative corrections, corrections 
to this approximation were made by Walker [TH] 151] as part of the improvements he 
made to the internal radiative corrections represented by 5i n t. 

5.6.2 Walker's Improved Internal Radiative Corrections 5i nt 

For a full description of the improvements made to the internal radiative corrections, 
the reader is referred to Walker f3J|53] where he discussed these corrections in detail. 
However, I will just summarize briefly these improvements. The term 5i n t is just 
an improvement on the precision to the original internal radiative correction. These 
improvements include correction for the AE(1 + ff^) "C E' approximation made 
by Tsai to the internal bremsstrahlung, inclusion of the qq and fi + fi~ contributions 
to the vacuum polarization or loop diagrams which were neglected previously, and 
correction for a sign error in Tsai's paper [32] for the Schwinger's correction [32] to 
the noninfrared divergent part of the soft photon emission cross section. 
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5.6.3 External Radiative Corrections 5 ext 

The external radiative corrections 5 ex t are determined based on the modification made 
by Walker ^3J|!3] and Ent [HI] to the work of Mo and Tsai [38J. The external radiative 
corrections are corrections applied to the cross section due to bremsstrahlung in the 
target material (major contribution), spectrometer material, and the small effects 
of the Landau tail of the ionization energy loss spectrum |lU2j . When the incident 
electron interacts with the Coulomb field of a nucleus other than the one involved 
primary scattering vertex, it will emit bremsstrahlung and hence lose energy. The 
cross section of an electron of initial energy E that scatters elastically from a proton 
and has a final energy E' = E e i ast — AE', when the electron is emitting bremsstrahlung 
photons with t{ and tf radiation lengths of material before and after the scattering, 
is: 

d 2 a(E,E') _ /RAE'\biti /RAE'\bfh 1 1 
dttdE' ~ \ E ) \ E elast ) r(l + biU) T(l + bft f ) 
rd*W hjlL , A* x ME-RAE> )Ml /gAgxi 
L dtt AE' \ E e i as t J dtt AE'^\ E )\ K J 

where the parameter b ~ 3/4 |lU3j . the parameter R represents the recoil of the 
proton and is ~ (-^) 2 , and the function <p ( R ^ E ' j = 0f~§) gives the shape of the 
bremsstrahlung spectrum which is normalized to one at ^ = zero where u is the lost 
energy of the electron after passing a thickness of t radiation lengths. The external 
radiative corrections 5 ex t can be then expressed in terms of the electron cross section 
equation (|5.12jl as: 



5 P ~t. = In 



dcj(E) 
dfl 



d?a{E,E') ■ 
if, -af> dYLdE' 



(5.13) 



where LT corr is the correction for the Landau tail as calculated from the Landau 
distribution |102j . Note that the integrand diverges as AE' — > zero due to the soft 
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photon bremsstrahlung spectrum. To handle this problem, the integral is usually 
calculated for AE' < 5 7 where 5 7 is the soft-photon cut. When calculating 5 ext , one 
needs an accurate geometrical model of the target along with precise knowledge of 
the materials of which the target is made. This information is built in the elastic e-p 
simulation code SIMC. 

5.6.4 E01-001 Specific Radiative Corrections 

As mentioned before, we use the radiative corrections procedure of Mo and Tsai as 
modified by Walker and Ent and implemented in the elastic e-p simulation code SIMC. 
Full coincidence (e,e'p) simulations are performed using the prescription of Ent [HI] 
taking into account bremsstrahlung from all three tails (bremsstrahlung from the in- 
cident electron, scattered electron, and scattered proton). The bremsstrahlung from 
these particles is calculated in the extended peaking approximation [HUED]- Further- 
more, only the contribution from the infrared divergent part of TPE is included. 

Based on equation (jl.34j) and because of the factor r that multiplies G 2 M , the cross 
section is dominated by G 2 Mp at large Q 2 . That is, ctr in the Born approximation has 
only a small e dependence at large Q 2 . Therefore, understanding the e dependence 
of the radiative correction S corr at large Q 2 becomes absolutely important as a few 
percent change in the e slope of or would have a sizable effect on — ■ 

The main e dependence to the radiative corrections comes from the internal 
and external bremsstrahlung corrections (and possibly the neglected TPE terms). 
Bremsstrahlung corrections are usually accounted for in the conventional radiative 
corrections procedure used and they enter differently depending on whether the elec- 
tron or proton are detected in the final state. On the other hand, several calculations 
in the 1950s and 1960s estimated the size of TPE to be extremely small (< 1%) 
|9T| EHl EM I100| 1101] and therefore was neglected in the conventional radiative cor- 
rections procedure used. An important point is that the infrared divergent part of 
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the TPE contribution cancels fully the corresponding infrared divergent contribution 
from the interference between the electron and proton bremsstrahlung. Furthermore, 
the finite part of TPE has always been neglected. Recently it has been found that 
the finite part of TPE has a significant e dependence and it should not be neglected 
|104U1051ll06| . Figure I3~B1 shows the radiative correction factor calculated for internal 
contributions only as a function of e for Q 2 = 2.64 GeV 2 . While the magnitude of the 
corrections is similar and both show an approximately linear dependence on e, the 
dependence is much smaller for protons, ~ —8%, than for the electron ~ 17%. Note 
that when detecting the electron, the e dependence of the bremsstrahlung correction 
exceeds the e dependence coming from Ge p which is not the case in the E01-001 
experiment. 

Previous Rosenbluth measurements [T3*l ITSI 14*7] estimated scale and random un- 
certainties in the radiative corrections of ~ 1.0% and ~ 0.60%, respectively. In these 
measurements, electrons rather than protons were detected. In our case, we assign 
the same scale uncertainty of 1.0% for both arms and take a 50% of the random un- 
certainty as quoted by the previous measurements or ~ 0.30%. This is due to the fact 
that detecting protons will yield a smaller e dependence bremsstrahlung correction as 
discussed above. We split the 0.30% uncertainty into a slope and random uncertain- 
ties with more emphasize on the slope contribution of 0.30% and random contribution 
of 0.20%. For the right arm, the same random uncertainty of 0.20% is assigned, while 
the slope uncertainty is expected to be small due to the small Ae range of 0.07. Note 
that for the left arm, the average Ae range of 0.70 is accounted for when the slope 
uncertainty is estimated. Therefore, we estimate the slope uncertainty for the right 
arm to be 3.0% (0.30% x 0.70/0.07) but use a smaller and more realistic figure of 
2.0% since the radiative corrections are expected to be small in a such limited Ae 
range. 

Finally and based on Figure 13.61 a relative ~ 30% radiative correction is needed 
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to bring the cross sections measured by detecting electrons and those measured by 
detecting protons into agreement. That is, before the radiative corrections are applied 
the electron and proton cross sections have an e dependence that differs by as much 
as ~ 30%. The extracted slope from the two measurements after the correction is 
applied agree to within uncertainty of 1-2%. If we assume an average uncertainty in 
the slope of 1.5%, this yields a fractional uncertainty in a 30% correction of 1.5%/30 
= 5%. For the proton, a 5% uncertainty in a 10% e dependence bremsstrahlung 
correction yields a ~ 0.50% uncertainty which is comparable to the combined 0.36% 
uncertainty estimated for the E01-001 measurements. It is important to mention that 
this is the first time a comparison between the electrons and protons measurements 
has been made. The excellent agreement in the extracted slope or G\ from the 
two measurements provides a strong test of the validity of the standard radiative 
corrections used. 

5.7 The e-p Reduced Cross Section <tr 
5.7.1 Overview 

In this section the procedure for extracting <7r is discussed. Figure 15.11 shows a flow 
chart of the procedure used. From Figures lB~7l and l5~Hl we see that the LH2 spectrum is 
dominated by the elastic e-p peak. In addition, there are backgrounds due to quasi- 
elastic scattering from the aluminum target windows and protons generated from 
photoreactions 7^ — > 7i°p and 7^ — > '-/p that contribute to the AP spectrum. There 
is no region where only one component, the elastic peak or inelastic backgrounds, is 
significant. Therefore, the spectra has to be unraveled and this is done with aid of 
the simulated spectra but first the simulated and measured spectra have to be put 
on the same energy scale as will be discussed below. 
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5.7.2 The Effective Dummy Thickness and Endcaps Subtrac- 
tion 

In this section the procedure of extracting the effective dummy thickness will be 
discussed. We start by subtracting the Q e ff corrected dummy spectrum using an 
initial dummy thickness factor ( endcap ra ^°) °^ ^-H from the Q e ff corrected LH 2 
spectrum. Figures 15.161 15.171 and 15.181 show this procedure for a low e left arm 
spectrum, a high e left arm spectrum, and a right arm spectrum, respectively. 

The resultant LH2 AP spectrum is then compared to that of elastic e-p simulation 
as generated by SIMC. The experimental e-p peak is shifted relative to that of the 
simulation. Therefore, we shift the position of the e-p peak from data to match that 
of simulation. Figures 15.191 15.201 and 15 . 2 ll illustrate this procedure. The contribution 
of the 7p — > ir°p and 7]? — > 7p backgrounds to the LH 2 AP spectrum is the difference 
between the dummy subtracted LH 2 AP spectrum and elastic e-p simulation. From 
Figure 15.211 we see that the contribution of these backgrounds to the right arm LH 2 
AP spectrum is almost negligible. The momentum acceptance of the spectrometers 
is ±5%. For the higher momentum protons detected in the left arm, the momentum 
acceptance is great enough to be flat over the region of interest. However, for the low 
momentum protons detected in the right arm, the falling acceptance at the lowest 
AP results in a distortion of the spectrum at the low momentum end. This region 
was not used in extracting the cross sections. However, because the shape of the 
right arm spectrum is almost the same at all the kinematics, this distortion would 
only introduce a scale uncertainty in the reduced cross sections. For the left arm, 
Figures 15.191 and I5.2UI show that the contribution of these backgrounds to the LH 2 
AP spectrum is significant. As the cross sections of the 7r° protons are forward 
peaked, the contribution to the LH 2 AP spectrum is larger for low e than for large 
e. 

After we shift the position of the e-p peak from data to match that of the simula- 
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Figure 5.16: The initial left arm dummy subtraction for kinematics b (low e). Top: 
the left arm LH2 AP spectrum (red) and dummy AP spectrum (blue). Bottom: the 
dummy subtracted LH2 AP spectrum. An initial dummy thickness factor of 4.11 was 
used. 



189 





Figure 5.17: The initial left arm dummy subtraction for kinematics m (high e). Top: 
the left arm LH 2 AP spectrum (red) and dummy AP spectrum (blue). Bottom: the 
dummy subtracted LH 2 AP spectrum. An initial dummy thickness factor of 4.11 was 
used. 
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Figure 5.18: The initial right arm dummy subtraction for kinematics b. Top: the 
right arm LH2 AP spectrum (red) and dummy AP spectrum (blue). Bottom: the 
dummy subtracted LH2 AP spectrum. An initial dummy thickness factor of 4.11 was 
used. 
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Figure 5.19: The left arm dummy subtracted LH2 AP spectrum and elastic e-p 
simulation for kinematics b (low e). Top: the left arm unshifted dummy subtracted 
LH2 AP spectrum (red) and elastic e-p simulation (blue). Bottom: the left arm 
shifted dummy subtracted LH 2 AP spectrum (red) and elastic e-p simulation (blue). 
This data are from kinematics b, using the nominal incident energy and scattering 
angle. The final version after applying the angular shift to the scattering angle (see 
section yields much smaller shift. 
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Figure 5.20: The left arm dummy subtracted LH 2 AP spectrum and elastic e-p 
simulation for kinematics m (high e). Top: the left arm unshifted dummy subtracted 
LH 2 AP spectrum (red) and elastic e-p simulation (blue). Bottom: the left arm 
shifted dummy subtracted LH 2 AP spectrum (red) and elastic e-p simulation (blue). 
This data are from kinematics m, using the nominal incident energy and scattering 
angle. The final version after applying the angular shift to the scattering angle (see 
section yields much smaller shift. 
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Figure 5.21: The right arm dummy subtracted LH 2 AP spectrum and elastic e-p 
simulation for kinematics b. Top: the right arm unshifted dummy subtracted LH 2 
AP spectrum (red) and elastic e-p simulation (blue). Bottom: the right arm shifted 
dummy subtracted LH 2 AP spectrum (red) and elastic e-p simulation (blue). This 
data are from kinematics b. 
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tion, we compare the two AP spectra over a region where the elastic e-p simulation 
AP bin content is larger than 10% of that of the peak. The number of LH 2 events, 
Nlh 2 i an d number of e-p elastic, N e _ p , for that window is used to determine the 
scaling factor 7 e „ p = (NLH 2 /N e -p) needed to scale the elastic e-p simulation to match 
that of LH 2 . Figures f5.22l and 15.231 illustrate this procedure. 

Knowing that the contribution from the elastic LH 2 events to the LH 2 AP spec- 
trum at large AP is small, the shape of the endcaps spectrum can be reliably deter- 
mined by normalizing the difference between the LH 2 AP and that of the scaled elas- 
tic e-p simulation "extracted dummy" to that of the original Q e ff corrected dummy 
spectrum in the range of AP>30 MeV. Figures 15.241 15.251 and 15.261 illustrate this 
procedure. 

The number of events in the extracted dummy spectrum, N extracted , and the num- 
ber of events in the original dummy spectrum, N origina i, for that window is used to 
determine the correction factor ^ d ummy = 'original /N extr acted) needed to correct the 
initial default value of the dummy thickness and hence define the new effective dummy 
thickness: 

. Default Dummy Thickness 4.11 
Effective Dummy Thickness = = , (5.14) 

^dumray ^duramy 

which is used to correct Q e ff for the dummy and hence scale the dummy target 
spectrum for the final subtraction. The original and extracted dummy spectra are 
shown for both arms in the bottom right hand side plots of Figures 15.241 15.251 and 
15.261 Figure 15. 2 71 shows the effective dummy thickness for all kinematics and for both 
arms. For the left arm, the effective dummy thicknesses are consistent with random 
fluctuation of (1-2)% from the average. The average effective dummy thickness is 
0.22% thicker than the default dummy thickness used initially. For the right arm, the 
effective dummy thicknesses are consistent with random fluctuation of (2-3)% from 
the average, but the average effective dummy thickness is 5.15% less than the default 
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Figure 5.22: Right arm AP spectrum for LH 2 data (black), elastic e-p simulation 
(red), and scaled elastic e-p simulation (blue) for kinematics b. 
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Figure 5.23: Left arm AP spectrum for LH 2 data (black), elastic e-p simulation (red), 
and scaled elastic e-p simulation (blue) for kinematics b. 
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Figure 5.24: The left arm effective dummy target thickness extraction for kinematics 
b. Top left: The left arm dummy subtracted LH2 AP spectrum (red) and scaled 
elastic e-p simulation (blue). Top right: the left arm dummy subtracted LH 2 AP 
spectrum minus scaled e-p AP spectrum. Bottom left: The difference between the 
left arm dummy subtracted LH 2 AP spectrum and scaled e-p AP spectrum (black) 
is then normalized to the left arm original dummy AP spectrum (red) in the range 
of 30-140 MeV. Bottom right: The left arm original dummy spectrum (red) and the 
new extracted dummy spectrum (blue). 
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Figure 5.25: The left arm effective dummy target thickness extraction for kinematics 
m. Top left: The left arm dummy subtracted LH 2 AP spectrum (red) and scaled 
elastic e-p simulation (blue). Top right: the left arm dummy subtracted LH 2 AP 
spectrum minus scaled e-p AP spectrum. Bottom left: The difference between the 
left arm dummy subtracted LH2 AP spectrum and scaled e-p AP spectrum (black) 
is then normalized to the left arm original dummy AP spectrum (red) in the range 
of 30-140 MeV. Bottom right: The left arm original dummy spectrum (red) and the 
new extracted dummy spectrum (blue). 
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Figure 5.26: The right arm effective dummy target thickness extraction for kinematics 
b. Top left: The right arm dummy subtracted LH 2 AP spectrum (red) and scaled 
elastic e-p simulation (blue). Top right: the right arm dummy subtracted LH 2 AP 
spectrum minus scaled e-p AP spectrum. Bottom left: The difference between the 
right arm dummy subtracted LH2 AP spectrum and scaled e-p AP spectrum (black) 
is then normalized to the right arm original dummy AP spectrum (red) in the range 
of 30-55 MeV. Bottom right: The right arm original dummy spectrum (red) and the 
new extracted dummy spectrum (blue). 
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dummy thickness value used initially. 

5.7.3 Subtracting the jp —> 7r°p and jp —> jp Backgrounds 

The next step is to subtract the protons generated from photoreactions 7p — > n p and 
IP IP from the LH2 AP spectrum. When the endcaps background and the scaled 
elastic e-p simulation events are subtracted from the LH 2 AP spectrum, the resultant 
spectrum represents the sum of the jp — > n°p and 7p — > ^p backgrounds. For the left 
arm, the sum of 7p — > 7r°p and ^p — > 7p spectra from simulations is normalized to the 
(LH2 - endcaps - scaled e-p simulation) spectrum in a window of -50< AP <-20 MeV. 
The normalized simulation is then subtracted. Figures T5.28l 15.251 and 15.301 show this 
procedure for three representative spectra. While the uncertainty in extracting the 
backgrounds is larger for large e, where the backgrounds contribution is small, even 
a large uncertainty in this small background yields a small uncertainty in the final 
extracted cross sections. The contribution of the 7^ — > n p and jp — > ^p backgrounds 
to the right arm LH2 spectrum is almost negligible although this contribution has a 
threshold that is close to the elastic peak as can be seen clearly from the bottom right 
hand plot of Figure 15.301 On the other hand, the contribution of the 7p — > n°p and 
IP ~~* IP backgrounds to the left arm LH 2 spectrum is much larger but with a lower 
threshold for the low e shown in the bottom right hand plots of Figures 15.281 

and 15.291 It is important to mention that for the right arm the same procedure for 
extracting the 7p — > 7c°p and 7p — ► jp backgrounds was used, and the background 
was found to be consistent with zero, with upper limit contribution of <<0.10%. 
Thus in the analysis, we do not make any correction for these processes and assign a 
0.05% random uncertainty. 

The procedure for unpealing the spectra can be iterated if necessary. However, 
the residual spectra after backgrounds and normalized simulated e-p are subtracted 
off show that this is not necessary. 
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Figure 5.27: The left and right arms effective dummy thickness. Top: The left arm 
dummy/endcap ratio plotted as a function of increasing e for the Q 2 = 2.64 GeV 2 
kinematics (solid red circles), Q 2 = 3.20 GeV 2 kinematics (solid blue circles), and Q 2 
= 4.10 GeV 2 kinematics (solid green circles). The solid black line is the unweighted 
average dummy/endcap ratio and the red dashed line is the initial default value used of 
4.11. Bottom: The right arm dummy/endcap ratio plotted as a function of increasing 
e for all the Q 2 = 0.5 GeV 2 kinematics. The same color was assigned on data so the 
kinematics on the right arm can be associated with their left arm counterpartner. 
Again, the solid black line is the average dummy/endcap ratio and the red dashed 
line is the initial default value used of 4.11. The uncertainty in the data is statistical 
only. 
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Figure 5.28: The left arm *yp — > n°p and 7p — * 7p backgrounds subtraction for 
kinematics 6 (low e). Top left: the LH 2 AP (black), extracted dummy (red), and 
scaled elastic e-p simulation (blue). Top right: the sum of the 7p — > 7r p and -yp — > -yp 
data backgrounds remaining after the subtraction of the extracted dummy and scaled 
elastic e-p simulation spectra from the LH2 AP (black), and the sum of the ^p — > 7r°p 
and 7p — > 7J9 from simulations (cyan). Bottom left: the sum of the 7p — > 7T°p 
and 7p — > 7p as in top right above (black) and simulations (cyan) from -50 to -20 
MeV. This window is used for normalization. Bottom right: the LH 2 AP (black), 
extracted dummy (red), scaled elastic e-p simulation (blue), and normalized sum of 
the 7p — > 7T°p and ^p — > 7J9 from simulations (cyan). 
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Figure 5.29: The left arm jp — ► 7i°p and yp — >■ 7£> backgrounds subtraction for 
kinematics m (high e). Top left: the LH 2 AP (black), extracted dummy (red), and 
scaled e-p elastic simulation (blue). Top right: the sum of the ■yp — > 7r°p and 7p — > 7p 
data backgrounds remaining after the subtraction of the extracted dummy and scaled 
e-p elastic simulation spectra from the LH2 AP (black), and the sum of the jp — ► 7r°p 
and 7p — > 7p from simulations (cyan). Bottom left: the sum of the jp — > 7r°p and 
IP IP as a result of the subtraction of the extracted dummy and scaled e-p elastic 
simulation spectra from the LH 2 AP (black) and from simulations (cyan) from -50 to 
-20 MeV. This window is used for normalization. Bottom right: the LH2 AP (black), 
extracted dummy (red), scaled e-p elastic simulation (blue), and normalized sum of 
the 7p — > 7c°p and yp — > yp from simulations (cyan). 
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Figure 5.30: The right arm yp — > 7r°p and yp — > 7p backgrounds subtraction for 
kinematics b. Top left: the LH 2 AP (black), extracted dummy (red), and scaled e-p 
elastic simulation (blue). Top right: the sum of the yp — > n°p and yp — > •yp data 
backgrounds remaining after the subtraction of the extracted dummy and scaled e-p 
elastic simulation spectra from the LH2 AP (black), and the sum of the 7p — ► 7i°p 
and 7p — > 7p from simulations (cyan). Bottom left: the sum of the jp — > n°p and 
IP — > 7P as a result of the subtraction of the extracted dummy and scaled e-p elastic 
simulation spectra from the LH 2 AP (black) and from simulations (cyan) from -50 to 
-20 MeV. This window is used for normalization. Bottom right: the LH 2 AP (black), 
extracted dummy (red), scaled e-p elastic simulation (blue), and normalized sum of 
the 7£> — > tt°p and yp — > yp from simulations (cyan). 
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5.7.4 Extracting the Reduced Cross Section ctr 

The top left plots of Figures I5.<U| 15. 3 2^ and 15.331 show the individual contributions 
to the AP spectrum for the three representative spectra. In order to extract <tr, 
the ratio of elastic e-p data events to unsealed elastic e-p simulation events R = 
(^elastic-data/^ e -p) * s determined in a AP window around the elastic peak. In 
determining the range of the AP window used for each kinematics, the AP window 
cut was varied to simultaneously minimize the ^-dependence of the dummy and pions 
subtraction, as well as the elastic acceptance correction. The top right plots of Figures 
I5.31| 15. 32^ and 15.331 show the individual contributions in this AP window cut while 
the magenta solid line represents the sum of all contributions. The number of elastic 
events -^elastic-data * n ^ na ^ window is determined as: 

^elastic-data = ^lh 2 ~ Ndummy — N 7p ^ n o p — N 1P -* 1P . (5.15) 



Tables 15.31 and 15.41 list the range of the AP window used and the extracted ratio 
R for each kinematics for both arms. Having determined the ratio R, the unsealed 
elastic e-p simulation is scaled as N^l p = R N e _ p and then added to the remaining 
contributions: 

"^total ^e—p -^dummy -^yp— »7T°p -^yp— >7 > (5.16) 

to generate the sum spectrum shown which reproduces the LH 2 AP spectrum very 
well. The reduced cross section ctr is taken to be the value of the reduced cross section 
from the Bosted fit, a R (Bosted)' usec ^ m ^ ne simulation scaled by the ratio R or: 



a R -Rcr R ( Bos ted) - R 



rG 2 Mp (Q 2 ) + eG 2 Ep (Q 2 ) , (5.17) 
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Figure 5.31: All contributions to the left arm LH2 AP spectrum for kinematics b 
(low e). Top left: the LH2 AP data (black circles), unsealed elastic e-p simulation 
(blue dashed line), extracted dummy (red stars), normalized sum of the 7^ — > 7T°p 
and 7p — > 7p backgrounds (solid green). Top right: the same as in Top left but 
for a smaller AP window cut to show the final AP window used for the final ratio 
extraction. The magenta line through the data is the sum of all contributions. See 
text for more details. Bottom left: the same as in Top right but for a larger AP 
window. Bottom right: the difference between the LH 2 AP data and the sum of all 
contributions. 
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Figure 5.32: All contributions to the left arm LH2 AP spectrum for kinematics m 
(high e). Top left: the LH2 AP data (black circles), unsealed elastic e-p simulation 
(blue dashed line), extracted dummy (red stars), normalized sum of the 7^ — > 7T°p 
and 7p — > 'jp backgrounds (solid green). Top right: the same as in Top left but 
for a smaller AP window cut to show the final AP window used for the final ratio 
extraction. The magenta line through the data is the sum of all contributions. See 
text for more details. Bottom left: the same as in Top right but for a larger AP 
window. Bottom right: the difference between the LH 2 AP data and the sum of all 
contributions. 
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Figure 5.33: All contributions to the right arm LH2 AP spectrum for kinematics 
b. Top left: (black circles) LH2 AP data, (blue dashed line) unsealed elastic e- 
p simulation, (red stars) the modified target walls background, (solid green) the 
modified sum of the jp — > ir°p and 7p — > 7p backgrounds. Top right: the same as 
in Top left but for a smaller AP window cut to show the final AP window needed 
for the final ratio extraction. The magenta line through the data is the sum of all 
contributions. See text for more details. Bottom left: the same as in Top right but 
for a larger AP window cut. Bottom right: the difference between the LH2 AP data 
and the sum of all contributions. 
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Q 2 


Kinematics 


Right AP Window 


Right Ratio 


(GeV 2 ) 




From : To 


R±6 R 






(MeV) 




0.50 





-19.0 


+11.0 


0.9970 ± 0.0024 




a 


-19.0 


+11.0 


1.0021 ± 0.0020 




i 


-18.0 


+12.0 


0.9881 ± 0.0022 




Q 


-18.0 


+12.0 


0.9919 ± 0.0021 




I 


-18.0 


+12.0 


0.9896 ± 0.0028 


0.50 


b 


-19.0 


+11.0 


1.0010 ± 0.0021 




j 


-18.0 


+12.0 


0.9883 ± 0.0019 




V 


-18.0 


+12.0 


0.9921 ± 0.0022 




m 


-18.0 


+12.0 


0.9884 ± 0.0026 


0.50 


k 


-18.0 


+12.0 


0.9865 ± 0.0020 




r 


-18.0 


+12.0 


0.9923 ± 0.0020 




n 


-18.0 


+12.0 


0.9838 ± 0.0029 



Table 5.3: The AP window range used and the extracted ratio R for the right arm 
kinematics. The listed uncertainties in the ratios are statistical only. 



Q' 2 


Kinematics 


Left AP Window 


Left Ratio 


(GeV 2 ) 




From : To 


R±5r 






(MeV) 




2.64 





-12.0 


+2.0 


1.0156 ± 0.0030 




a 


-12.0 


+6.0 


1.0084 ± 0.0026 




i 


-12.0 


+9.0 


1.0054 ± 0.0033 




1 


-13.0 


+11.0 


1.0038 ± 0.0029 




I 


-14.0 


+12.0 


1.0022 ± 0.0026 


3.20 


b 


-12.0 


+3.0 


1.0236 ± 0.0036 




j 


-12.0 


+8.0 


1.0241 ± 0.0033 




V 


-12.0 


+12.0 


1.0152 ± 0.0033 




m 


-14.0 


+14.0 


1.0234 ± 0.0029 


4.10 


k 


-12.0 


+5.0 


1.0296 ± 0.0049 




r 


-15.0 


+10.0 


1.0368 ± 0.0037 




n 


-17.0 


+13.0 


1.0454 ± 0.0040 



Table 5.4: The AP window range used and the extracted ratio R for the left arm 
kinematics. The listed uncertainties in the ratios are statistical only. 
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where the form factors Ge p {Q 2 ) and Gm p (Q 2 ) at a given Q 2 point are determined 
from the Bosted parameterization equations (|2.2j) and (|2.3|) . Tables 1531 and 1531 list 
the reduced cross sections and their uncertainties for each kinematics for both arms. 
Note that the statistical uncertainty in <jr includes the statistical uncertainty in all of 
the yields in equation 15. 151 as well as statistical uncertainty in the extracted dummy 
normalization factor. 

To estimate the size of the uncertainty in the dummy subtraction, we plot the 
ratio of N dummy to N LH2 events or "Dummy Subtraction" in the AP window cut used 
as a function of kinematics number and for both arms. Figures 15.341 and 15.351 show 
the results for the right and left arm, respectively. The dummy subtraction is a ~ 
10% correction for both arms and with a ~ 2% slope dependence for the left arm 
and ~ 0.5% for the right arm. We assign a 5% systematic uncertainty in the dummy 
subtraction correction which yields a slope uncertainty of 0.10% for the left arm 
and 0.36% (0.025%/0.07) for the right arm when we consider the Ae range of 0.07. 
Furthermore, a 5% uncertainty in a 10% correction yields a scale uncertainty of 0.5% 
for the both arms. The random uncertainty is set to 0.0% to avoid double counting 
the random uncertainty in the dummy subtraction since the statistical uncertainty 
in the effective dummy thickness as shown in Figure 15.271 is taken into account when 
calculating the statistical uncertainty 5r in the ratio R for both arms. See Tables 1531 
and 1531 

Similarly, to estimate the size of the uncertainty in the pion subtraction, we plot 
the ratio of ( N 1P ^o p + N 1P ^ 1P ) to Nlh 2 events "Pion Subtraction" in the AP window 
cut used as a function of kinematics number and for both arms. Figure 15". 361 shows the 
results for the left arm. The pion subtraction is typically a ~ 0.5% correction for the 
left arm with a large scatter and e dependence. We assign a ~ 50% uncertainty in the 
pion subtraction correction or 0.25% with almost equal contribution from the slope 
(0.20%) and random (0.15%) uncertainties. The slope and random uncertainties are 
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Q 2 


Kinematics 


e 


^ Stat 


^Random 






(GeV 2 ) 














0.50 





0.914 


0.054E-02 


0.104E-02 


0.117E-02 


0.22856 




a 


0.939 


0.046E-02 


0.149E-02 


0.156E-02 


0.23253 




i 


0.962 


0.051E-02 


0.106E-02 


0.118E-02 


0.23187 




Q 


0.979 


0.051E-02 


0.107E-02 


0.119E-02 


0.23468 




I 


0.986 


0.066E-02 


0.108E-02 


0.127E-02 


0.23491 


0.50 


b 


0.939 


0.048E-02 


0.148E-02 


0.156E-02 


0.23231 




j 


0.962 


0.045E-02 


0.106E-02 


0.115E-02 


0.23192 




V 


0.979 


0.052E-02 


0.107E-02 


0.119E-02 


0.23492 




m 


0.986 


0.063E-02 


0.108E-02 


0.125E-02 


0.23462 


0.50 


k 


0.962 


0.047E-02 


0.105E-02 


0.115E-02 


0.23150 




r 


0.979 


0.047E-02 


0.107E-02 


0.117E-02 


0.23477 




n 


0.986 


0.070E-02 


0.107E-02 


0.127E-02 


0.23353 



Table 5.5: The elastic e-p reduced cross section cr R for the right arm kinematics. 
In addition, the statistical uncertainty 5stat, random systematic uncertainty SR an dom, 
and total uncertainty 5 aR (combined 5stat and 8R an dom) in &r are also given. 



Q 2 


Kinematics 


e 


<5 Stat 


^Random 


S <TR 


ctr 


(GeV 2 ) 














2.64 





0.117 


0.40E-04 


0.51E-04 


0.65E-04 


0.13402E-01 




a 


0.356 


0.35E-04 


0.82E-04 


0.89E-04 


0.13823E-01 




i 


0.597 


0.47E-04 


0.55E-04 


0.72E-04 


0.14301E-01 




1 


0.782 


0.44E-04 


0.57E-04 


0.72E-04 


0.14677E-01 




I 


0.865 


0.38E-04 


0.58E-04 


0.69E-04 


0.14832E-01 


3.20 


b 


0.131 


0.30E-04 


0.42E-04 


0.52E-04 


0.08686E-01 




j 


0.443 


0.29E-04 


0.35E-04 


0.45E-04 


0.09075E-01 




P 


0.696 


0.30E-04 


0.36E-04 


0.47E-04 


0.09305E-01 




m 


0.813 


0.27E-04 


0.37E-04 


0.46E-04 


0.09525E-01 


4.10 


k 


0.160 


0.23E-04 


0.22E-04 


0.31E-04 


0.04726E-01 




r 


0.528 


0.17E-04 


0.23E-04 


0.29E-04 


0.04973E-01 




n 


0.709 


0.19E-04 


0.24E-04 


0.31E-04 


0.05121E-01 



Table 5.6: The elastic e-p reduced cross section &r for the left arm kinematics. In 
addition, the statistical uncertainty Sstat, random uncertainty 5 Random-, and total un- 
certainty 5 aR (combined 5 S tat and 5 Random ) in a R are also given. 
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Figure 5.34: The ratio of N dummy to -/V £ # 2 events or "Dummy Subtraction (%)" for the 
right arm in the AP window cut used as a function of kinematics number. The points 
are sorted according to the kinematics of the left arm. Kinematics 1-5 correspond 
to Q 2 = 2.64 GeV 2 and shown in red, while kinematics 6-9(10-12) correspond to 
3.20(4.10) GeV 2 and shown in black(blue). For each Q 2 value, the points are sorted 
by e (low to high). 
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Figure 5.35: The ratio of N dummy to N LH2 events or "Dummy Subtraction (%)" for the 
left arm in the AP window cut used as a function of kinematics number. Kinematics 
1-5 correspond to Q 2 = 2.64 GeV 2 and shown in red, while kinematics 6-9(10-12) 
correspond to 3.20(4.10) GeV 2 and shown in black(blue). For each Q 2 value, the 
points are sorted by e (low to high). 
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doubled for Q 2 = 4.10 GeV 2 to account for the larger observed e dependence. For the 
right arm, the pion subtraction correction is set to zero and a random uncertainty 
of 0.05% is assigned. It is crucial to mention that when determining the ratio R 
for the right arm, the number of elastic events -^elastic-data as gi ven by equation 
(I5.15|) is now given by Elastic-data = Nlh 2 ~ N dummy . This is equivalent to set the 
normalization factor or the 7r°-factor needed to normalize the sum of jp —>■ n°p and 
IP ~ * IP from simulations to the resultant (LH 2 - endcaps - scaled e-p simulation) 
spectrum to zero. The statistical uncertainty in the 7r°-factor is taken into account 
when calculating the statistical uncertainty 5r in the ratio R for the left arm. See 
section 15.7.31 for details. 
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Figure 5.36: The ratio of {N lp ^o p + N 1P ^ W ) to Nlh 2 events or "Pion Subtraction 
(%)" for the left arm in the AP window cut used as a function of kinematics number. 
Kinematics 1-5 correspond to Q 2 = 2.64 GeV 2 and shown in red, while kinematics 
6-9(10-12) correspond to 3.20(4.10) GeV 2 and shown in black(blue). For each Q 2 
value, the points are sorted by e (low to high). 

To estimate the uncertainty due to the AP cut used, we plot the percentage of 
the N e _ p events or "Elastic Peak Acceptance" within the AP window cut used as a 
function of kinematics number and for both arms. Figures 15.371 and 15.381 show the 
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results for the right and left arm, respectively. The elastic acceptance is ~ 85% and 
~ 95% for the left and right arm, respectively. Of the 15% events lost on the left arm, 
~ 14% are events lost in the radiative tail and ~ 1% are events lost in the smeared 
non-gaussian tails of the elastic peak. The uncertainty in the AP cut used due to 
the correction that results from the events lost in the radiative tail is accounted for 
in the radiative corrections uncertainties. See section 15.61 for details. The smeared 
non-gaussian tails are known to within ~ 50%, and therefore, a 1% correction due to 
events lost in the non-gaussian tails yields a 0.50% scale uncertainty for the left arm. 
We vary the AP window cut used on the left arm by ± 2 MeV which changes the ratio 
R by ~ 0.20%. The 0.20% uncertainty in R is then broken down equally as 0.14% 
random uncertainty and 0.14% slope uncertainty. However, if we consider the average 
Ae range of 0.70, that will modify the slope uncertainty to 0.20% (0.14%/0.70). 

The right arm is more complicated than the left arm since the HRS momentum 
acceptance cuts off the AP peak. The loss due to a non-gaussian tails is ~ 1%, and 
we take the same error estimates as for the left arm but account for the Ae range of 
0.07 which modifies the slope uncertainty to 2.0% (0.14%/0.07). 

It is crucial to realize that although the elastic peak is broadened with increasing e 
for all the kinematics on the left arm, it is well within the momentum acceptance 5 of 
the HRS and thus no elastic events are lost after applying the momentum acceptance 
5 cut. On the other hand, this is not the case for the right arm since the elastic peak is 
broader than that of the left arm and extends a bit beyond the edges of the momentum 
acceptance. Therefore, applying the § cut will get rid of some elastic events. The 
elastic acceptance for the right arm shows a ~ 1.5% e dependence as can be seen 
from Figure 15.371 Varying the 5 cut changes the correction by ~ 30%. Therefore, 
a 30% uncertainty in a 1.5% correction yields a ~ 0.5% uncertainty. However, the 
30% uncertainty is not simply a scale, random, or slope uncertainty and so we break 
the 0.5% uncertainty equally among the scale, random, and slope uncertainties and 



215 




2 4 6 8 10 
Kinematics Number 



Figure 5.37: The percentage of the events or "Elastic Peak (%)" for the right arm in 
the AP window cut used as a function of kinematics number. The points are sorted 
according to the kinematics of the left arm. Kinematics 1-5 correspond to Q 2 = 2.64 
GeV 2 and shown in red, while kinematics 6-9(10-12) correspond to 3.20(4.10) GeV 2 
and shown in black(blue). For each Q 2 value, the points are sorted by e (low to high). 
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Figure 5.38: The percentage of the N e _ p events or "Elastic Peak %" for the left 
arm in the AP window cut used as a function of kinematics number. Kinematics 
1-5 correspond to Q 2 = 2.64 GeV 2 and shown in red, while kinematics 6-9(10-12) 
correspond to 3.20(4.10) GeV 2 and shown in black(blue). For each Q 2 value, the 
points are sorted by e (low to high). 
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assign a ~ 0.30% uncertainty for each. Finally, if we consider the Ae range of 0.07 
for the right arm, that will modify the slope uncertainty to 4.28%. 



Chapter 6 Results and Two-Photon 
Exchange 



6.1 Overview 

In this chapter I will discuss the procedure used to extract the electromagnetic form 
factors of the proton, Ge p and Gmp, and their ratio M ^P Ep . Rosenbluth extractions 
of the proton form factors and their ratio using the reduced cross sections of each 
arm separately will be discussed. Rosenbluth extractions using the left arm reduced 
cross sections were recently published |lU7j . Initially, we planned to use the Super- 
Rosenbluth ratio, Rsuper, defined as the left arm reduced cross section normalized to 
the right arm ratio at each kinematics point or: 

Left ctr 

Rsuper = RighTR ' (0) 

to extract the form factor ratios. Note that by using Rsuper, we take advantage of 
the right arm spectrometer which served as a luminosity monitor to remove any un- 
certainties due to beam charge, current, target length, raster size, and target density 
fluctuations. See chapter 0] for a detailed discussion of the corrections and efficiencies 
applied to the measured beam charge, and Tables 1%. 181 and fa.lDl for a summary of the 
systematic uncertainties in these corrections and efficiencies. However, we decided 
not to use Rsuper to extract the form factors and their ratio since it is estimated that 
by dividing by the luminosity monitor, that will introduce more uncertainty in Rsuper 
than it will eliminate. By looking at Table 1531 we see that the right arm cross sections 

have ~ 0.50% random uncertainty, while the luminosity-related random uncertainty 
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is < 0.20% as can be seen from Table E~TS1 Instead we use the right arm results as a 
consistency check of the assumed uncertainties. 

The form factors and their ratio extracted from this work are discussed and com- 
pared to previous Rosenbluth extractions and recent polarization transfer results. 
See section 12.21 and subsections therein for a full discussion of previous Rosenbluth 
extractions and recent polarization transfer results. An introduction to the two- 
photon-exchange (TPE) and Coulomb distortion corrections as possible sources for 
the discrepancy between the Rosenbluth extractions and polarization transfer results 
is given. The linearity in the Born approximation of the Rosenbluth separation plots 
(L-T plots) as a function of e is examined for this work because any deviation from 
linearity would provide a clear signature of the effects beyond traditional radiative 
corrections, and would provide information on the nonlinear component of the TPE. 
Finally, as the results of this work raised the interest in the physics of the TPE and 
Coulomb distortion corrections and laid down the foundations for new experiments 
aimed at measuring the size of the TPE corrections, I conclude with a brief discussion 
of the TPE calculations and presentations of these future experiments. 



6.2 Form Factors Extraction 

In this section I will describe the general procedure used to extract Ge p , Gm p , and 
^p Gep by using the reduced cross sections or, as listed in Tables 1531 and lo~o1 

6.2.1 Form Factors Extraction Using ctr (Single Arm) 

Equation (jl.H4j) defines the elastic e-p reduced cross section ctr. A linear fit of ctr 
to e at a fixed Q 2 gives tG 2 Mp (Q 2 ) as the intercept and G 2 Ep (Q 2 ) as the slope. If the 
goal is to extract Gm p and Ge p , then a linear fit of ctr to e is performed using the 
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following form: 

a R = T A 2 + B 2 e , (6.2) 

where A and B are the parameters of the fit and represent Gm p and Ge p , and the 
uncertainties in A and B represent SGm p and SGe p , respectively. On the other hand, 
if the goal is to extract 77^-, then a linear fit of (Tr to e is performed using: 



*« tA-\ 1.0 + B 2 ^ ) . (6.3) 



where A and B represent Gm p and off-, and the uncertainties in A and B represent 



Gep 

SGm p and ^off-, respectively. Note that the ratio and its uncertainty can be 
determined first by taking the ratio of the individual form factors as extracted using 
equation ()6.2|) and then propagating the (correlated) errors on Ge p and Gm p , however, 

Gep 
Gm p 



equation ()6.2|) is used since it directly provides the ratio and its uncertainty. 



(a) Ge p and Gm p Extraction Using Equation 16.21 

The classification of the systematic uncertainties into scale, random, and slope uncer- 
tainties and their effects on the extracted form factors has been discussed in section 
14.11 To determine the individual form factors, we start by adding in quadrature the 
total random uncertainty in <tr, 5n an dom, to the statistical uncertainty in each <jr, 
Sstat, as listed in Tables 1531 and I5~o1 or: 



5 *R = V 6 Stat + ^Random " ( 6 ' 4 ) 



A linear fit of to e at a fixed Q 2 using equation (J6.2|) yields a value of 
Gm p = Gj^i^G]^ and Ge p = ±5G^ where the contribution from the random 
uncertainty to the total uncertainty in Gm p and Ge p is given by <5Gf/™ dom = SG^ 
and 5G^ ndom = 5G$ p . Figures EU E3 ESI and EM show such fits done at Q 2 = 
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0.50, 2.64, 3.20, and 4.10 GeV 2 , respectively. 

To determine the contribution of the slope uncertainty in or to the total uncer- 
tainty in Gmp and Ge p , we vary ctr by the uncertainties that vary linearly with e or 
the total slope uncertainty, 6si ope , as given at the bottom of Tables l4~T51 and 14. 191 



a 



Slope 
R 



G R 1.0 + E5, 



Slope 



(6.5) 



then, a linear fit of a s R lope to e using equation ()6.2|) . with uncertainties in a R ope given 
by 8 c sio P e = 5 aR (1.0 + e5 siope), yields a value of G Mp = Gfj p ±5G^ p and G Ep = G% p ± 

(2) 

SG Ep . Therefore, the contribution from the slope uncertainty to the total uncertainty 
in Gm p and Ge p is given by SG^ ope 



^Mp 



^Mp 



and 5G E l ° pe 



^Ep 



r (2) 
^Ep 



Note that for the right arm we used 5 sio Pe 
measurements were taken near e=1.0. 



5or R ^1.0 + (e — l-0)Ssiopej since all the 
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Figure 6.1: A linear fit of the right arm or to e at Q 2 = 0.50 GeV 2 using equation 
(16. 2|) . The shown uncertainties are the combined statistical and random from equation 

dm. 
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Figure 6.2: A linear fit of the left arm or to e at Q 2 = 2.64 GeV 2 using equation ()6.2|) . 
The shown uncertainties are the combined statistical and random from equation (|6.4|) . 




Figure 6.3: A linear fit of the left arm cxr to e at Q 2 = 3.20 GeV 2 using equation ()6.2|) . 
The shown uncertainties are the combined statistical and random from equation ([6.4)1 . 
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Figure 6.4: A linear fit of the left arm cxr to e at Q 2 = 4.10 GeV 2 using equation ()6.2|) . 
The shown uncertainties are the combined statistical and random from equation (j6.4j) . 

Similarly, we determine the contribution of the scale uncertainty in <jr to the total 
uncertainty in Gm p and Ge p by varying a R by the total scale uncertainty, Ss ca ie, as 
given at the bottom of Tables 14.181 and 14.191 



a 



Scale 
R 



a R \ 1.0 + 5 



Scale 



(6.6) 



and a linear fit of af^ ale to e using equation (j6.2j) . with uncertainties in a§£ ale given by 
^scaie = ^(l.O+^ca/e), yields a value of G A / P = G^ p ±SG^ p and G Sp = G^> p ±8G^> p . 
The contribution from the scale uncertainty to the total uncertainty in Gm p and Ge p 
is therefore given by SGf^ le 



^Mp ^Mp 



.(3) 



and 5G s E c p ale 



G 



Finally, G Mp and G Bp are given by G Mp = G { ^ p ±5G\f p al and G Bp = G^itfG'^ 



yTotal 



yTotal 
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where 5G\° t p al and 5G^ al are given by: 



Ep 



XfiTotal 
0Lr Mp 



^QRandomy + (SG^) 2 



{5G S ^ 2 



(6.7) 



and 



\j (SG§p ndom ) 2 + (6G s J; p y + (SG% ale y . 



(6.8) 



The electric and magnetic form factors for the proton with their uncertainties 
were determined at Q 2 = 0.50, 2.64, 3.20, and 4.10 GeV 2 . Tables O and O list the 
results. 



Q 2 

(GeV 2 ) 


Gup 


X/^Random 
0Ljr Mp 


cs-iSlope 
0Lr Mp 




0Lr Mp 


± 

S( ° Mp ) 


0.50 


1.03340 


0.06726 


0.01247 


0.01702 


0.07049 


1.07467 
± 

0.07330 


2.64 


0.1325 


0.0003466 


6.988E-06 


0.001913 


0.001944 


1.05632 
± 

0.01549 


3.20 


0.09691 


0.0003204 


7.398E-06 


0.001398 


0.001435 


1.05233 
± 

0.01558 


4.10 


0.06291 


0.0002822 


4.701E-06 


0.0009097 


0.0009525 


1.03390 
± 

0.01565 



Table 6.1: The magnetic form factor for the proton Gm p at Q 2 = 0.50, 2.64, 3.20, 
and 4.10 GeV 2 as determined by fitting a R to e using equation ()6.2|) . In addition, the 
contributions from the random, slope, and scale uncertainties to the total uncertainty 
in Gmp is also given. 
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Q 2 

(GeV 2 ) 


Gep 


^^yRandom 


5G E l ; pe 


X/^ Scale 
O^Ep 


xr^Total 
0<^Ep 


± 


0.50 


0.2906 


0.03589 


0.02682 


0.003182 


0.04492 


0.8440 
i 

± 
0.1304 


2.64 


0.04385 


0.001276 


0.0009129 


0.0006389 


0.001694 


0.9761 
± 

0.03772 


3.20 


0.03436 


0.001377 


0.0007532 


0.0004959 


0.001646 


1.04198 
± 

0.04994 


4.10 


0.02666 


0.001489 


0.0007684 


0.0003775 


0.001718 


1.2238 
± 

0.07886 



Table 6.2: The electric form factor for the proton Ge p at Q 2 = 0.50, 2.64, 3.20, and 
4.10 GeV 2 as determined by fitting a R to e using equation (|6.2|) . In addition, the 
contributions from the random, slope, and scale uncertainties to the total uncertainty 
in Ge p is also given. 



G 

(b) G Mp and g| Extract ion Using Equation 16.3 



A linear fit of or to e at a fixed Q 2 using equation ()6.3|) yields Gm p and -§ff-- The 
same procedure outlined in subsection (a) is used with the following changes: 

• replaces Ge p whenever it occurs or: Ge p — > ~§ff~- 

Therefore G Mp = G { l} p ± and = (§fj) {1) ± <K§J)™ where 

and 5 (-^ L ) Total are given by: 



SG T M t l = \ {5GZ dom ) 2 + (SG s ^y + (5Gfr p le ) 2 , (6.9) 
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and 



6( 



Ge p sTotal 



G 



Mp 




£l_ \ Random j 



v t*Mp J v ^Mp 



\ Scale 



(6.10) 



The magnetic form factor and the ratio of electric to magnetic form factor for the 
proton with their uncertainties were determined at Q 2 = 0.50, 2.64, 3.20, and 4.10 
GeV 2 . The results obtained for Gup are identical to those listed in Table I6~T1 Table 

G 

lists the results obtained for the ratio -p^ 2 -. 

M p 



Q 2 

(GeV 2 ) 




r/ Ge p \ Random 
^ Gm p ' 


XI Ge p \ Slope 
\ Gm p ' 


5(c Ep ) Scale 


5 ,G Ep , Total 
y Gm p i 


P^pGEp 
Gm p 
± 

y GMp ' 


0.50 


0.2803 


0.05303 


0.03039 


0.000655 


0.06112 


0.7828 

± 
0.1707 


2.64 


0.3308 


0.01040 


0.006930 


6.202E-05 


0.01249 


0.9240 
± 

0.03490 


3.20 


0.3545 


0.01529 


0.007770 


4.100E-05 


0.01715 


0.9900 
± 

0.04789 


4.10 


0.4238 


0.02540 


0.01222 


6.601E-05 


0.02819 


1.1837 
± 

0.07873 



Table 6.3: The ratio of electric to magnetic form factor for the proton at Q 2 = 
0.50, 2.64, 3.20, and 4.10 GeV 2 as determined by fitting a R to e using equation (|6.3|) . 
In addition, the contributions from the random, slope, and scale uncertainties to the 
total uncertainty in ^ is also given. 
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6.3 Discussion of the Results 
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Figures 16.51 shows the world data on ^q^ p ratio for the proton from previous Rosen- 
bluth measurements along with the results of E01-001 single arm ( Table l6~3j) . In 
addition, Figures EH shows the ratio from the E01-001 single arm extractions 

along with the results of recoil polarization. The global analysis of previous cross sec- 
tion data by Arrington is also shown. By comparing the results of the EO 1-001 
with those of previous Rosenbluth separations and recoil polarization measurements, 
we conclude that: 

1. The values of — ratio from the single arm extractions are in 
good agreement with the previous Rosenbluth measurements but with 
much smaller statistical and systematic uncertainties than the pre- 
vious best measurements [15]. The total uncertainty is comparable 
to those of recoil polarization. 

2. The high precision achieved in the E01-001 extractions: 

• Leaves little room for doubting that the values of ratio ex- 
tracted from Rosenbluth separations are inconsistent with those 
reported from the high-Q 2 recoil polarization measurements. 

• Makes it clear that the problem is not simply an experimental 
error in the Rosenbluth measurements or technique. 

• Confirms the reliability of the elastic e-p cross sections ex- 
tracted from previous Rosenbluth separations. 

3. (This point will be discussed below): 

Provides a significant test of the validity of the commonly used ra- 
diative corrections (bremsstrahlung corrections). 
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Figure 6.5: The world data of Rosenbluth separation determination of \i v Ge v IGm p 
ratio compared to that of E0 1-001 (single arm extraction) (Table 163)1 . 
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Figure 6.6: The values of h p Ge p /Gm p for the proton from E01-001 (single arm ex- 
traction) (Table I6.3)l and recoil polarization. 
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In summary, we see that the results of the E01-001 are in agreement with the previ- 
ous Rosenbluth data. The experiment provided systematic uncertainties much smaller 
than the best previous Rosenbluth measurements, ^H], and comparable to those of 
the recoil polarization clearly establishing the discrepancy between the Rosenbluth 
separations and recoil polarization results. These results confirmed that the discrep- 
ancy is not due to experimental error in the Rosenbluth measurements and provided 
a strong test of the conventional radiative corrections used. 

A possible source for the discrepancy is the effect of higher-order corrections to 
the reduced cross sections. Suggested sources for the discrepancy will be discussed 
in section 16.41 The reduced cross sections as defined by equation (|1.34|) and listed in 
Tables 15.51 and 15.61 are derived in the Born approximation and have been corrected 
for many of the higher-order radiative corrections terms. The form factors and thus 
their ratio have been extracted using these reduced cross sections without including 
fully the effects of these higher-order corrections on the reduced cross sections for 
consistency with previous Rosenbluth measurements. 

The standard radiative correction procedure used was discussed in section 15.61 
We correct the data for bremsstrahlung, vertex corrections, and loop diagrams, with 
almost all of the e dependence coming from bremsstrahlung. Figure 13.61 shows the 
calculated radiative correction factor (internal corrections only) as a function of e for 
Q 2 = 2.64 GeV 2 . Again, while the magnitude of the corrections is similar and both 
show an approximately linear dependence on e, the dependence is much smaller for 
protons, ~ —8%, than for the electron ~ 17%, and of opposite sign. Note that when 
detecting the electron, the e dependence of the bremsstrahlung correction exceeds the 
e dependence coming from Ge p which is not the case in the EO 1-001 experiment. 
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6.4 Possible Sources for the Discrepancy 



The results of the EO 1-001 experiment have clearly established and confirmed the 
discrepancy between the Rosenbluth separations and recoil polarization results. The 
questions that need to be answered are: why do the two techniques disagree? and 
which form factors are the correct ones? 

As noted above, we have essentially ruled out the possibility that the discrep- 
ancy is due to an experimental error in the Rosenbluth separations measurements or 
technique. That confirms that previous elastic e-p cross sections are reliable which is 
crucial since these cross sections are used for normalization or as an input to the anal- 
ysis of different experiments. In addition, reanalysis of the existing form factor data 
by Arrington [HH], see section l2~H for details, has confirmed that all the different data 
from previous Rosenbluth separations measurements are consistent. This reanalysis 
combined with the new results from E01-001 rule out any possibility for a bad data 
sets or incorrect normalization in the combined Rosenbluth analysis. It should be 
noted that while the E01-001 experiment provided an independent check 
on previous Rosenbluth measurements there has been, to the best of my 
knowledge, no independent check of the high-Q 2 recoil polarization results. 

As for which form factors are the correct ones to use, again, if the cross section 
measurements and hence the Rosenbluth extractions are incorrect, that still will not 
solve the problem since the recoil polarization technique provides the ratio of Ge p 
to Gmp arid not the actual values for the individual form factors. While the form 
factors extracted by Rosenbluth separations may possibly not represent the true form 
factors for the proton, they do indeed provide the correct elastic e-p cross sections and 
presumably provide the best parameterization when elastic cross sections are used for 
normalization or as an input to the analysis of different experiments. 
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6.4.1 Two-Photon-Exchange (TPE) Correction 

The difference in the ratio from Rosenbluth separations and recoil polariza- 

tions results can be explained by a common e-dependent systematic error in the cross 
section measurements of (5-8)% j22 HUH DUE] ■ The establishment of the discrepancy 
has led people to suggest that the discrepancy is due to a missing corrections in the 
cross section due to a higher-order processes such as two-photon-exchange (TPE). 
The standard radiative correction procedure used was discussed in section 15.61 Fig- 
ure 15.151 shows the full contribution of the elastic and inelastic radiative corrections 
to the elastic e-p cross section including the first-order QED (Born) radiative correc- 
tions or single-photon-exchange. Figure l5~T5l (e). crossed-box, and (f), box, are the 
contributions from TPE to the elastic e-p cross section which are not fully included 
in the standard radiative correction procedures. See |109t 11041 1105] for details. It is 
believed that inclusion of TPE contributions may remove the discrepancy. 

If we assume that the discrepancy is in fact due to TPE and the recoil polarization 
results are approximately correct, then the ratio becomes small enough as 

Q 2 gets large (note the smaller slope predicated by recoil polarization) that the e 
dependence of the reduced cross section is dominated by TPE. Figures 16.71 and 16.81 
show the contribution of the TPE, A2 7 , to e dependence of the L-T plots at Q 2 
= 2.64 and 4.10 GeV 2 , respectively. At large Q 2 where Ge p is suppressed and the 
reduced cross section in the Born approximation has almost no e dependence, the 
discrepancy implies cross section difference of 5-10%. In this region, almost all of 
the e dependence would come from TPE and the e dependence of the reduced cross 
section will represent the full contributions of TPE. 

The high precision provided by the E01-001 experiment on <jr and hence 
on is of great importance as it constitutes a precise measurement of 

the discrepancy between the Rosenbluth separations and recoil polariza- 
tion results, thus allowing for a precise test as to whether TPE effects can 
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Figure 6.7: The e dependence of the reduced cross section at Q 2 = 2.64 GeV 2 as 
measured in the EO 1-001 experiment (blue solid circles), and as predicted from the 
recoil polarization (red dashed line). The black solid line is a linear fit to the data. If 
the recoil polarization results represent the true form factors (note the smaller slope 
or G\ v predicted by the recoil polarization and shown in red double arrow), then TPE 
or Ao"2 7 shown as blue double arrow will yield roughly about half of the e dependence 
at Q 2 = 2.64 GeV 2 . 
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Figure 6.8: The e dependence of the reduced cross section at Q 2 = 4.10 GeV 2 as 
measured in the EO 1-001 experiment (blue solid circles), and as predicted from the 
recoil polarization (red dashed line). The black solid line is a linear fit to the data. If 
the recoil polarization results represent the true form factors (note the smaller slope 
or G\ predicted by the recoil polarization and shown in red double arrow), then 
TPE or Ao"2 7 shown as blue double arrow will yield about 85% of the e dependence 
at 4.10 GeV 2 . 
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explain the difference between Rosenbluth and recoil polarization mea- 
surements. The reduced cross section in the Born approximation varies linearly 
with e. Therefore, at low Q 2 where Ge p contributes most, any deviation from linear- 
ity would have to come from higher-order terms that are not included in the standard 
radiative correction procedures. Such deviation from linearity would provide a clean 
signature of TPE and would provide information about the nonlinear component of 
TPE. If the TPE correction is to explain the discrepancy, it would have to introduce 
(5-8)% e-dependent correction to the cross sections. Such correction would have to 
be roughly linear since large nonlinearities are not observed. 

An important point to emphasize is that the standard radiative corrections 
have an e dependence that is comparable to the slope brought about by the 
form factors, (see Figure l3.6|) . In addition, if missing radiative correction terms 
are the reason for such discrepancy, then the form factors extracted from Rosenbluth 
separations and to lesser extent those from recoil polarization do not represent the 
true form factors for the proton since we are not correctly isolating the single-photon- 
exchange process. Therefore, it is absolutely crucial to understand the radiative 
corrections and in particular higher-order processes such as TPE theoretically and 
experimentally before we can be confident in our knowledge of the proton form factors. 

Several calculations in the 1950s and 1960s tried to estimate the size of the TPE 
contributions to the unpolarized elastic e-p cross sections |9T| 1951 1991 1100(1101] . While 
some calculations used only the proton intermediate state [101 j . other included the 
excited intermediate states of the proton as well |97|I581I59|I1UU| . The TPE corrections 
estimated from these calculations were extremely small (< 1%) and were not included 
in the standard radiative correction procedures. 

Motivated in part by the results of the E0 1-001 experiment, the interest in the 
physics of TPE and Coulomb distortion corrections and their impact on electron 
scattering observables has recently risen both experimentally and theoretically. In 
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addition, the results of the EO 1-001 experiment have laid down the foundation for 
new experiments aimed at measuring the size of the TPE corrections (to be discussed 
in section 16. 5 j) . All of the new calculations have predicted a nonlinearities in the 
e dependence of TPE. A brief description of these recent calculations will be given 
below. In addition, their effect on the ratio from the EO 1-001 extractions will 

be examined where possible. 

6.4.2 Recent Calculations of the TPE Corrections 

• Calculation by Blunden, Melnitchouk, and Tjon et al. \104)l : 

This is a low energy hadronic model that takes into account the proton inter- 
mediate state and neglects excited intermediate states. It includes only the elastic 
contributions to the TPE correction or the box and crossed-box diagrams shown in 
Figure UTTol (f) and (e), respectively. Their initial calculations yielded an e depen- 
dence of ~ 2% with small nonlinearities at small e and insignificant Q 2 dependence. 
These corrections became larger with the inclusion of an improved form factors 105J . 
Figure 16.91 shows the effect of such a calculation on the ratio ^ Ggp extracted from 
the EO 1-001 experiment. In addition, the effect of Coulomb correction by Sick and 
Arrington |110j (see section 16.4.31 for details) is also shown. While the discrepancy is 
largely resolved for Q 2 = 2-3 GeV 2 , this calculation fails to resolve the discrepancy for 
Q 2 > 3 GeV 2 . Recently, the authors have included the A resonance as an intermedi- 
ate excited state to their elastic box and crossed-box calculations of TPE corrections 
106 . The A contribution is smaller in magnitude than the contribution when the 
proton is in the intermediate state with a small modification to the e dependence of 
the TPE corrections. 
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Figure 6.9: The \i v Ge v IGm p ratio for the proton from the E01-001 experiment (red 
open circles) and recoil polarization (black open diamonds). The EO 1-001 \i v Ge v IGm v 
corrected for TPE correction by Blunden et al 1 1 04] (blue solid circles). The E01- 
001 UpGep/Gmp ratio corrected for TPE correction by Blunden et al including the 
Coulomb correction by Sick and Arrington |110j (green crosses). 



• Calculation by Chen et al. [Ill]: 

This is a high energy model at the quark-parton level. It calculates the TPE 
correction using a generalized parton distribution (GPD) for the quark distribution 
to describe the emission and re-absorption of the partons by the nucleon as shown 
in Figures [6.101 (the dominant contribution) and 16 .111 For large Q 2 , the calculation 
shows a significant e dependence and nonlinearity to the correction with a weak Q 2 
dependence. Since this is a high energy model, it is not expected to be valid at low 
Q 2 values. Figure 16.121 shows the effect of such a calculation on the ratio from 

the E0 1-001 experiment. In addition, the effect of Coulomb correction by Sick and 
Arrington is also shown. Their correction predicts ~ 50% of what is needed to fully 
resolve the discrepancy. Their recent result |112j also shows the model dependence 
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due to the use of different GPDs. Using a different GPD model can significantly 
enlarge their corrections. 




Figure 6.10: Feynman diagram known as the handbag diagram used for the partonic 
calculation of TPE contribution to elastic e-p scattering by Chen et al The 
electron (green) interacts with the proton via two virtual photons (red), while the 
proton (blue) shown as solid line and circle to represent the proton GPDs emits quark 
(black) which interacts with the two virtual photons. The quark gets re-absorbed 
again by the proton. 




Figure 6.11: Feynman diagram known as the cat's ears diagram used for the partonic 
calculation of TPE contribution to elastic e-p scattering by Chen et al |lllj . The 
electron (magenta) interacts with the proton via two virtual photons (red), while 
the proton (blue) shown as solid line and circle to represent the proton GPDs emits 
two quarks qi and qi (black and green solid lines) which interact with the two virtual 
photons. The two quarks get re-absorbed again by the proton (black and green dashed 
lines) . 
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Figure 6.12: The u p Ge p /Gmp ratio for the proton from the E01-001 experiment (red 
open circles) and recoil polarization (black open diamonds). The EO 1-001 \i v Ge v IGm v 
ratio corrected for TPE correction by Chen et al |lllj (blue solid circles). The EO 1-001 
H p Gep/Gmp ratio corrected for TPE correction by Chen et al including the Coulomb 
correction by Sick and Arrington (green crosses). 



• Calculation by Afanasev et al. [113]: 

This calculation is done at the quark-parton level in the double logarithm approx- 
imation. This calculation predicts only the e and Q 2 dependence but not the overall 
magnitude, and therefore cannot be applied as a correction to the data without ad- 
ditional assumptions. In addition, this calculation yields large nonlinearities at large 
e where the e dependence of the TPE corrections behaves in the opposite direction 
to that of Blunden fTUIj. 

• Predictions by Rekalo and Tomasi-Gustafsson et al. [114l : 

This work predicts the TPE correction based on symmetry arguments. The re- 
sults show large nonlinearities at large e. They predict that TPE corrections should 



depend on x = y/1 + ej y/1 — e, but do not calculate the size or Q 2 dependence, 
calculation predicts a similar e dependence to that of Afanasev |112j . 
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Clearly, if TPE explains the discrepancy it must be due to an e dependence mod- 
ification of the reduced cross sections and thus the overestimate of Ge p in the Rosen- 
bluth separations. But this is not the whole story since TPE can have an effect on 
Gm p as well. Most of the previously described calculations predicated that the largest 
effect on the cross section would occur at small e. This would reduce <7r at e = 0.0 
(tG 2 Mjj ) by 3-5% at large Q 2 with small Q 2 dependence. Therefore, it is crucial to 
know the exact e dependence of the TPE correction and in particular the size of the 
nonlinearity as e — >-0 in order to extract an accurate value of Gm p - 

6.4.3 Multiple Soft Photon Exchange (Coulomb Distortion) 

It has been reported that inclusion of the Coulomb distortion to the elastic e-p cross 
section can slightly reduce the discrepancy between the Rosenbluth separations and 
recoil polarization results |1 15| IllOj . Higinbotham |115j evaluated the effect in the 
Effective Momentum Approximation (EMA) where the incident electron with energy 
E approaches the nucleus is accelerated due to its Coulomb interaction with the 
nucleus. Consequently, the four momentum transfer Q gets modified and the nucleus 
is probed at a slightly larger Q, ^Modified 1 



will affect the reduced cross sections such as it will raise up the reduced cross sections 
more at low e values than at high e. This in turn will reduce the Rosenbluth extracted 




(6.11) 



where R eq is the hard sphere equivalent radius of the proton. Therefore, ^Modified 
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slope, G E , and hence lower the form factors ratio. The Coulomb distortion correction 
to the elastic e-p cross section exhibits a small e dependence of 0.5% which is ~ 
10% of the effect needed to bring the Rosenbluth and recoil polarization results into 
agreement. 

The Coulomb distortion and in particular the effect of the proton charge on the 
ingoing and outgoing electron waves has been examined using the second order Born 
approximation based on the approach used in electron-deuteron scattering |110j . Such 
distortions have usually been ignored for Z = 1, but have a significant effect on the 
proton rms radius. An exponential charge distribution for the proton was assumed for 
the calculation using the fact that Ge p ~ Gjj where Gn is the dipole form factor. The 
Coulomb distortion correction to the elastic e-p cross section exhibits an e dependence 
of (1.0-2.5)% and behaves approximately like The correction is maximum near 
Q 2 = 1 GeV 2 and then begins to decrease by increasing Q 2 . The Coulomb distortion 
explains a portion of the discrepancy but is a smaller effect than TPE. At large Q 2 , 
Ge p is suppressed and the effect of the Coulomb distortion on Ge p is larger than one 
might expect. 

6.4.4 Search For Nonlinear it ies 

Due to the linearity of the reduced cross section with e in the Born approximation, 
any deviation from linearity would have to come from higher-order terms that are 
not included in the standard radiative correction procedures. At low Q 2 , such a 
nonlinearity would provide a clean signature of TPE and give information about the 
nonlinear component of TPE. To search for a deviation from linearity in the reduced 
cross section, we fit the measured cross sections to a second-order degree polynomial 
of the form: 




(6.12) 
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where P 2 is the curvature parameter and provides a simple measure of the size of the 
nonlinear term relative to the cross section at e — 0.0. The uncertainty in P 2 , 5P2, 
sets a limit on the e 2 term. Figures 16.131 and 16.141 show such fit done for the SLAC 
NE11 HS|, and E01-001 experiments at Q 2 = 2.50 and 2.64 GeV 2 , respectively. In 
addition, the value of the extracted P2 and 5P 2 are also shown. We would like to 
compare the value of 5P 2 from the two experiments at Q 2 = 2.64 GeV 2 , however, 
since SLAC NE11 did not take any measurements at this Q 2 , we had to compare to 
the closest Q 2 available which is Q 2 = 2.50 GeV 2 . For the SLAC NE11, we get P 2 = 
0.0031 with 5P 2 = ±0.12, while the E01-001 yields P 2 = 0.0154 with 5P 2 = ±0.0445 
which is clearly a much better limit on P 2 . Moreover, as e — > 0, the variation of Pq 
between the linear and quadratic fits can help estimate an upper limit on the TPE 
contribution to 5{tG 2 m ). 
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Figure 6.13: Search for nonlinearity in the SLAC NE11 a R at Q 2 = 2.50 GeV 2 . The 
solid black line is the linear fit, while the dashed red lines are quadratic fits with 
P 2 = ± 0.120. 
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Figure 6.14: Search for nonlinearity in the E01-001 or at Q 2 = 2.64 GeV 2 . The solid 
black line is the linear fit, while the dashed red lines are quadratic fits with P2 = ± 
0.0445. 

The E01-001 experiment provides lower uncertainty in rG 2 Mp than SLAC NE11 
which is the best previous Rosenbluth measurements. Similar analysis was done at 
other Q 2 values for the E0 1-001 experiment. Table IB~4l lists the results. 

In order to be fully sensitive to nonlinearity that occurs at low e, measurements 
of the reduced cross section must be made at very low e. In fact, for an improved 
measurement of nonlinearities, it is crucial that one covers the maximum possible e 
range at both the low and high e values, and to take enough e points in the linear re- 
gion to work as a baseline for deviations from linearity since most of TPE calculations 
predicted that nonlinearities are likely to appear as e — > or e — > 1. 
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Parameter 


NE11 


E01-001 


E01-001 


E01-001 




g 2 =2.50 GeV 2 


g 2 =2.64 GeV 2 


g 2 =3.20 GeV 2 


g 2 =4.10GeV 2 


P2 


0.003 


0.015 


0.013 


0.057 


SP 2 


±0.120 


±0.044 


±0.056 


±0.12 


6(tG 2 Mv ) (%) 


3.27 


0.87 


1.10 


2.19 



Table 6.4: The curvature parameter P 2 , its uncertainty 5P 2 , and the uncertainty in 
rG 2 Mp or S(rG 2 Mp ) as extracted for the SLAC NE11 experiment at Q 2 = 2.50 GeV 2 
and for the E01-001 experiment at Q 2 = 2.64, 3.20, and 4.10 GeV 2 . 



6.5 Future Experiments and TPE 
6.5.1 Introduction 

It has been pointed out previously that the results of the E0 1-001 experiment have 
clearly established and confirmed the discrepancy between the Rosenbluth separations 
and recoil polarization results. Consequently, the E01-001 results have laid down the 
foundation for new experiments aimed at measuring the size of the TPE corrections 
at Jefferson Lab. 

Based on the formalism of Guichon and Vanderhaeghen |l()8j , which provides the 
connection between the cross section and polarization observables and goes beyond 
the Born approximation, the e-p cross section and the recoil polarization results can 
be expressed in terms of three generalized form factors or amplitudes, Ge p , Gm p , and 
F 3 which are complex functions that depend on both Q 2 and e: 

p< = _ , H l-e)\G Mr? , + x 

V T CTr V / 

Pi = V(l + s)(l-s) 1 -^ (l + ^ V 27 ) , (6.14) 
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a R 1 , e tj2 , n J-, , R 



l + -R 2 + 2e 1 + - )Y 2y , (6.15) 
\G Mp \ 2 r V t) 

where <tr is the reduced e-p cross section, Y 2l is a dimensionless variable replacing F 3 
as it will be discussed below, and R is given by: 

R=^l. (6.16) 
I Gmp I 

Note that in the Born approximation, these amplitudes are real and depend only 
on Q 2 or: G Ep — ► Ge p , Gm p — » Gm p , and F 3 — > zero. In order to separate the Born 
and higher order corrections, the generalized form factors are broken into the Born 
and TPE contributions or: 

G Ep (e, Q 2 ) = G Ep (Q 2 ) + AG Ep (e, Q 2 ) , (6.17) 



G Mp (e, Q 2 ) = G Mp {Q 2 ) + AG Mp (e, Q 2 ) , (6.18) 

F 3 (e,Q 2 ) = F 3 (e,Q 2 ) . (6.19) 
Therefore and following |lU8j . we can replace F 3 with a dimensionless variable 

Y 2l = Rer-^l), (6.20) 



^ M*\G Mp \ > 

where v = M 2 ^/(l + e)/(l - s)^t{1 + r). That way we have the two real Born 
amplitudes or G Ep and Gm p , and three two-photon amplitudes or AG Ep , AGm p , and 
Y 2t Note that AG Ep and AGm p are complex and expected to be small compared to 
the Born amplitudes. They interfere with the Born amplitudes, which are real and 
much larger, so the imaginary part of two-photon amplitudes should be negligible. 
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Therefore, the imaginary part of the two- photon amplitudes will be dropped and Ge p 
and Gm p will refer to the real part of the amplitudes. 

In the generalized formalism, and if the TPE contributions are e- independent, the 
measured ratio of the form factors can be written in terms of these generalized form 
factors for the recoil polarization as: 



From equations (|6.21j) and (J6.22)) we see that only the Y27 term leads to a difference 
between the Rosenbluth separations and recoil polarization form factor ratio and such 
difference is proportional to Y 2l . 

6.5.2 Experiment E04-019 

Experiment E04-019 |116j will measure the e dependence of the recoil polarization 
extractions of -7^-. The recoil polarization will be measured using the big Focal 



Plane Polarimeter (FPP) in Hall C and it will be done at a fixed Q 2 value of 3.20 
GeV 2 . The value Q 2 = 3.20 GeV 2 is chosen to overlap the high-precision cross section 
measurements of the Super-Rosenbluth E01-001 experiment. Measurements will be 
taken at 4 e points and the kinematics are chosen to be the same as those in the 
Super-Rosenbluth experiment so that the high-precision cross sections measured by 
the Super-Rosenbluth experiment can be used in the analysis. Based on equations 
(|6.13|) . (|6.14|) . and (|6.15|) and by using the measured values of J^, and or at a given e 
point, one can solve for Gm p , Y 2j , and R. Note that because the E04-019 experiment 




(6.21) 



and for the Rosenbluth separations (L-T) as: 




(6.22) 
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will measure the e dependence and not the size of Y 2l , it will not by itself provide 
sufficient information to correct the recoil polarization for TPE effects. Experiment 
E04-019 is scheduled to run along with the coming Hall C experiment G P E — III of 
JLAB in 2006. 



6.5.3 Experiment E05-017 

Experiment E05-017 |117j is a high-precision Rosenbluth separation measurement 
similar to the Super-Rosenbluth E01-001 experiment where protons will be detected 
with an improved precision. It aims to mapping the e and Q 2 dependence of TPE 
contributions to the elastic e-p cross section by providing high-precision Rosenbluth 
separations for 1< Q 2 <6 GeV 2 by taking more e and Q 2 points over a wide range. 
Such high precision measurements will allow for precise extraction of TPE effects from 
the difference between Rosenbluth and recoil polarization measurements. Moreover, 
it will make improved measurements of any nonlinearities in the e dependence of the 
cross sections than those reported in this work. By providing improved Rosenbluth 
separations measurements of , E05-017 will determine the Q 2 dependence of 
TPE. The overall goal is to increase the precision on / ^ Ggp by a factor of 2-3 over 
the entire Q 2 range compared to a global Rosenbluth extractions which will allow 
for precise comparison with the recoil polarization results. Experiment E05-017 is 
approved to run in Hall C of JLAB. 



6.5.4 Experiment E04-116 

Experiment E04-116 |118j we will measure the e dependence of the charge asymme- 
try R = f(f=T) where cr(e + ) and cr(e~) are the elastic cross sections of positron- and 
electron-proton scattering, respectively. Such measurements will allow for direct de- 
termination of TPE contribution to the elastic e-p scattering cross section and will 
be done for 0.5< Q 2 <2.0 and 0.1< e <0.9. 
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The crucial point here is that the radiative corrections can be classified into two 
categories: those whose their sign depends on the charge of the incident lepton, and 
those whose their sign does not depend on the charge of the incident lepton. The 
standard radiative corrections are all independent of the sign of the lepton except 
for the interference term between the electron(positron) bremsstrahlung and proton 
bremsstrahlung that is regularized by the infrared term from TPE. In other words, 
the standard radiative corrections developed by Mo and Tsai [3H| include only the 
infrared divergent contributions from TPE diagrams, which necessarily cancel the 
infrared divergent contributions from the interference term between electron and pro- 
ton bremsstrahlung, and neglect contributions from multiple soft photons exchange 
(Coulomb corrections). 

The relation of R to TPE and hence to the elastic e-p scattering cross section 
is given below: The amplitude of the elastic e-p scattering up to a 2 including only 
amplitudes that contribute to the charge asymmetry can be written as: 

A e p^ep &e&pA-Born ~\~ & e ^p-^e—bremsst &e&pAp—bremsst -\- C £ CpA2^f • (6.23) 

where e e (e p ) is the charge of the electron (proton) and the amplitudes Asorn, A e _bremsst, 
Ap-bremsst, and A 2l are the Born(single-photon exchange), electron-bremsstrahlung, 
proton-bremsstrahlung, and two-photon exchange amplitudes, respectively. By squar- 
ing equation (|6.23|) and keeping only corrections up to a that have odd powers of 
electron charge, we get: 

|^4ep-^ep| = 6 e 6p ^| AB orn | -|- 26 e 6p^4^; orn i?e(j42^,) ~\~ 2,C e 6pR6^A e —bremsstAp_j )rernss f} S j , 

(6.24) 

where Re stands for the real part of the amplitude. 

The term 2e e e p Re(A e -b remss tA* bremsst ) in equation ()6.24|) describes the interfer- 
ence between the electron and proton bremsstrahlung and its infrared divergence 
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contribution cancels exactly the corresponding infrared divergence contribution from 
2e e e p A Born Re(A2 1 ) . Therefore, after correcting the elastic e-p scattering cross section 
for this interference term, the only radiative correction that leads to a charge asym- 
metry would have to come from TPE. The effects of TPE and Coulomb corrections 
have the opposite sign for electrons and positrons or = o"s orn (l =F<5 27 ) where 5 27 

is the TPE correction yielding a charge asymmetry R = ~ 1 — 25 27 - By plotting 
R as a function of e, any deviation of R from 1.0 is a model-independent measure of 
TPE in elastic e-p scattering. Experiment E04-116 will run in Hall B of JLAB. It is 
approved for engineering run to check backgrounds with secondary beam. 

A similar e + -e~ comparison is proposed to run at VEPP-3. The experiment will 
measure the ratio R = at low e and moderate Q 2 , where no precision data exist. 
For more details, the experiment is described in reference |119j . 



Chapter 7 Summary and Conclusion 



High precision measurements of the elastic e-p scattering cross sections were made 
at the Hall A of the Thomas Jefferson National Accelerator Facility. The Super- 
Rosenbluth experiment EO 1-001 ran in May 2002 where an unpolarized electron beam 
from the Continuous Electron Beam Accelerator Facility in the range of 1.912< E < 
4.702 GeV was used and directed on a 4-cm-long unpolarized liquid-hydrogen target. 
Protons were detected, in contrast to previous measurements where the scattered 
electrons were detected, to dramatically decrease any e dependence systematic un- 
certainties and corrections. The left arm spectrometer was used to measure three 
Q 2 points of 2.64, 3.20, and 4.10 GeV 2 . Simultaneously, measurements at Q 2 = 0.5 
GeV 2 were carried out using the right arm spectrometer which served as a luminosity 
monitor to remove any uncertainties due to beam charge, current, and target density 
fluctuations. A total of 12 e points (5 e points for Q 2 = 2.64 GeV 2 , 4 e points for Q 2 
= 3.20 GeV 2 , and 3 e points for Q 2 = 4.10 GeV 2 ) were measured covering an angular 
range of 12.52°< 9i <38.26° for the left arm, while the right arm was at Q 2 = 0.5 
GeV 2 , and used to simultaneously measure 5 e points covering an angular range of 
58.29°< 9r <64.98°, allowing for high precision L-T separation of the proton electric 
and magnetic form factors and hence their ratio — ■ The measured cross sections 
were corrected for internal and external radiative corrections using updated proce- 
dure based on Mo and Tsai |3H1 ECU HH] • The absolute uncertainty (scale uncertainty) 
in the measured cross sections is approximately 3% for both arms and with relative 
uncertainties, random and slope, below 1% for the left arm, and below 1% random 
and 6% slope for the right arm. The extracted form factors, and ^ , were 
determined for the left arm at the three Q 2 points to 4%-7% and 1.5%, respectively. 
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Moreover, the ratio was determined at the three Q 2 points to 4%-7%, and found 

to approximate form factor scaling or M ^P Ep ~ 1.0. 

The results of this work are in agreement with the previous Rosenbluth data and 
inconsistent with the high-Q 2 recoil polarization results. The E01-001 experiment 
provided systematic uncertainties much smaller than the best previous Rosenbluth 
measurements ^S]; and comparable to those of the recoil polarization, clearly estab- 
lishing the discrepancy between the Rosenbluth separations and recoil polarization 
results. Furthermore, the high precision of the results confirmed that the discrep- 
ancy is not an experimental error in the Rosenbluth measurements or technique, 
confirmed the reliability of the elastic e-p scattering cross sections extracted from 
previous Rosenbluth separations, provided a strong test of the conventional radiative 
corrections used, and constituted a precise measurement of the discrepancy. 

A possible source for the discrepancy is the effect of a missing correction to the 
reduced cross sections due to two-photon exchange (TPE) corrections. Therefore, it 
is absolutely crucial to understand the radiative corrections and in particular higher- 
order processes such as TPE theoretically and experimentally before we can be confi- 
dent in our knowledge of the proton form factors. Consequently, and motivated by the 
results of the E01-001 experiment, the interest in the physics of TPE and Coulomb 
distortion corrections and their impact on electron scattering observables has risen 
experimentally and theoretically. In addition, the results of the E0 1-001 experiment 
have laid down the foundation for new experiments aimed at measuring the size of 
the TPE corrections. 

"I have noticed that nobody writes a book in his day, but is certain to 
declare on the morrow, that if this were to be altered it would be better, 
if this were to be added it would be more desirable, if this were to be 
moved forward it would be preferable and if this were to be omitted it 
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would be more pleasing. This saying is one of the wisest maxims and it 
is evidence that the whole human race is seized by imperfection." 
— Al-Imad Al-Asfahani (Al-Katib), 1125-1200 A.D.: A historian and one of the major 
writers. He accompanied Salah Al-Din Al-Ayyubi (Saladin) and recorded his history. 
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